Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,395,812 G→C 5.3% L29V (CTG→GTG)  nuoJ ← NADH:quinone oxidoreductase subunit J

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,395,8120GC 5.3% 260.1 / 4.2 76L29V (CTG→GTG) nuoJNADH:quinone oxidoreductase subunit J
Reads supporting (aligned to +/- strand):  ref base G (44/28);  new base C (2/2);  total (46/30)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

ACCGCCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTTCTCTGCTTCGCC  >  NC_000913/2395663‑2395961
                                                                                                                                                      |                                                                                                                                                     
aCCGCCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTAca‑gcag                                                                                                                                                       <  1:1581808/150‑1 (MQ=255)
  cGCCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTg                                                                                                                                                     <  1:1442323/150‑1 (MQ=255)
        ggTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTa                                                                                                                                               >  2:913249/1‑150 (MQ=255)
             aGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGa                                                                                                                                          >  2:875899/1‑150 (MQ=255)
             aGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGa                                                                                                                                          >  2:627635/1‑150 (MQ=255)
                     caccacGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATg                                                                                                                                  <  1:139115/150‑1 (MQ=255)
                      accacGAACACAAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCCGGTACACGCACTGC‑TGTACCTGAttatttcgc                                                                                                                                 <  1:411908/150‑8 (MQ=38)
                       cctcgaagacgaacaGCACCTTAATGGCACCCGCGTAGAAGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCGAGTGAGAAGAACACCGCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGGAGCGGATTGGTATggg                                                                                                                                <  2:1596846/147‑1 (MQ=255)
                       ccacGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCCGGTACACGCACTGC‑TGTACCTGAttatttcgc                                                                                                                                 >  1:1022676/1‑142 (MQ=38)
                       ccacGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCCGGTACACGCACTGC‑TGTACCTGAttatttcgc                                                                                                                                 <  2:1022676/149‑8 (MQ=38)
                       ccacGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATggg                                                                                                                                <  2:85766/150‑1 (MQ=255)
                       acacgaaaacgaaAAGCACCATAAAGGCACACGCGTACACTATAATATCAAGCGCACCTACGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATggg                                                                                                                                >  2:563313/2‑150 (MQ=255)
                                      gCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAgg                                                                                                                 >  1:1426274/1‑150 (MQ=255)
                                      gCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAgg                                                                                                                 >  2:636414/1‑150 (MQ=255)
                                      gCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAgg                                                                                                                 >  1:569585/1‑150 (MQ=255)
                                      gCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAgg                                                                                                                 <  2:1462773/150‑1 (MQ=255)
                                       cACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTAGA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGt                                                                                                                >  1:745488/1‑150 (MQ=255)
                                         ccATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCg                                                                                                              >  1:371447/1‑150 (MQ=255)
                                            tAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCaa                                                                                                           <  2:284946/150‑1 (MQ=255)
                                            tAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCaa                                                                                                           <  1:1422776/150‑1 (MQ=255)
                                              aTGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGt                                                                                                         >  2:61435/1‑150 (MQ=255)
                                                ggCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTAt                                                                                                       <  2:371447/150‑1 (MQ=255)
                                                  cACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATgg                                                                                                     >  2:592844/1‑150 (MQ=255)
                                                  cACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATgg                                                                                                     >  1:896231/1‑150 (MQ=255)
                                                      cgcgTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTAt                                                                                                 <  1:913249/150‑1 (MQ=255)
                                                            gACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGcc                                                                                           <  1:1250403/150‑1 (MQ=255)
                                                              cGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCGAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCcac                                                                                         >  2:977506/1‑150 (MQ=255)
                                                              cGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCcac                                                                                         >  1:738872/1‑150 (MQ=255)
                                                              cGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCcac                                                                                         >  1:1620745/1‑150 (MQ=255)
                                                                   aTTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGata                                                                                    <  2:569585/150‑1 (MQ=255)
                                                                         aGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATaaaaa                                                                              >  2:782002/1‑150 (MQ=255)
                                                                          gcgcACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAg                                                                             >  2:1312040/1‑150 (MQ=255)
                                                                             cACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGa                                                                          <  2:1356605/150‑1 (MQ=255)
                                                                               ccGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAAc                                                                        >  2:831176/1‑150 (MQ=255)
                                                                                cGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GAAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACt                                                                       >  1:1302292/1‑150 (MQ=255)
                                                                                 ggCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTc                                                                      >  1:936181/1‑150 (MQ=255)
                                                                                 ggCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGACTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTc                                                                      >  2:882582/1‑150 (MQ=255)
                                                                                  gCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTcc                                                                     <  1:1048122/150‑1 (MQ=255)
                                                                                      aGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTg                                                                 <  1:719691/150‑1 (MQ=255)
                                                                                       gTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGc                                                                <  1:366676/150‑1 (MQ=255)
                                                                                         aaGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGccc                                                              <  2:792843/150‑1 (MQ=255)
                                                                                           gcgcCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGAGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCct                                                            >  2:1523230/1‑150 (MQ=255)
                                                                                            cgcCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCctc                                                           >  1:392969/1‑150 (MQ=255)
                                                                                                cAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCt                                                       >  2:1179519/1‑150 (MQ=255)
                                                                                                    gagaAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACg                                                   >  2:85127/1‑150 (MQ=255)
                                                                                                     agaagaACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACgg                                                  <  2:19046/150‑1 (MQ=255)
                                                                                                      gaagaaCACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGt                                                 >  1:1439047/1‑150 (MQ=255)
                                                                                                      gaagaaCACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGt                                                 >  1:1444937/1‑150 (MQ=255)
                                                                                                         gaaCACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAAc                                              >  2:1240632/1‑150 (MQ=255)
                                                                                                              ccccGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGct                                         >  1:440285/1‑150 (MQ=255)
                                                                                                               cccGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGctc                                        >  1:18835/1‑150 (MQ=255)
                                                                                                               cccGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGctc                                        >  2:424303/1‑150 (MQ=255)
                                                                                                                ccGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGctct                                       <  2:1537912/150‑1 (MQ=255)
                                                                                                                ccGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGctct                                       <  1:1312040/150‑1 (MQ=255)
                                                                                                                ccGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGctct                                       >  2:265652/1‑150 (MQ=255)
                                                                                                                     aaTCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACg                                  <  1:763922/150‑1 (MQ=255)
                                                                                                                     aaTCCCCAGCAGCTAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACg                                  >  2:93262/1‑150 (MQ=255)
                                                                                                                      aTCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGt                                 >  2:1567019/1‑150 (MQ=255)
                                                                                                                      aTCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGt                                 <  1:57608/150‑1 (MQ=255)
                                                                                                                           cagcagCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATa                            >  1:880800/1‑150 (MQ=255)
                                                                                                                           cagcagCGAAATAATCAGGTACA‑GCAGTGAGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATa                            <  1:1523230/150‑1 (MQ=255)
                                                                                                                             gcagcGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAgg                          <  1:563313/150‑1 (MQ=255)
                                                                                                                              cagcGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggc                         >  1:620683/1‑150 (MQ=255)
                                                                                                                                gcGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggctt                       >  1:887416/1‑150 (MQ=255)
                                                                                                                                gcGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggctt                       >  1:630308/1‑150 (MQ=255)
                                                                                                                                gcGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggctt                       >  2:154534/1‑150 (MQ=255)
                                                                                                                                  gAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggcttgg                     <  1:85127/150‑1 (MQ=255)
                                                                                                                                     ccaatCCGGTACACGCACTGC‑TGTACCTGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggcttggctt                  >  1:1606218/3‑150 (MQ=255)
                                                                                                                                     ataatCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggcttggctt                  <  2:392969/150‑1 (MQ=255)
                                                                                                                                      taatCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTc                 >  1:479585/1‑150 (MQ=255)
                                                                                                                                       aatCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCg                >  1:613299/1‑150 (MQ=255)
                                                                                                                                        atCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGt               >  2:347187/1‑150 (MQ=255)
                                                                                                                                        atCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGt               <  1:1240632/150‑1 (MQ=255)
                                                                                                                                         tCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGtt              <  2:1620745/150‑1 (MQ=255)
                                                                                                                                            ggTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTtctc           >  2:277206/1‑150 (MQ=255)
                                                                                                                                              tACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTtctact         >  1:976430/1‑147 (MQ=255)
                                                                                                                                              tACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTTCTCTg         <  1:93262/150‑1 (MQ=255)
                                                                                                                                                    ccgTGCTTGTACCGTATTGTTAGGGGGGAGACCTCGCAAGGCCGCACGAAGGGCTATCCGGCCAACGATAT‑AAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTTCTCTGCTTCGc   <  1:875899/148‑1 (MQ=255)
                                                                                                                                                      gTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTTCTCTGCTTCGcc  <  1:627635/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
ACCGCCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACA‑GCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTTCTCTGCTTCGCC  >  NC_000913/2395663‑2395961

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑