Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,061,935 G→A S46N (AGC→AAC)  torD → trimethylamine‑N‑oxide reductase‑specific chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,061,9350GA90.5% 101.5 / ‑5.5 42S46N (AGC→AAC) torDtrimethylamine‑N‑oxide reductase‑specific chaperone
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/20);  minor base C (0/4);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 1.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAGCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTG  >  NC_000913/1061695‑1062176
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
gAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAAccgccg                                                                                                                                                                                                                                         <  1:289492/251‑1 (MQ=255)
      gagaatgtaaacggaacggTGGAGCAGGTGACGGCTTTTAACGGCCCCGGGGAGTTGGTGGCGCAGTGCGAATATGTCCCCGCGTCGCAGGGGAATTCATGGCCACGCTGACAGCCCAACAGTTTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGGGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCGGAATGTTTTTCTTTTCTGAAAACCGAACCGCCGCTCACt                                                                                                                                                                                                                                   <  2:243515/242‑1 (MQ=255)
                     aCAGGGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTAAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTCCGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAg                                                                                                                                                                                                                    <  1:366752/251‑1 (MQ=255)
                      caggggAGCAGGTGACGGTGTTTACCGGCCCCGTGGGGATGGTGGGGCAGTGAGAATATGTCCCCGCGTCGCAGGTGAAATCATGACCACGCTGAAAGCAAAACGGATTGCGTGTTTTTACGCTTGGCTAGCGCATTTTGTCCCCCGTGAGCTGAACGAGGAACAACGGACGCAAACCGGAAGTGCGCAGATGGGTGAATGGGTTGGGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGc                                                                                                                                                                                                                   <  1:33049/247‑1 (MQ=255)
                            agccggTGGCGGCGTTTAACGGCCCCGTGGAGAGGGTGGCGCGGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCCCGCGGACAGCACAACAGATTGCTTGTGTTTACGCCTGTCTAGCGCAGTTTTTCACCCGGGGGCTGGCCGATGAAAAACGGACGCAAATCACCAGGGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGaaa                                                                                                                                                                                                             <  1:321706/247‑1 (MQ=255)
                              ggggTGAGGGGGTTTAACGGCCCCGTGGAGAGGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGGAAAATCATGCCCCCGCTGACAGCACAACAGTTTGCCGGTGTTTACGCCTGGCTAGCGCATTTGTTCTCCCGTGAGCGGGACGATGAACAACTGACGCAAATCGCCGGTGCGCAGATGGCTGAATGGTTTTGGTTGCTGAAAACCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAAc                                                                                                                                                                                                           <  2:134225/249‑1 (MQ=255)
                                     cGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGGGAACGAGCTGGAAAACCGTATTg                                                                                                                                                                                                    >  2:52653/1‑251 (MQ=255)
                                                     ggTGGAGATGGTGGCGAAGTGCGAATATGTTCCCGCGTCGCAGGGGAAATCAGGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCATTTGTTCTCCCGTGAGCTGGACGATAAACAACTGACGCAAATCGCCAGGGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACg                                                                                                                                                                                    <  1:383331/250‑1 (MQ=255)
                                                                   ggcggtgagAATATGTCCCCGTGTCGCAGGTAAAATCATGACCACGCTGAAAGCACACAAGTTGGCTTGTGTTTACGCCGGGTGAGCTCTGTTGTTCTCCGGGGAGCGGGAGGATAAACAACTGACGAAAATCGCCAGTGCGCAGTTGTCTGAGTGTTTTTGGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTc                                                                                                                                                                      <  2:271372/243‑1 (MQ=255)
                                                                    gcAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGCCCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCt                                                                                                                                                                     <  1:52653/251‑1 (MQ=255)
                                                                          cGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGGATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAAcc                                                                                                                                                               >  1:225719/1‑250 (MQ=255)
                                                                                           cGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAACTCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACACTGAAAGTACGTGACGATGCCCGTCTGGGAATGGCcgagggccttttcg                                                                                                                                              >  2:370836/1‑239 (MQ=255)
                                                                                            gCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGGGAACGAGCTGGCAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCgg                                                                                                                                             >  1:244937/1‑251 (MQ=255)
                                                                                                        aTGACCACGCTGACAGCACACCAGATTGCCTGTGTTTACGCCTGGCTAGCGCATTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCtgat                                                                                                                                 <  1:146170/251‑1 (MQ=255)
                                                                                                        aTGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGCTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCAATGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGAGCGTACGTGACGAGGCACGTCGGGCACTGGCCGCGGACTTTTGCGGCCTGtttttgcg                                                                                                                                 >  2:391903/1‑246 (MQ=255)
                                                                                                          gACCACGTGGACAGCACAACAGAGTGCTGTGGTTTCCGCCTGGCTAGCGCAGTTGTTCCCCCGTGAGCGGGACGATGAACAACTGACGCAAATCGCCGGTGCGCAGAGGGTGGAATGGTTTTCTTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCtgatg                                                                                                                                <  1:398159/250‑1 (MQ=255)
                                                                                                                  tGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACCGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCggcca                                                                                                                        >  2:28324/1‑246 (MQ=255)
                                                                                                                    aCAGCACAACAGATTGCATGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACACCTGACACAAATCCCCAGTGCGCAGAAGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGACACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGAACGacaacc                                                                                                                      >  1:31617/1‑248 (MQ=255)
                                                                                                                         ccaactGATGGCCTGTGTTCACGCCTGGGTAGCGCGTTTGTCCTCCCGTGGGCGGGACGATGAACAACGAACGAAAATGCCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAgcgg                                                                                                                <  2:360651/250‑1 (MQ=255)
                                                                                                                          caacaGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAgcggc                                                                                                               >  2:330779/1‑251 (MQ=255)
                                                                                                                           caacacatggctggtttgtgcGCCTGGCTGGCGCAGTTGTTCCCCCGGGGGCTGGACGAGGAAAAACTGACGCAAATCGCCGGTGGGCAGATGGTGGAAGGGTTTTCGTTGGTAAAAACCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACGTGTGCGGCCTGTTTCTGATGACCGACAAACAAgcggc                                                                                                               <  2:406134/231‑1 (MQ=255)
                                                                                                                                atttctgttgtTTACGCCTGTCTAGCGCATTTGTTCTCCCGTGAGCTGGACGAGGAACAACTGAGGCAAATCGCCGGTGCGCGGAGGGCTGAATGGTTTGCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGcc                                                                                                         <  1:165177/243‑1 (MQ=255)
                                                                                                                                         gtTTCCGCCTGTCTAGCGCAGTTGTTCTCCCGTGGGCTGGCCAAGGAACAACGGACGCAAATCGCCAGTGCGCAGAGGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATc                                                                                                <  1:176948/251‑1 (MQ=255)
                                                                                                                                            ttACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTAttcctcgg                                                                                              >  2:333645/1‑244 (MQ=255)
                                                                                                                                             tACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGc                                                                                             >  2:76816/1‑251 (MQ=255)
                                                                                                                                              aCGCCTGGCTAGGGCAGTTTTTCTCCCGTGGGCGGGATGATGAACAACTGCCGAAAATCGCAAGTGCGCAGTTGGCTAAAGTGTTTTTGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGcc                                                                                            <  1:262150/251‑1 (MQ=255)
                                                                                                                                                  ctagctgacgcAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGAAGCGAATCGCCAGTGCGCAGATGGCTGACTGGGTTTCGTTGCTTAACAACGAACCGCCGCTCACTGCGGCGGTATACGGGATGGAGAACAGTCGTGGCACACTGACAGTTCTTGACGAAGCCCCGCTGGACATGGCCGCGGGGATTTGGGGCCTATTTCTGCGGACCGACAAAAAAGGGGACCTGCCGTATGCATCCGccctc                                                                                         >  2:703/9‑247 (MQ=255)
                                                                                                                                                   tGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACttt                                                                                                                                                  >  2:15190/1‑192 (MQ=255)
                                                                                                                                                   tGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGAGGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACttt                                                                                                                                                  <  1:15190/192‑1 (MQ=255)
                                                                                                                                                                tgttCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCa                                                                           >  1:227194/1‑250 (MQ=255)
                                                                                                                                                                          ggAGCTGGACGATAAACAACTTACAAAAATCGCCGGTGCGCAGATGGGGAAATGGTTTTGGTTGATGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAAc                                                                <  2:348079/250‑1 (MQ=255)
                                                                                                                                                                          ggAGCGGGACGATGAACAACTGACGAAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAAc                                                                <  1:142688/250‑1 (MQ=255)
                                                                                                                                                                             gtgggtgtagaaaaaacgaacgaaaaTCGCATGTGCGCATGTGGCTGAATGTTTTTGGTGGCTGAAAACCACACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGTCGAGCAAGAGAGTAAACGCt                                                             <  2:240917/228‑1 (MQ=255)
                                                                                                                                                                               gggacgaagaacagctaccgaaaaTCCCCGGTGCGCAGAGGGCTGAATGGTTTTTTTTGCGGAAAAACGACCCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTg                                                           <  1:53950/230‑1 (MQ=255)
                                                                                                                                                                                ggACGTTGAACACCGGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTGCTTTGCTAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGTGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCtt                                                            <  1:276630/249‑1 (MQ=255)
                                                                                                                                                                                                     aaTCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCCGGGATGGaa                                      >  2:291720/1‑250 (MQ=255)
                                                                                                                                                                                                               gcgcAGATGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAagag                                                                       >  1:35026/1‑207 (MQ=255)
                                                                                                                                                                                                                     aTGGCTGAATGTTTTTCGTTGCGAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCaa                     <  1:52917/251‑1 (MQ=255)
                                                                                                                                                                                                                     aTGGCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAAACAGCGGCAAttttc                      >  2:52917/1‑248 (MQ=255)
                                                                                                                                                                                                                       ggCTGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACg                   <  1:45145/251‑1 (MQ=255)
                                                                                                                                                                                                                          tGAATGGTTTTCGTTGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAAc                >  1:389104/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tGCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCCCGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACCAACAAGCGGCGCCGCCGTATGCCTCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTTAGGCAGGGGGGGAAACCAGCGGCAATTTTAACGaaaccgcaggacctcct  >  2:23962/1‑236 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
GAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTTGCTGAAAAGCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACGTGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTG  >  NC_000913/1061695‑1062176

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: