Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,395,808:1 +C 5.3% coding (89/555 nt) 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,8081.C 5.3% 319.9 / 12.1 76coding (89/555 nt)nuoJNADH:quinone oxidoreductase subunit J
Reads supporting (aligned to +/- strand):  ref base . (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.31e-01

CTGAACCGCCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTTCTCTG  >  NC_000913/2395659‑2395954
                                                                                                                                                      |                                                                                                                                                  
cTGAACCGCCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTAca                                                                                                                                                     <  1:1202938/150‑1 (MQ=255)
 tGAACCGCCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTAcag                                                                                                                                                   >  2:763922/1‑150 (MQ=255)
  gAACCGCCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTAcagc                                                                                                                                                  <  2:387063/150‑1 (MQ=255)
   aaCCGCCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATATCAGGTAcagccg                                                                                                                                                >  2:1587811/1‑148 (MQ=255)
    aCCGCCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTAcagcag                                                                                                                                                <  1:1581808/150‑1 (MQ=255)
      cGCCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTg                                                                                                                                              <  1:1442323/150‑1 (MQ=255)
            ggTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTa                                                                                                                                        >  2:913249/1‑150 (MQ=255)
                 aGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGa                                                                                                                                   >  2:875899/1‑150 (MQ=255)
                 aGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGa                                                                                                                                   >  2:627635/1‑150 (MQ=255)
                         caccacGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATg                                                                                                                           <  1:139115/150‑1 (MQ=255)
                          accacGAACACAAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCCGGTACACGCACTGCTGTACCTGAttatttcgc                                                                                                                          <  1:411908/150‑8 (MQ=38)
                           cctcgaagacgaacaGCACCTTAATGGCACCCGCGTAGAAGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCGAGTGAGAAGAACACCGCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGGAGCGGATTGGTATggg                                                                                                                         <  2:1596846/147‑1 (MQ=255)
                           ccacGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCCGGTACACGCACTGCTGTACCTGAttatttcgc                                                                                                                          >  1:1022676/1‑142 (MQ=38)
                           ccacGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCCGGTACACGCACTGCTGTACCTGAttatttcgc                                                                                                                          <  2:1022676/149‑8 (MQ=38)
                           ccacGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATggg                                                                                                                         <  2:85766/150‑1 (MQ=255)
                           acacgaaaacgaaAAGCACCATAAAGGCACACGCGTACACTATAATATCAAGCGCACCTACGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATggg                                                                                                                         >  2:563313/2‑150 (MQ=255)
                                          gCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAgg                                                                                                          >  2:636414/1‑150 (MQ=255)
                                          gCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAgg                                                                                                          >  1:1426274/1‑150 (MQ=255)
                                          gCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAgg                                                                                                          <  2:1462773/150‑1 (MQ=255)
                                          gCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAgg                                                                                                          >  1:569585/1‑150 (MQ=255)
                                           cACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTAGAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGt                                                                                                         >  1:745488/1‑150 (MQ=255)
                                             ccATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCg                                                                                                       >  1:371447/1‑150 (MQ=255)
                                                tAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCaa                                                                                                    <  2:284946/150‑1 (MQ=255)
                                                tAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCaa                                                                                                    <  1:1422776/150‑1 (MQ=255)
                                                  aTGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGt                                                                                                  >  2:61435/1‑150 (MQ=255)
                                                    ggCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTAt                                                                                                <  2:371447/150‑1 (MQ=255)
                                                      cACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATgg                                                                                              >  2:592844/1‑150 (MQ=255)
                                                      cACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATgg                                                                                              >  1:896231/1‑150 (MQ=255)
                                                          cgcgTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTAt                                                                                          <  1:913249/150‑1 (MQ=255)
                                                                gACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGcc                                                                                    <  1:1250403/150‑1 (MQ=255)
                                                                  cGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCGAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCcac                                                                                  >  2:977506/1‑150 (MQ=255)
                                                                  cGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCcac                                                                                  >  1:1620745/1‑150 (MQ=255)
                                                                  cGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCcac                                                                                  >  1:738872/1‑150 (MQ=255)
                                                                       aTTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGata                                                                             <  2:569585/150‑1 (MQ=255)
                                                                             aGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATaaaaa                                                                       >  2:782002/1‑150 (MQ=255)
                                                                              gcgcACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAg                                                                      >  2:1312040/1‑150 (MQ=255)
                                                                                 cACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGa                                                                   <  2:1356605/150‑1 (MQ=255)
                                                                                   ccGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAAc                                                                 >  2:831176/1‑150 (MQ=255)
                                                                                    cGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGAAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACt                                                                >  1:1302292/1‑150 (MQ=255)
                                                                                     ggCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTc                                                               >  1:936181/1‑150 (MQ=255)
                                                                                     ggCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGACTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTc                                                               >  2:882582/1‑150 (MQ=255)
                                                                                      gCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTcc                                                              <  1:1048122/150‑1 (MQ=255)
                                                                                          aGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTg                                                          <  1:719691/150‑1 (MQ=255)
                                                                                           gTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGc                                                         <  1:366676/150‑1 (MQ=255)
                                                                                             aaGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGccc                                                       <  2:792843/150‑1 (MQ=255)
                                                                                               gcgcCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGAGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCct                                                     >  2:1523230/1‑150 (MQ=255)
                                                                                                cgcCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCctc                                                    >  1:392969/1‑150 (MQ=255)
                                                                                                    cAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCt                                                >  2:1179519/1‑150 (MQ=255)
                                                                                                        gagaAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACg                                            >  2:85127/1‑150 (MQ=255)
                                                                                                         agaagaACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACgg                                           <  2:19046/150‑1 (MQ=255)
                                                                                                          gaagaaCACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGt                                          >  1:1439047/1‑150 (MQ=255)
                                                                                                          gaagaaCACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGt                                          >  1:1444937/1‑150 (MQ=255)
                                                                                                             gaaCACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAAc                                       >  2:1240632/1‑150 (MQ=255)
                                                                                                                  ccccGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGct                                  >  1:440285/1‑150 (MQ=255)
                                                                                                                   cccGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGctc                                 >  1:18835/1‑150 (MQ=255)
                                                                                                                   cccGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGctc                                 >  2:424303/1‑150 (MQ=255)
                                                                                                                    ccGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGctct                                <  2:1537912/150‑1 (MQ=255)
                                                                                                                    ccGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGctct                                <  1:1312040/150‑1 (MQ=255)
                                                                                                                    ccGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGctct                                >  2:265652/1‑150 (MQ=255)
                                                                                                                         aaTCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACg                           <  1:763922/150‑1 (MQ=255)
                                                                                                                         aaTCCCCAGCAGCTAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACg                           >  2:93262/1‑150 (MQ=255)
                                                                                                                          aTCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGt                          <  1:57608/150‑1 (MQ=255)
                                                                                                                          aTCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGt                          >  2:1567019/1‑150 (MQ=255)
                                                                                                                               cagcagCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATa                     >  1:880800/1‑150 (MQ=255)
                                                                                                                               cagcagCGAAATAATCAGGTACAGCAGTGAGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATa                     <  1:1523230/150‑1 (MQ=255)
                                                                                                                                 gcagcGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAgg                   <  1:563313/150‑1 (MQ=255)
                                                                                                                                  cagcGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggc                  >  1:620683/1‑150 (MQ=255)
                                                                                                                                    gcGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggctt                >  1:887416/1‑150 (MQ=255)
                                                                                                                                    gcGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggctt                >  2:154534/1‑150 (MQ=255)
                                                                                                                                    gcGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggctt                >  1:630308/1‑150 (MQ=255)
                                                                                                                                      gAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggcttgg              <  1:85127/150‑1 (MQ=255)
                                                                                                                                         ccaatCCGGTACACGCACTGCTGTACCTGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggcttggctt           >  1:1606218/3‑150 (MQ=255)
                                                                                                                                         ataatCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAggcttggctt           <  2:392969/150‑1 (MQ=255)
                                                                                                                                          taatCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTc          >  1:479585/1‑150 (MQ=255)
                                                                                                                                           aatCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCg         >  1:613299/1‑150 (MQ=255)
                                                                                                                                            atCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGt        >  2:347187/1‑150 (MQ=255)
                                                                                                                                            atCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGt        <  1:1240632/150‑1 (MQ=255)
                                                                                                                                             tCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGtt       <  2:1620745/150‑1 (MQ=255)
                                                                                                                                                ggTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTtctc    >  2:277206/1‑150 (MQ=255)
                                                                                                                                                  tACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTtctact  >  1:976430/1‑147 (MQ=255)
                                                                                                                                                  tACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTTCTCTg  <  1:93262/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CTGAACCGCCCAGGTTGAGCATCATCACCACGAACACGAACAGCACCATAATGGCACCCGCGTAGACGATAATTTCCAGCGCACCGGCGAAGTAAGCGCCCAGTGAGAAGAACACCCCGGAAATCGCCAGCAGCGAAATAATCAGGTACAGCAGTGCGTGTACCGGATTGGTATGGGTGATCACTCGCAAGGTCGCAAGTATGGCTATCAGGCCACAGATATAAAAAGCGAACTCCATTGCCCCTCTCCTTACGGTAACAGGCTCTTGACGTCGATAGGCTTGGCTTCGTTCTCTG  >  NC_000913/2395659‑2395954

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: