Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,021,409 C→G 5.6% *314S (TGA→TCA)  PP_3547 ← NAD(P)‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,021,4090CG 5.6% 323.5 / 8.2 107*314S (TGA→TCA) PP_3547NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (55/45);  new base G (4/2);  total (60/47)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 111 total aligned reads displayed.

CTCAAGGTGCTGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGTCA  >  NC_002947/4021260‑4021558
                                                                                                                                                      |                                                                                                                                                     
cTCAAGGTGCTGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTc                                                                                                                                                       >  1:1070131/1‑150 (MQ=255)
 tCAAGGTGCTGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCa                                                                                                                                                      >  2:741291/1‑150 (MQ=255)
 tCAAGGTGCTGTTGACCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCa                                                                                                                                                      >  1:1352289/1‑150 (MQ=255)
          tGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTGAGcgcccgtc                                                                                                                                             >  2:1747499/1‑148 (MQ=255)
          tGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTGAGcgcccgtc                                                                                                                                             <  1:1747499/150‑3 (MQ=255)
          tGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGcgcccgcc                                                                                                                                             <  1:1635665/150‑1 (MQ=255)
            ttGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGGCGGGCGCTGAGcgcccgtcaa                                                                                                                                           <  2:1784672/150‑5 (MQ=255)
                tCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAg                                                                                                                                       <  1:1349575/150‑1 (MQ=255)
                tCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCAGCTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAg                                                                                                                                       >  1:1345367/1‑150 (MQ=255)
                 cACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAgg                                                                                                                                      >  1:953278/1‑150 (MQ=255)
                  aCCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGt                                                                                                                                     >  2:2027444/1‑150 (MQ=255)
                   ccGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGACGTTAGGACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTc                                                                                                                                    >  1:885989/1‑150 (MQ=255)
                        tACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGa                                                                                                                               <  2:1352289/150‑1 (MQ=255)
                        tACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGa                                                                                                                               <  2:1221433/150‑1 (MQ=255)
                        tACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGa                                                                                                                               <  1:491957/150‑1 (MQ=255)
                        tACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGa                                                                                                                               >  1:1821079/1‑150 (MQ=255)
                        tACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGAGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGa                                                                                                                               >  1:1023100/1‑150 (MQ=255)
                        tACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACCTTGGTCTCTCCGAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGa                                                                                                                               <  1:24817/150‑1 (MQ=255)
                           gATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGc                                                                                                                            >  1:110358/1‑150 (MQ=255)
                            aTGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCt                                                                                                                           >  2:280551/1‑150 (MQ=255)
                                gagTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCACCGCCCCGGTCGATGAGGCTGCCg                                                                                                                       >  2:381619/1‑150 (MQ=255)
                                   tGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTg                                                                                                                    <  1:741291/150‑1 (MQ=255)
                                   tGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTg                                                                                                                    <  1:715174/150‑1 (MQ=255)
                                   tGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTg                                                                                                                    <  2:1117052/150‑1 (MQ=255)
                                      ggcgTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTAACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCATCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCa                                                                                                                 >  1:178309/1‑150 (MQ=255)
                                               cGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATagg                                                                                                             <  1:586062/145‑1 (MQ=255)
                                               cGACGGCCAACCGCTGCAAGGTACGTTGGTGGTGACCAAGCCGATCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGAGGAGGCTGCCGGTGGCATAGGAGGcc                                                                                                        >  2:372778/1‑149 (MQ=255)
                                               cGACGGCCAACCGCTGAAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATagg                                                                                                             >  2:586062/1‑145 (MQ=255)
                                                gACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCtt                                                                                                       >  1:1612548/1‑150 (MQ=255)
                                                gACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGACCAGCAATTCCTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCtt                                                                                                       >  1:1861946/1‑150 (MQ=255)
                                                 aCGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAATCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCttt                                                                                                      <  1:1347455/150‑1 (MQ=255)
                                                 aCGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCttt                                                                                                      <  1:202325/150‑1 (MQ=255)
                                                   ggagaaCCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGt                                                                                                    >  2:222186/5‑150 (MQ=255)
                                                    gCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTc                                                                                                   >  2:1274692/1‑150 (MQ=255)
                                                      cAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCttt                                                                                                      >  2:1782840/1‑145 (MQ=255)
                                                      cAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCttt                                                                                                      <  1:1782840/145‑1 (MQ=255)
                                                      cAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCg                                                                                                 >  2:1557794/1‑150 (MQ=255)
                                                        aCCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGaa                                                                                               >  1:467439/1‑150 (MQ=255)
                                                         ccGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGaaa                                                                                              <  1:1274692/150‑1 (MQ=255)
                                                         ccGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGaaa                                                                                              <  2:2421/150‑1 (MQ=255)
                                                         ccGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGAGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGaaa                                                                                              <  2:1023100/150‑1 (MQ=255)
                                                         ccGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTAGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGaaa                                                                                              >  2:1328452/1‑150 (MQ=255)
                                                                  ggTACGCTGGTCTGGCCCAAGCCGTTCGAACTGCAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCaga                                                                                     <  2:662124/150‑1 (MQ=255)
                                                                     aCGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGa                                                                                  <  2:690820/150‑1 (MQ=255)
                                                                       gTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATg                                                                                <  2:151146/150‑1 (MQ=255)
                                                                           gtgtTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCt                                                                            >  2:2124331/1‑150 (MQ=255)
                                                                            tgtTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTc                                                                           >  2:857993/1‑150 (MQ=255)
                                                                                 aCCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCa                                                                      >  1:2031520/1‑150 (MQ=255)
                                                                                          ttCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCATCGGCCACGTCCTCgg                                                             >  2:545718/1‑150 (MQ=255)
                                                                                               aCTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGc                                                        <  1:857993/150‑1 (MQ=255)
                                                                                                 tGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCgg                                                                         >  1:1815071/1‑131 (MQ=255)
                                                                                                 tGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCgg                                                                         <  2:1815071/131‑1 (MQ=255)
                                                                                                  ggAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGATAACCGGGCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGcccc                                                     >  2:1796934/1‑150 (MQ=255)
                                                                                                   gAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTc                                                                                                   <  2:205792/103‑1 (MQ=255)
                                                                                                   gAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTc                                                                                                   >  1:205792/1‑103 (MQ=255)
                                                                                                      gCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTc                                                 >  2:861119/1‑150 (MQ=255)
                                                                                                        aCTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCAc                                                                     <  1:1802761/128‑1 (MQ=255)
                                                                                                        aCTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCAc                                                                     >  2:1802761/1‑128 (MQ=255)
                                                                                                        aCTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCgg                                               >  1:382160/1‑150 (MQ=255)
                                                                                                          tGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGcc                                             >  2:134660/1‑150 (MQ=255)
                                                                                                            gCATAGCGGGTTCGCGAGCTGCTGGAACCTTGACTGCTGCGCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCa                                           >  1:2410938/1‑150 (MQ=255)
                                                                                                            gCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGTCACTTCCTCGGGACGCCCCCCTCGGCCCa                                           >  2:1803141/1‑150 (MQ=255)
                                                                                                             cAGAGCGGGTTCGCGAGCTGTTGGAACCTTGACGGGCGCTCATCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCAt                                          <  2:178309/150‑1 (MQ=255)
                                                                                                               gagCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATgg                                        <  2:1889610/150‑1 (MQ=255)
                                                                                                               gagCGGGTTCGCGAGCTGCTGGAACCTGGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATgg                                        <  1:506286/150‑1 (MQ=255)
                                                                                                                 gCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGc                                      >  1:515868/1‑150 (MQ=255)
                                                                                                                 gCGGGTTCGCGAGATGATGGAAACTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGc                                      >  2:627870/1‑150 (MQ=255)
                                                                                                                        cgcgAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAg                               <  1:1895699/150‑1 (MQ=255)
                                                                                                                          cgAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCt                                                                                                        <  2:333442/75‑1 (MQ=255)
                                                                                                                          cgAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCt                                                                                                        >  1:333442/1‑75 (MQ=255)
                                                                                                                          cgAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGt                             <  2:885989/150‑1 (MQ=255)
                                                                                                                          cgAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGCCCGAAAGCAACGATAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGt                             <  1:2124331/150‑1 (MQ=255)
                                                                                                                           gAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGtt                            <  2:515868/150‑1 (MQ=255)
                                                                                                                           gAGCTGCTGGAACCTTGACGAGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGtt                            >  2:892021/1‑150 (MQ=255)
                                                                                                                            aGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTc                           <  1:2027444/150‑1 (MQ=255)
                                                                                                                             gctgctGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCg                          >  1:354611/1‑150 (MQ=255)
                                                                                                                              ctgctgGACCCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTAGCTGGGGTTGACGGTGTGATATGATGCGTTGTCAGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGa                         <  1:280551/150‑1 (MQ=255)
                                                                                                                               tgctgGCACCTTGACGGGCGCTCAGCGCCCGCCGCCAATCACGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAg                        <  1:1705854/150‑1 (MQ=255)
                                                                                                                               tgctgGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAg                        >  1:2305335/1‑150 (MQ=255)
                                                                                                                                 ctgGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGcc                                                       <  1:638521/117‑1 (MQ=255)
                                                                                                                                 ctgGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGcc                                                       >  2:638521/1‑117 (MQ=255)
                                                                                                                                   ggAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCt                                                                                                        >  1:60824/1‑66 (MQ=255)
                                                                                                                                    gAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCTAGCCAGCTTCGAGCTTGc                   >  2:1262285/1‑150 (MQ=255)
                                                                                                                                    gAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGc                   >  2:1921938/1‑150 (MQ=255)
                                                                                                                                        ccTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGa                >  2:170547/1‑150 (MQ=255)
                                                                                                                                        ccTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGa                <  2:110358/150‑1 (MQ=255)
                                                                                                                                         cTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGAc               >  1:109982/1‑150 (MQ=255)
                                                                                                                                         cTTGACGGGCGCCCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGAc               <  1:2006765/150‑1 (MQ=255)
                                                                                                                                          ttGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGAcg              >  1:452674/1‑150 (MQ=255)
                                                                                                                                           tGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGAcgc             <  1:1557794/150‑1 (MQ=255)
                                                                                                                                           gggcGGGCGCTGAGCGCCCGTCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGAcgc             >  2:2158788/4‑150 (MQ=255)
                                                                                                                                             gcGGGCGCTGAGCGCCCGCCGCCAAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTAATAGCAACACCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCgg           >  2:1638405/2‑150 (MQ=255)
                                                                                                                                               gggCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTc         >  1:328342/1‑150 (MQ=255)
                                                                                                                                               gggCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTc         <  2:417217/150‑1 (MQ=255)
                                                                                                                                               gggCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTc         >  2:2128655/1‑150 (MQ=255)
                                                                                                                                                  cgcTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCgg      <  1:627870/150‑1 (MQ=255)
                                                                                                                                                     tCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGTc   >  2:150213/1‑150 (MQ=255)
                                                                                                                                                      cAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGTCa  <  1:1328452/150‑1 (MQ=255)
                                                                                                                                                      cAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGTCa  <  1:170547/150‑1 (MQ=255)
                                                                                                                                                      cAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGTCa  <  2:1861946/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTCAAGGTGCTGTTGATCACCGGCTACGATGAGAGTGCGGCGTTCAGCGACGGCCAACCGCTGCAAGGTACGTTGGTGTTGACCAAGCCGTTCGAACTGGAAGCACTGGCAGAGCGGGTTCGCGAGCTGCTGGAACCTTGACGGGCGCTCAGCGCCCGCCGCCCAGGTCGATGAGGCTGCCGGTGGCATAGGAGGCTTTGTCCGAAAGCAACCAGAGGATGGCCTCGGCCACTTCCTCGGGACGCCCCCCTCGGCCCATGGGCAAGCCAGGTTCGAGCTTGCTGACGCGGTCCGGGTCA  >  NC_002947/4021260‑4021558

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