Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,316,890 C→G 7.7% 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,8900CG 7.7% 452.8 / 22.3 141V213L (GTG→CTG) nusAtranscription termination/antitermination L factor
Reads supporting (aligned to +/- strand):  ref base C (69/61);  new base G (5/6);  total (74/67)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 144 total aligned reads displayed.

CACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTCGCGCA  >  NC_000913/3316741‑3317039
                                                                                                                                                     |                                                                                                                                                     
  ccGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACtt                                                                                                                                                     >  1:1818240/1‑150 (MQ=255)
        gAACACGCGCGCCACGCATACCTACGCAAGCCCCTACCGGATCGATACGTTTATCGTTGGTTTTCCCCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCCATCACTTCTTCGCCGATTTCTGGCACTTCAATAc                                                                                                                                               >  2:1402871/1‑150 (MQ=255)
          acacGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACgg                                                                                                                                             <  2:993883/150‑1 (MQ=255)
          acacGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTAATCGTTGGATTTCCCCGCGATTTTCGCACGAGAAACCTGATCGCGAGCCGCTGCTTTAATTACAATCACTTCATCGCCGATATCTGGCACTTCATTACgg                                                                                                                                             >  2:2047652/1‑150 (MQ=255)
             cgcgcgCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAAc                                                                                                                                          <  1:1982876/150‑1 (MQ=255)
              gcgcgcCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTATTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACa                                                                                                                                         <  2:1101097/150‑1 (MQ=255)
               cgcgcCATGCATACCTACGCAAGCACCTACTGGAGCGATCCGTTTCTCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAg                                                                                                                                        <  1:1471516/150‑1 (MQ=255)
               cgcgcCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAg                                                                                                                                        <  1:37081/150‑1 (MQ=255)
                     aCGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGAt                                                                                                                                  <  2:1742998/150‑1 (MQ=255)
                        cATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAg                                                                                                                               <  2:2008001/150‑1 (MQ=255)
                         aTACCTCCGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGc                                                                                                                              <  2:867633/150‑1 (MQ=255)
                         aTACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGc                                                                                                                              <  1:684203/150‑1 (MQ=255)
                          tACCTACGCAAGCACCTACCTGATCGATACGTTTATCGTTGGTGTTCACCGCGATTTGCGCACGAGAACCCGGATCGCGAGCCGCGGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGACCAGTTCGATCAGCa                                                                                                                             <  1:115748/150‑1 (MQ=255)
                          aaCCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCa                                                                                                                             <  1:656276/149‑1 (MQ=255)
                            ccTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCAtt                                                                                                                           >  2:841154/1‑150 (MQ=255)
                            ccTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCAtt                                                                                                                           >  2:1720282/1‑150 (MQ=255)
                              tACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGCTCAGCATTTc                                                                                                                         >  2:1768827/1‑150 (MQ=255)
                              tACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTc                                                                                                                         >  2:948360/1‑150 (MQ=255)
                              tACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTc                                                                                                                         >  1:5267/1‑150 (MQ=255)
                              tACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTc                                                                                                                         >  1:2021520/1‑150 (MQ=255)
                              tACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTAGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTc                                                                                                                         >  2:1309678/1‑150 (MQ=255)
                                cGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCg                                                                                                                       <  2:468199/150‑1 (MQ=255)
                                cGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCg                                                                                                                       >  2:747978/1‑150 (MQ=255)
                                  cAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGc                                                                                                                     >  2:1545410/1‑150 (MQ=255)
                                  cAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGc                                                                                                                     >  2:1229064/1‑150 (MQ=255)
                                    aGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAAGTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCtt                                                                                                                   >  1:1246064/1‑150 (MQ=255)
                                         cTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAAc                                                                                                              >  1:533849/1‑150 (MQ=255)
                                           aCCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGATTTAATTTCAATCACTTCTTCGCCGATTTCTGGCTCTTCAATACGGAACAGTTCGATAAGCATTTGCGGCTTGGAACGa                                                                                                            >  1:557676/1‑150 (MQ=255)
                                            ccGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAg                                                                                                           <  2:2041198/150‑1 (MQ=255)
                                             cGGATCGTTACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGt                                                                                                          <  1:655804/150‑1 (MQ=255)
                                             cGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGt                                                                                                          <  2:2154287/150‑1 (MQ=255)
                                                 tCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACg                                                                                                      >  1:487142/1‑150 (MQ=255)
                                                  cGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGa                                                                                                     >  2:1471765/1‑150 (MQ=255)
                                                     tACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACa                                                                                                  >  1:33155/1‑150 (MQ=255)
                                                     tACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGACTTGGAACGAGTGACGAACa                                                                                                  >  2:1867106/1‑150 (MQ=255)
                                                      aCGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAg                                                                                                 <  1:1664440/150‑1 (MQ=255)
                                                       cGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGt                                                                                                >  1:184778/1‑150 (MQ=255)
                                                        gTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGtt                                                                                               >  1:938008/1‑150 (MQ=255)
                                                                 tGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAc                                                                                      <  2:990341/150‑1 (MQ=255)
                                                                  ggatttCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCAcg                                                                                     <  2:1766521/147‑1 (MQ=255)
                                                                       tCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCtt                                                                                >  2:1791945/1‑150 (MQ=255)
                                                                        cACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTc                                                                               >  1:1135170/1‑150 (MQ=255)
                                                                        cACCGCGATTTTCGCACGAGAACCCGGATAGCTAGCCGCTGCTTTAATTTCAATCACATCTTCGCCGATTTCTGGCACTTCAATACTGAACAGTTCGATCAGCATTTCCGGATTGGAACGAGTGACTAACAGTTGCGCGCCACTCGCtaa                                                                               >  1:1576369/1‑148 (MQ=255)
                                                                         aCCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTcc                                                                              >  1:34527/1‑150 (MQ=255)
                                                                         aCCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTcc                                                                              <  1:1313694/150‑1 (MQ=255)
                                                                          ccGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCCCGCGCATccc                                                                             >  1:501290/1‑149 (MQ=255)
                                                                          ccGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCg                                                                             >  2:909329/1‑150 (MQ=255)
                                                                          ccGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCg                                                                             >  1:2165555/1‑150 (MQ=255)
                                                                          ccGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCg                                                                             >  2:1020720/1‑150 (MQ=255)
                                                                            gcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCggg                                                                           <  2:1544324/150‑1 (MQ=255)
                                                                            gcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGATCAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCggg                                                                           <  2:258697/150‑1 (MQ=255)
                                                                            gcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCggg                                                                           <  1:263643/150‑1 (MQ=255)
                                                                            gcgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCGTAAATTTCACTCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGAGCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTGCCggg                                                                           <  2:2080566/150‑1 (MQ=255)
                                                                            gcgATTTTCGCACGAGAACCCGGAGCGCGAGCCGCTGCTTTAATTTCAATCACTTCATCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCggg                                                                           <  2:1945922/150‑1 (MQ=255)
                                                                             cgATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTTTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGc                                                                          <  2:1229615/150‑1 (MQ=255)
                                                                              gATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATATCAATCACTTCTTCGCCGATTTCTGGCACTTCAATCCGGAACAGTGCGATCAGCATTGCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCg                                                                         <  1:192041/150‑1 (MQ=255)
                                                                                    cGCACGAGAACCCGGATCGCGCGCCGCTGCTATAATTTCAATCACTTCTTCGCCGATTTCTGGCACTACCATACGGAACAGTTCGAGCAGCATTTCCGGCTTGGAACGAGTGACGAACAGATGCGCGCCACGCGCTTCCGGGCGAACGGa                                                                   >  1:399377/1‑150 (MQ=255)
                                                                                    cGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGa                                                                   >  2:2115774/1‑150 (MQ=255)
                                                                                      cACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAAt                                                                 >  2:1647731/1‑150 (MQ=255)
                                                                                      cACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAAt                                                                 >  1:2094612/1‑150 (MQ=255)
                                                                                        cGAGAACCCGGATCTCGAGCCGCTTCTTTAATTTCAATCACTTCTTCGCCTATTTCTGGCACTTCAATACGGAACAGTTCGATCTGCATTTCCTGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATag                                                               <  2:1436539/150‑1 (MQ=255)
                                                                                        cGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATag                                                               <  2:5267/150‑1 (MQ=255)
                                                                                             aCCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCAGGCACTACAAAACGGAACAGTTCGATCACCATATCCGGCTTGCAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAGCGGAATAGAGCAc                                                          >  2:2039892/1‑150 (MQ=255)
                                                                                              cccGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCCCGCGCTTCCGGGCGAACGGAATAGAGCACg                                                         >  2:244892/1‑150 (MQ=255)
                                                                                                cGGATCGCGAGCCGCTGCTTTAATATCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGCACCGTTCGATCAGCATTTCCGGCGTGGAACGAGTGAGGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCCCGcc                                                       >  1:172212/1‑150 (MQ=255)
                                                                                                cGGAGCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTCCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGcc                                                       <  2:256294/150‑1 (MQ=255)
                                                                                                 ggATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCa                                                      <  2:650675/150‑1 (MQ=255)
                                                                                                   aTCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACg                                                    <  2:538912/150‑1 (MQ=255)
                                                                                                    tCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGa                                                   >  1:2062741/1‑150 (MQ=255)
                                                                                                      gcgAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGCCCAGGTCGATCAGCATTTCCGGCTTGGAACGAGTGACGCACAGTTGCGCGCCCCTCGCTTCCGGGCGAAAGGAATAGAGCACGCCACGAAc                                                 >  2:375610/1‑150 (MQ=255)
                                                                                                      gcgAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGATCGAGTGACGAACAGTTGCGAGCCCCGCGCTTCCGTGCGAACGGAATAGAGCACGCCACGAAc                                                 >  1:1017076/1‑150 (MQ=255)
                                                                                                      gcgAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAAc                                                 >  1:2076432/1‑150 (MQ=255)
                                                                                                       cgAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAAcg                                                <  2:2063658/150‑1 (MQ=255)
                                                                                                             gctgctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCAGTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCg                                          <  2:1855284/150‑1 (MQ=255)
                                                                                                             gctgctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCg                                          >  1:1798140/1‑150 (MQ=255)
                                                                                                             gatgctTTAATATCAATCACTTCTTCGCCGATTTCTGTAACTTCAATCCGGAACAGATCGATCAGAATTTCAGGCTTGGAACGTGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGATCACGCCACGAACGCGGTCg                                          <  2:684522/148‑1 (MQ=255)
                                                                                                               tgctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGcc                                        <  2:1303878/150‑1 (MQ=255)
                                                                                                                 ctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGTAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACTAACGCGGTCGCCAg                                      >  1:1332876/1‑150 (MQ=255)
                                                                                                                 ctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAg                                      <  1:1229064/150‑1 (MQ=255)
                                                                                                                 ctTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAg                                      >  1:2147402/1‑150 (MQ=255)
                                                                                                                  tttACTTCCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAgg                                     <  2:1623169/150‑1 (MQ=255)
                                                                                                                     aaTTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCg                                  <  1:1436761/150‑1 (MQ=255)
                                                                                                                     aaTTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGGGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCg                                  <  1:772466/150‑1 (MQ=255)
                                                                                                                      aTTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCCCGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCgg                                 >  2:381215/1‑150 (MQ=255)
                                                                                                                       tttCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGa                                >  1:696335/1‑150 (MQ=255)
                                                                                                                       tttCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGACACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGa                                <  1:790001/150‑1 (MQ=255)
                                                                                                                         tCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAg                              >  2:1781983/1‑150 (MQ=255)
                                                                                                                         tCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAg                              <  1:1853211/150‑1 (MQ=255)
                                                                                                                         tCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGATTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGATTAGAGCCCGCCCCGAACGCGGGCGCCAGGGCGGAAg                              >  2:1790645/1‑150 (MQ=255)
                                                                                                                              cACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTc                         >  1:1369387/1‑150 (MQ=255)
                                                                                                                              cACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTc                         >  2:577657/1‑150 (MQ=255)
                                                                                                                                     tcGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgc                  >  1:639207/1‑150 (MQ=255)
                                                                                                                                     tcGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgc                  <  1:624574/150‑1 (MQ=255)
                                                                                                                                      cGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGATTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGCACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgca                 >  1:1266820/1‑150 (MQ=255)
                                                                                                                                       gCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCACTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGgcag                >  1:1275745/1‑150 (MQ=255)
                                                                                                                                          gATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCat             >  1:1613495/1‑150 (MQ=255)
                                                                                                                                           aTTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCata            <  1:2115774/150‑1 (MQ=255)
                                                                                                                                               cTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCtt        <  1:34538/150‑1 (MQ=255)
                                                                                                                                                     cTTCAATACGGAACAGTTCGCTCAGCCTTTCCGGCTGGGAACGCGTGCCGAACAGTTGCGCGCCACGCGCGTCCGGGCGAACGGAATAGAGCCCGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATAGCTTCGCGCa  <  1:1395351/150‑1 (MQ=255)
                                                                                                                                                     cTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGGCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTCGCGCa  <  1:1402871/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGAACAGTTCGATCAGCATTTCCGGCTTGGAACGAGTGACGAACAGTTGCGCGCCACGCGCTTCCGGGCGAACGGAATAGAGCACGCCACGAACGCGGTCGCCAGGGCGGAAGTTTTCACGCGGCAGCATATCTTCGCGCA  >  NC_000913/3316741‑3317039

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