Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,668,461 G→A 5.5% G96D (GGC→GAC)  nuoM → NADH‑quinone oxidoreductase subunit M

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,668,4610GA 5.5% 212.5 / 2.1 72G96D (GGC→GAC) nuoMNADH‑quinone oxidoreductase subunit M
Reads supporting (aligned to +/- strand):  ref base G (38/30);  new base A (2/2);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

TGGGGTACCGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGACCTGTTCCTGTTCTTCTTCTTCTGGGA  >  NC_002947/4668313‑4668608
                                                                                                                                                    |                                                                                                                                                    
tggGGTACCGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCCGc                                                                                                                                                     <  2:379808/150‑1 (MQ=255)
    gTACCGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGc                                                                                                                                                 >  1:554187/1‑150 (MQ=255)
        cGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCTAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCgg                                                                                                                                             >  2:1310987/1‑150 (MQ=255)
         ggCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGgt                                                                                                                                            <  2:462462/150‑1 (MQ=255)
            gACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGtgtg                                                                                                                                         >  2:353203/1‑150 (MQ=255)
            gACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGtgtg                                                                                                                                         <  2:255315/150‑1 (MQ=255)
               tACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCAGCTCACCGGCCTGCTCGGTGTGCTg                                                                                                                                      >  2:1304071/1‑150 (MQ=255)
                          tCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCTAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACtgtg                                                                                                                           <  1:1310987/150‑1 (MQ=255)
                          tCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACtgtg                                                                                                                           >  2:1276703/1‑150 (MQ=255)
                          tCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACtgtg                                                                                                                           >  1:870631/1‑150 (MQ=255)
                            ccGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGgtgt                                                                                                                                          <  2:1351707/133‑1 (MQ=255)
                            ccGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGgtgt                                                                                                                                          >  1:1351707/1‑133 (MQ=255)
                             cGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTc                                                                                                                        <  1:326342/150‑1 (MQ=255)
                                cccGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTg                                                                                                                     >  2:607089/1‑150 (MQ=255)
                                 ccGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCCGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTgg                                                                                                                    >  2:243430/1‑150 (MQ=255)
                                  cGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGa                                                                                                                   >  1:237514/1‑150 (MQ=255)
                                        ccGAACCGGCCTGGGCTGTTGAAGACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCATCTCCAGCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAgaga                                                                                                             <  1:799104/150‑1 (MQ=255)
                                         cGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATGCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGAt                                                                                                            >  1:1358418/1‑150 (MQ=255)
                                           aaCCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATcc                                                                                                          <  2:239243/150‑1 (MQ=255)
                                             ccGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAg                                                                                                        <  1:773051/150‑1 (MQ=255)
                                                   tGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCAc                                                                                                  >  1:729391/1‑150 (MQ=255)
                                                       ctctCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCCGGACGGCCTGTCGCTGCTGATGATCCTGCTCCCCGGCCTGCTAGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCg                                                                                              <  1:243430/150‑1 (MQ=255)
                                                          tcGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGct                                                                                           <  2:870631/150‑1 (MQ=255)
                                                             aaTACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCAGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGct‑tct                                                                                       <  1:1304071/150‑1 (MQ=255)
                                                               tACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCt‑tcttc                                                                                     >  1:294379/1‑150 (MQ=255)
                                                               tACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCt‑tcttc                                                                                     >  2:1065410/1‑150 (MQ=255)
                                                                aCAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTcc                                                                                    <  2:294379/150‑1 (MQ=255)
                                                                  aaaGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCAc                                                                                  >  1:533342/1‑150 (MQ=255)
                                                                             gATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGAt                                                                       >  2:889787/1‑150 (MQ=255)
                                                                               tCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGAtgc                                                                     >  1:950444/1‑149 (MQ=255)
                                                                               tCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGACCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAgaga                                                                                                             <  2:867908/111‑1 (MQ=255)
                                                                               tCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGACCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAgaga                                                                                                             >  1:867908/1‑111 (MQ=255)
                                                                                   gCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTgtgt                                                                                                                          >  2:50561/1‑94 (MQ=255)
                                                                                   gCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTgtgt                                                                                                                          <  1:50561/94‑1 (MQ=255)
                                                                                     gTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATcc                                                               >  2:1405970/1‑150 (MQ=255)
                                                                                     gTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGCTCTGTCGGTACTGTGTTCCTAGAAAAAGATCAAGCTCAACGTCTGCT‑TCTTCCACCTCAACCTGATGTGGATcc                                                               >  1:94985/1‑150 (MQ=255)
                                                                                           gCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTAGGTGTGGTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTcggcg                                                         >  2:547377/1‑150 (MQ=255)
                                                                                              tCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGgcggcg                                                      >  1:414928/1‑150 (MQ=255)
                                                                                                 gCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTgg                                                   <  2:664360/150‑1 (MQ=255)
                                                                                                  cATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGt                                                  >  2:2830/1‑150 (MQ=255)
                                                                                                      cACCTGGCCCTGGACGGCCTGTCGGTGCTGATGATCCTGCTCACTGGCCTGTTCGGTGTGCTGTAGGTACTGTGTTCCTGGAAAGAGTTCCAGCGCCACGTCGCCT‑TCTTCCACCTCAACATGATGTGGGTCCTCGGGGCCGTGGTCGgt                                              >  2:153822/1‑150 (MQ=255)
                                                                                                      cACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGgt                                              <  2:414928/150‑1 (MQ=255)
                                                                                                      cACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGgt                                              <  2:729391/150‑1 (MQ=255)
                                                                                                       cccTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGt                                                  >  1:984923/2‑145 (MQ=255)
                                                                                                       aCCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTAGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGt                                                  <  2:984923/145‑1 (MQ=255)
                                                                                                       aCCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACAGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGgtg                                             >  2:1033932/1‑150 (MQ=255)
                                                                                                        ccTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGgtgt                                            >  2:1177539/1‑150 (MQ=255)
                                                                                                           ggCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑GCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGtt                                         <  2:554187/150‑1 (MQ=255)
                                                                                                              ccTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCt                                      >  1:143430/1‑150 (MQ=255)
                                                                                                              ccTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCt                                      >  1:1325624/1‑150 (MQ=255)
                                                                                                              ccTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGATGTCGATACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCt                                      >  2:126473/1‑150 (MQ=255)
                                                                                                               cTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTg                                     <  2:1358418/150‑1 (MQ=255)
                                                                                                                     ggCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTc                                                                                               <  2:1302664/87‑1 (MQ=255)
                                                                                                                     ggCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTc                                                                                               >  1:1302664/1‑87 (MQ=255)
                                                                                                                       ccTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCa                                                                                                   >  2:437491/1‑81 (MQ=255)
                                                                                                                       ccTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCa                                                                                                   <  1:437491/81‑1 (MQ=255)
                                                                                                                       ccTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGa                             >  2:1250565/1‑150 (MQ=255)
                                                                                                                        cTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGAc                            <  1:1250565/150‑1 (MQ=255)
                                                                                                                           tCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGAcctg                         >  2:72375/1‑150 (MQ=255)
                                                                                                                              ctgctgAGGATCCTGCTCACCGACCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGAcctgttc                      >  1:836990/1‑150 (MQ=255)
                                                                                                                                      gatCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGACCTGTTCCTGttctt              <  1:72375/150‑1 (MQ=255)
                                                                                                                                       atCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGACCTGTTCCTGttcttc             <  2:237514/150‑1 (MQ=255)
                                                                                                                                           tgctactcGGTCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGCGATATCAAACCGCCT‑GGTTGTGGTTTGCCTCACCCTGATTCGAATCCTCGGCGGCGTGGTTGGTGTGTGTTTTGCCCTTGACCTGTTCCTGttcttcttct         <  1:1033932/143‑1 (MQ=255)
                                                                                                                                           tGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCGCTGGACCTGTTCATGttcttcttct         >  1:297100/1‑150 (MQ=255)
                                                                                                                                               cACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGCGGTCGGTGTGTTCCTGGCCCTGGACCTGTTCCTGATATTCTCctgctg     <  2:1348430/150‑1 (MQ=255)
                                                                                                                                                aCCGACCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGACCTGTTCCTGTTCTTCTTCTTCTgg    <  2:836990/150‑1 (MQ=255)
                                                                                                                                                 ccGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGACCTGTTCCTGTTCTTCTTCTTCTggg   <  1:607089/150‑1 (MQ=255)
                                                                                                                                                  cGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCgg                                               <  2:1124007/105‑1 (MQ=255)
                                                                                                                                                  cGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCgg                                               >  1:1124007/1‑105 (MQ=255)
                                                                                                                                                  cGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGACCTGTTCCTGTtcttcttcttct      >  2:71241/1‑146 (MQ=255)
                                                                                                                                                  cGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGACCTGTTCCTGTtcttcttcttct      <  1:71241/146‑1 (MQ=255)
                                                                                                                                                  cGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGAGCCACGTCGGCT‑TCTTCCACCTCAACCTGATCTGGACCCTCGTCGGCGGGGTCCGGGTGTTCGTGGCCGCGGGCTTGTTCCCGTTCCTCTTCTTCTgggt  >  2:1403385/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGGGGTACCGGCGACTACACCCTTGCTCCCGCCCCGGGCGCCGAACCGGCCTGGGCTCTCGAATACAAAGTCGAGTGGATCAAGCGTTTCGGCATCAGCATCCACCTGGCCCTGGACGGCCTGTCGCTGCTGATGATCCTGCTCACCGGCCTGCTCGGTGTGCTGTCGGTACTGTGTTCCTGGAAAGAGATCCAGCGCCACGTCGGCT‑TCTTCCACCTCAACCTGATGTGGATCCTCGGCGGCGTGGTCGGTGTGTTCCTGGCCCTGGACCTGTTCCTGTTCTTCTTCTTCTGGGA  >  NC_002947/4668313‑4668608

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑