Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,316,890 C→G 22.9% V213L (GTG→CTG)  nusA ← transcription termination/antitermination L factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,316,8900CG22.9% 654.7 / 213.3 270V213L (GTG→CTG) nusAtranscription termination/antitermination L factor
Reads supporting (aligned to +/- strand):  ref base C (116/92);  new base G (29/33);  total (145/125)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 273 total aligned reads displayed.

CACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTCGCGCA  >  NC_000913/3316741‑3317039
                                                                                                                                                     |                                                                                                                                                     
cacCGCCTGCACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAc                                                                                                                                                       <  2:88498/150‑1 (MQ=255)
 acCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACt                                                                                                                                                      >  2:1069711/1‑150 (MQ=255)
  ccGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCACTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACtt                                                                                                                                                     <  1:1052374/150‑1 (MQ=255)
    gCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCa                                                                                                                                                   >  1:851417/1‑150 (MQ=255)
     ccTGAACACGCGCGCCACTCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCaa                                                                                                                                                  <  2:548922/150‑1 (MQ=255)
         aaCACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACg                                                                                                                                              >  1:1472127/1‑150 (MQ=255)
          acacGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACgg                                                                                                                                             <  2:1932756/150‑1 (MQ=255)
            acGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGaa                                                                                                                                           >  1:1449782/1‑150 (MQ=255)
            acGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGaa                                                                                                                                           <  1:1595629/150‑1 (MQ=255)
            acGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGaa                                                                                                                                           >  1:605227/1‑150 (MQ=255)
             cgcgcgCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAAc                                                                                                                                          >  2:454594/1‑150 (MQ=255)
             cgcgcgCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAAc                                                                                                                                          >  2:1186941/1‑150 (MQ=255)
             cgcgcgCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGacc                                                                                                                                          >  2:2612143/1‑148 (MQ=255)
             cgcgcgCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAAc                                                                                                                                          >  2:2612258/1‑150 (MQ=255)
                     aCGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGAt                                                                                                                                  <  2:1271911/150‑1 (MQ=255)
                     aCGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGAt                                                                                                                                  <  1:349074/150‑1 (MQ=255)
                      cGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATc                                                                                                                                 <  2:2316388/150‑1 (MQ=255)
                       gCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCa                                                                                                                                >  1:749384/1‑150 (MQ=255)
                       gCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCa                                                                                                                                <  1:1789251/150‑1 (MQ=255)
                         aTACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGc                                                                                                                              <  2:768963/150‑1 (MQ=255)
                         aTACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGc                                                                                                                              <  1:937779/150‑1 (MQ=255)
                         aTACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGc                                                                                                                              <  1:101709/150‑1 (MQ=255)
                         aTACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGc                                                                                                                              <  1:1157911/150‑1 (MQ=255)
                             cTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCAttt                                                                                                                          >  1:1305029/1‑150 (MQ=255)
                              tACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTc                                                                                                                         <  2:2809840/150‑1 (MQ=255)
                                    aGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCtt                                                                                                                   <  1:1535178/150‑1 (MQ=255)
                                      cACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTgg                                                                                                                 >  2:1572292/1‑150 (MQ=255)
                                        ccTACCGGATCGATACGTTTATCGTTGGTTTTCCCCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGaa                                                                                                               >  2:358579/1‑150 (MQ=255)
                                         cTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAAc                                                                                                              >  2:1960038/1‑150 (MQ=255)
                                         cTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAAc                                                                                                              >  1:2253881/1‑150 (MQ=255)
                                          tACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACg                                                                                                             >  2:279333/1‑150 (MQ=255)
                                          tACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACg                                                                                                             >  2:2237296/1‑150 (MQ=255)
                                              ggATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTg                                                                                                         >  1:465034/1‑150 (MQ=255)
                                               gATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGa                                                                                                        <  2:210129/150‑1 (MQ=255)
                                               gATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGa                                                                                                        >  2:1793462/1‑150 (MQ=255)
                                               gATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGa                                                                                                        <  1:1830856/150‑1 (MQ=255)
                                               gATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGa                                                                                                        <  2:2380337/150‑1 (MQ=255)
                                                   gATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGaa                                                                                                    <  2:1472127/150‑1 (MQ=255)
                                                    aTACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAAc                                                                                                   >  2:2444171/1‑150 (MQ=255)
                                                     tACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACa                                                                                                  <  1:2069285/150‑1 (MQ=255)
                                                      aCGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAg                                                                                                 >  2:398434/1‑150 (MQ=255)
                                                         tttATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTg                                                                                              >  1:2118766/1‑150 (MQ=255)
                                                         tttATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTg                                                                                              >  1:2285891/1‑150 (MQ=255)
                                                          ttATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTgc                                                                                             >  2:262051/1‑150 (MQ=255)
                                                            aTCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTgcgc                                                                                           >  2:1113030/1‑150 (MQ=255)
                                                             atcgcggTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAGCAGTTgcgcg                                                                                          >  2:2587498/6‑150 (MQ=255)
                                                              cGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTgcgcgc                                                                                         <  1:2334542/150‑1 (MQ=255)
                                                               gTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGcc                                                                                        >  1:1660827/1‑150 (MQ=255)
                                                               gTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGcc                                                                                        >  1:2662451/1‑150 (MQ=255)
                                                               gTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGcc                                                                                        >  2:1189640/1‑150 (MQ=255)
                                                                ttGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCa                                                                                       >  1:2559479/1‑150 (MQ=255)
                                                                 tGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAc                                                                                      <  1:425654/150‑1 (MQ=255)
                                                                    ttttCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAcgcg                                                                                   >  2:2818589/1‑150 (MQ=255)
                                                                     tttCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACgcgc                                                                                  >  2:289347/1‑150 (MQ=255)
                                                                     tttCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACgcgc                                                                                  <  2:1665680/150‑1 (MQ=255)
                                                                       tCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCtt                                                                                <  2:1030032/150‑1 (MQ=255)
                                                                          ccGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCg                                                                             <  2:2486550/150‑1 (MQ=255)
                                                                          ccGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGAACAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCg                                                                             <  2:1864507/150‑1 (MQ=255)
                                                                            gcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCggg                                                                           <  2:2568140/150‑1 (MQ=255)
                                                                            gcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCggg                                                                           >  1:2922625/1‑150 (MQ=255)
                                                                                     gCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGaa                                                                  <  2:1325840/150‑1 (MQ=255)
                                                                                     gCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGaa                                                                  <  2:1598358/150‑1 (MQ=255)
                                                                                      cACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAAt                                                                 >  2:476397/1‑150 (MQ=255)
                                                                                       aCGAGAACCCGGCTCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATa                                                                <  2:2935021/150‑1 (MQ=255)
                                                                                       aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATa                                                                <  1:2745630/150‑1 (MQ=255)
                                                                                       aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATa                                                                <  2:1719645/150‑1 (MQ=255)
                                                                                       aCGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTACGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATa                                                                >  1:1229934/1‑150 (MQ=255)
                                                                                        cGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTAAATCACTTCTTCGCCGATTTATGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGGGTGACGAACAGTTGCGCGCCACGCGCTGCCGGGCGAACGGAATag                                                               >  1:2212883/1‑150 (MQ=255)
                                                                                         gagaACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATaga                                                              <  2:67768/150‑1 (MQ=255)
                                                                                            aaCCCGGATCGCGAGCCGCTGGTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCa                                                           >  1:2131675/1‑150 (MQ=255)
                                                                                               ccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGc                                                        >  2:1807894/1‑150 (MQ=255)
                                                                                               ccGGATCGCGAGCCGCGGCTTTAATTTCAATCAATTCTTCGCCGATTTCTGGCACTTCAATAAGGTACAGTTCGATAAGCATTTACGGCGTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCATCCGGGCGCACGGAATAGAGCACGc                                                        >  1:2029094/1‑150 (MQ=255)
                                                                                                 ggATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCa                                                      <  1:890433/150‑1 (MQ=255)
                                                                                                 ggATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCa                                                      >  2:962723/1‑150 (MQ=255)
                                                                                                             gctgctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCg                                          >  2:1292740/1‑150 (MQ=255)
                                                                                                             gctgctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCg                                          <  1:645537/150‑1 (MQ=255)
                                                                                                              ctgctTTAATTTCAATCACTTCGTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTcac                                         <  2:2168215/150‑3 (MQ=255)
                                                                                                               tgctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGcc                                        >  2:631827/1‑150 (MQ=255)
                                                                                                               tgctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGcc                                        <  2:1565048/150‑1 (MQ=255)
                                                                                                                gctTTTTTTTCAATCACTTCTTCGCCGTTTTCTGGCTCTTCAATTAGGATCAGTTCGATCTGGATTTCCTGCTTTGAAGGAGTGACGTTCAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCa                                       <  2:2429483/150‑1 (MQ=255)
                                                                                                                gctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCa                                       <  1:914781/150‑1 (MQ=255)
                                                                                                                  tttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAgg                                     >  1:667235/1‑150 (MQ=255)
                                                                                                                   ttAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAggg                                    <  2:1530883/150‑1 (MQ=255)
                                                                                                                    tAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGc                                   <  1:1207004/150‑1 (MQ=255)
                                                                                                                     aaTTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCg                                  <  1:2508926/150‑1 (MQ=255)
                                                                                                                       tttCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGa                                <  2:644723/150‑1 (MQ=255)
                                                                                                                           aaTCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGtt                            >  2:2507782/1‑150 (MQ=255)
                                                                                                                             tCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGtttt                          <  1:1957759/150‑1 (MQ=255)
                                                                                                                             tCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGtttt                          <  1:1154843/150‑1 (MQ=255)
                                                                                                                              cACTTCTTCGCCGATTTCTGGCACTTCAATCCGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTc                         <  2:2638084/150‑1 (MQ=255)
                                                                                                                                 ttcttcGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCAcg                      >  1:2427927/1‑150 (MQ=255)
                                                                                                                                  tcttcGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCAcgc                     <  1:1281906/150‑1 (MQ=255)
                                                                                                                                       gCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgcag                >  1:2839360/1‑150 (MQ=255)
                                                                                                                                        ccGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgcagc               <  2:2112700/150‑1 (MQ=255)
                                                                                                                                        ccGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgcagc               >  2:2024923/1‑150 (MQ=255)
                                                                                                                                          gATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCat             <  1:1113030/150‑1 (MQ=255)
                                                                                                                                           aTTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCata            >  1:2638622/1‑150 (MQ=255)
                                                                                                                                            tttCTGGCACTTCAATCCGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCatat           <  2:1391292/150‑1 (MQ=255)
                                                                                                                                                tGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTc       >  2:1841945/1‑150 (MQ=255)
                                                                                                                                                     cTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTCGCGCa  >  1:2246277/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTCGCGCA  >  NC_000913/3316741‑3317039

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