Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 105,046 A→C 5.2% Q355H (CAA→CAC)  ftsA → cell division protein FtsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913105,0460AC 5.2% 534.2 / 14.8 168Q355H (CAA→CAC) ftsAcell division protein FtsA
Reads supporting (aligned to +/- strand):  ref base A (85/74);  new base C (7/2);  total (92/76)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

AACCTGGTCAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCATCAG  >  NC_000913/104897‑105194
                                                                                                                                                     |                                                                                                                                                    
aaCCTGGTCAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCaa                                                                                                                                                      >  1:2017309/1‑150 (MQ=255)
  ccTGGTCAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGt                                                                                                                                                    >  2:2207191/1‑150 (MQ=255)
   cTGGTCAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTg                                                                                                                                                   >  1:1894980/1‑150 (MQ=255)
     ggTCAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCg                                                                                                                                                 <  2:2374006/150‑1 (MQ=255)
       tCAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTa                                                                                                                                               >  1:1986151/1‑150 (MQ=255)
        cAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTAt                                                                                                                                              >  1:2440407/1‑150 (MQ=255)
         aaCGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATc                                                                                                                                             <  1:112571/150‑1 (MQ=255)
         aaCGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATc                                                                                                                                             >  2:1242873/1‑150 (MQ=255)
           cGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCgg                                                                                                                                           >  1:216654/1‑150 (MQ=255)
           cGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCgg                                                                                                                                           >  1:647250/1‑150 (MQ=255)
             aaGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCACTTGCGTAtgaaac                                                                                                                                         <  2:2194551/150‑6 (MQ=255)
              agagATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGgcgc                                                                                                                                        <  1:2113347/150‑1 (MQ=255)
                   tatTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGcgccgc                                                                                                                                   >  1:1336120/1‑150 (MQ=255)
                         agttgcagGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACa                                                                                                                             <  1:1384540/150‑1 (MQ=255)
                          gttgcagGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACAt                                                                                                                            <  1:2150158/150‑1 (MQ=255)
                                ggAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCgg                                                                                                                      <  2:1751782/150‑1 (MQ=255)
                                ggAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCgg                                                                                                                      >  1:1616542/1‑150 (MQ=255)
                                ggAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCgg                                                                                                                      <  1:946730/150‑1 (MQ=255)
                                ggAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCgg                                                                                                                      >  1:424836/1‑150 (MQ=255)
                                ggAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCgg                                                                                                                      <  1:91343/150‑1 (MQ=255)
                                ggAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCgg                                                                                                                      >  1:110475/1‑150 (MQ=255)
                                  aaaaGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGtt                                                                                                                    >  1:467014/1‑150 (MQ=255)
                                    aaGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTa                                                                                                                  <  1:2189361/150‑1 (MQ=255)
                                     aGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTaa                                                                                                                 >  1:358955/1‑150 (MQ=255)
                                       cTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACg                                                                                                               <  1:1752160/150‑1 (MQ=255)
                                           gCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGAtt                                                                                                           >  1:1643202/1‑150 (MQ=255)
                                             caacaaGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTAt                                                                                                         <  1:1173420/150‑1 (MQ=255)
                                              aacaaGGGGTTAAACATCACCTGGCGGCAGTCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAATGGATTATg                                                                                                        >  2:1584431/1‑150 (MQ=255)
                                                 aaGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTc                                                                                                     <  2:2064252/150‑1 (MQ=255)
                                                 aaGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTc                                                                                                     <  2:567775/150‑1 (MQ=255)
                                                  aGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCa                                                                                                    <  2:358955/150‑1 (MQ=255)
                                                  aGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCa                                                                                                    >  1:1022024/1‑150 (MQ=255)
                                                           aCATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGta                                                                                           >  1:832595/1‑150 (MQ=255)
                                                            cATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGtat                                                                                          <  2:1665789/150‑1 (MQ=255)
                                                              tCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGtatta                                                                                        >  1:651986/1‑150 (MQ=255)
                                                                 ccTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTc                                                                                     >  1:980880/1‑150 (MQ=255)
                                                                 ccTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGACGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTc                                                                                     <  1:596128/150‑1 (MQ=255)
                                                                  cTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGTGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCg                                                                                    <  2:859344/150‑1 (MQ=255)
                                                                      cggcAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGAcgg                                                                                >  2:151133/1‑150 (MQ=255)
                                                                      cggcAGGCATTGTATTAACCGGT‑GCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGAcggc                                                                               <  1:454286/150‑1 (MQ=255)
                                                                           ggcaTTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTTCGTATCGGCTCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTg                                                                           <  2:879557/150‑1 (MQ=255)
                                                                           ggcaTTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTg                                                                           >  2:2333032/1‑150 (MQ=255)
                                                                            gcaTTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTgg                                                                          <  2:1022024/150‑1 (MQ=255)
                                                                                    ttAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCtt                                                                  <  2:961763/150‑1 (MQ=255)
                                                                                    ttAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCtt                                                                  <  1:953986/150‑1 (MQ=255)
                                                                                      aaCCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCa                                                                >  1:858037/1‑150 (MQ=255)
                                                                                      aaCCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCa                                                                >  1:1236710/1‑150 (MQ=255)
                                                                                      aaCCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCa                                                                >  1:1294051/1‑150 (MQ=255)
                                                                                       aCCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCAc                                                               >  1:507402/1‑150 (MQ=255)
                                                                                       aCCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCAc                                                               >  2:1629503/1‑150 (MQ=255)
                                                                                         cGGTGGCCCAGCGCAGATCTAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTa                                                             >  2:770606/1‑150 (MQ=255)
                                                                                             ggCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATggg                                                         <  2:891650/150‑1 (MQ=255)
                                                                                             ggCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATggg                                                         <  1:2081556/150‑1 (MQ=255)
                                                                                             ggCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATggg                                                         <  2:1986151/150‑1 (MQ=255)
                                                                                               cgcagcgcagATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGaa                                                       >  1:1486398/1‑150 (MQ=255)
                                                                                               cgcagcgcagATAGAACGTCTTGCAGCCTGTGATCACCGAGTGTTTCATACGCAAGTGCGTATCGGAGCGCCGATGAACATTACAGGTTTAACGCATTATGCTCAGGAGCCGTATTATTAGACGGCGGTGGGATTGCTTCACTATGGGaa                                                       >  2:1196726/1‑150 (MQ=255)
                                                                                                   gcgcagATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAagag                                                   >  2:1752035/1‑150 (MQ=255)
                                                                                                   gcgcagATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAagag                                                   <  2:1384612/150‑1 (MQ=255)
                                                                                                     gcagATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCAc                                                               >  1:346110/1‑136 (MQ=255)
                                                                                                     gcagATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCAc                                                               <  2:346110/136‑1 (MQ=255)
                                                                                                     gcagATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTc                                                 >  1:2014802/1‑150 (MQ=255)
                                                                                                      cagATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTca                                                <  1:2051577/150‑1 (MQ=255)
                                                                                                        gATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTcaca                                              <  1:931505/150‑1 (MQ=255)
                                                                                                           cGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCt                                           <  2:2273048/150‑1 (MQ=255)
                                                                                                           cGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGCGAAAGAGTCACATCt                                           >  1:1955950/1‑150 (MQ=255)
                                                                                                            gAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCtt                                          <  2:467014/150‑1 (MQ=255)
                                                                                                              aGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTaa                                        >  2:1430697/1‑150 (MQ=255)
                                                                                                               ggTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAAc                                       >  2:1863646/1‑150 (MQ=255)
                                                                                                                 tCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGTTTCACTATGGGAAAGAGTCACATCTTAACgg                                     >  2:1107495/1‑150 (MQ=255)
                                                                                                                   ttGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTg                                   <  1:1629503/150‑1 (MQ=255)
                                                                                                                   ttGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTg                                   <  1:1863646/150‑1 (MQ=255)
                                                                                                                    tGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGa                                  >  1:1594310/1‑150 (MQ=255)
                                                                                                                    tGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGa                                  >  2:2094367/1‑150 (MQ=255)
                                                                                                                    tGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGa                                  >  1:278845/1‑150 (MQ=255)
                                                                                                                    tGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATATTAACGGTGa                                  >  1:1594314/1‑150 (MQ=255)
                                                                                                                     gCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGaa                                 >  2:1809902/1‑150 (MQ=255)
                                                                                                                      cAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAg                                >  1:2438244/1‑150 (MQ=255)
                                                                                                                          cTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGa                            >  1:1447119/1‑150 (MQ=255)
                                                                                                                             tgCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGt                         <  2:1336120/150‑1 (MQ=255)
                                                                                                                             tgCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGt                         <  2:424836/150‑1 (MQ=255)
                                                                                                                              gCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTa                        >  2:91908/1‑150 (MQ=255)
                                                                                                                               cTCAGCGCGTGTTTCATACGTAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAg                       >  1:406493/1‑150 (MQ=255)
                                                                                                                               cTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAg                       <  2:2440407/150‑1 (MQ=255)
                                                                                                                                  aGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGaaa                    >  2:47317/1‑150 (MQ=255)
                                                                                                                                  aGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGaaa                    >  2:802480/1‑150 (MQ=255)
                                                                                                                                    cgcgTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGaaaaa                  <  1:1242873/150‑1 (MQ=255)
                                                                                                                                       gtgtTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACgt               <  2:1177552/150‑1 (MQ=255)
                                                                                                                                       cgcgccgaTACGCACTTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACgt               >  1:2234036/8‑150 (MQ=255)
                                                                                                                                         gtTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACgtgt             <  1:1673570/150‑1 (MQ=255)
                                                                                                                                         gtTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACgtgt             >  1:1009120/1‑150 (MQ=255)
                                                                                                                                          gccgaTACGCACTTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGtt            >  1:2318741/5‑150 (MQ=255)
                                                                                                                                           ttCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTa           >  2:112368/1‑150 (MQ=255)
                                                                                                                                           ccgaTACGCACTTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTAACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTa           >  2:1570620/4‑150 (MQ=255)
                                                                                                                                            tCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTAc          >  2:1177908/1‑150 (MQ=255)
                                                                                                                                             cATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGa                      >  2:647585/1‑137 (MQ=255)
                                                                                                                                              aTACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAg        >  1:794976/1‑150 (MQ=255)
                                                                                                                                               tACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGc       <  2:1486398/150‑1 (MQ=255)
                                                                                                                                                 cGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCAt     <  2:1527560/150‑1 (MQ=255)
                                                                                                                                                 cGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCAt     <  2:2162833/150‑1 (MQ=255)
                                                                                                                                                  gCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCATc    <  1:1039104/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACCTGGTCAACGAAGAGATATTGCAGTTGCAGGAAAAGCTTCGCCAACAAGGGGTTAAACATCACCTGGCGGCAGGCATTGTATTAACCGGTGGCGCAGCGCAGATCGAAGGTCTTGCAGCCTGTGCTCAGCGCGTGTTTCATACGCAAGTGCGTATCGGCGCGCCGCTGAACATTACCGGTTTAACGGATTATGCTCAGGAGCCGTATTATTCGACGGCGGTGGGATTGCTTCACTATGGGAAAGAGTCACATCTTAACGGTGAAGCTGAAGTAGAAAAACGTGTTACAGCATCAG  >  NC_000913/104897‑105194

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