Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,400,966 A→C I41S (ATT→AGT)  mreB ← cell wall structural complex MreBCD, actin‑like component MreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,400,9660AC98.8% 303.6 / ‑3.5 80I41S (ATT→AGT) mreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (37/42);  total (37/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01

CATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGG  >  NC_000913/3400818‑3401112
                                                                                                                                                    |                                                                                                                                                  
cATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTg                                                                                                                                                  <  1:845339/151‑1 (MQ=255)
 aTTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTgg                                                                                                                                                 >  2:226863/1‑151 (MQ=255)
    ttttCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGcc                                                                                                                                               >  2:1285322/1‑150 (MQ=255)
    ttttCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGcc                                                                                                                                               >  1:76873/1‑150 (MQ=255)
    ttttCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGcc                                                                                                                                               >  1:1239625/1‑150 (MQ=255)
     tttCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGccac                                                                                                                                             >  1:1082496/1‑151 (MQ=255)
     tttCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGcca                                                                                                                                              <  1:751399/150‑1 (MQ=255)
       tCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGccacc                                                                                                                                            >  2:393654/1‑150 (MQ=255)
         aGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACg                                                                                                                                         >  2:536747/1‑151 (MQ=255)
          gTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACgg                                                                                                                                        >  2:628303/1‑151 (MQ=255)
            cACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGaa                                                                                                                                      <  2:666884/151‑1 (MQ=255)
            cACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGa                                                                                                                                       >  1:298180/1‑150 (MQ=255)
              cGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAgg                                                                                                                                    >  2:293181/1‑151 (MQ=255)
               gaagaaGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGc                                                                                                                                   <  2:9856/151‑1 (MQ=255)
                 agaagTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTc                                                                                                                                 >  1:816036/1‑151 (MQ=255)
                 agaagTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTc                                                                                                                                 >  2:556185/1‑151 (MQ=255)
                    agTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCAtt                                                                                                                              <  1:101162/151‑1 (MQ=255)
                    agTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCAtt                                                                                                                              <  2:317559/151‑1 (MQ=255)
                        ggCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAAt                                                                                                                          <  2:980252/151‑1 (MQ=255)
                         gCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTTAAt                                                                                                                          >  1:143352/1‑150 (MQ=255)
                         gCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAAt                                                                                                                          >  1:1053368/1‑150 (MQ=255)
                              aaCGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATg                                                                                                                    >  2:324179/1‑151 (MQ=255)
                              aaCGCCGTCGTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCTCGGAAGGCTCATTCAATACGATg                                                                                                                    <  1:1186666/151‑1 (MQ=255)
                                 gCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCt                                                                                                                 <  2:212351/151‑1 (MQ=255)
                                 accGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGcc                                                                                                                  >  1:667783/2‑150 (MQ=255)
                                         tCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCtt                                                                                                          <  2:562875/150‑1 (MQ=255)
                                             tGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTAc                                                                                                      <  1:1175927/150‑1 (MQ=255)
                                              gggCGAATGGCAGCAATATTGCCCGGCCTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTGTCGGTGAACCGGCACGAGCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACAt                                                                                                    <  1:680942/151‑1 (MQ=255)
                                                   aaTGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATg                                                                                               >  2:1154171/1‑151 (MQ=255)
                                                    aTGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATg                                                                                               >  2:967600/1‑150 (MQ=255)
                                                     tGGCAGCAATATTGCCCGGCGTACGGCCCAGCACCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGAACCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGa                                                                                              <  2:170371/150‑1 (MQ=255)
                                                        cagcaATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGt                                                                                          >  2:969111/1‑151 (MQ=255)
                                                             atatTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCg                                                                                     >  1:598315/1‑151 (MQ=255)
                                                               atTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCa                                                                                   <  2:421272/151‑1 (MQ=255)
                                                                 tGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACa                                                                                                     <  1:248700/131‑1 (MQ=255)
                                                                 tGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACa                                                                                                     >  2:248702/1‑131 (MQ=255)
                                                                 tGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAg                                                                                  <  2:200036/150‑1 (MQ=255)
                                                                  gCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGt                                                                                 >  1:183500/1‑150 (MQ=255)
                                                                      ggCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTAccc                                                                             >  1:797519/1‑150 (MQ=255)
                                                                      ggCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCa                                                                            <  2:716659/151‑1 (MQ=255)
                                                                       gCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTAc                                                                               <  1:611731/147‑1 (MQ=255)
                                                                        cGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAgg                                                                          >  2:1245365/1‑151 (MQ=255)
                                                                         gTACCGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCCTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGt                                                                         >  1:897336/1‑151 (MQ=255)
                                                                              gCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAAt                                                                     >  1:236736/1‑150 (MQ=255)
                                                                                ccAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATgg                                                                   >  1:1300976/1‑150 (MQ=255)
                                                                                   gCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGACGGTATTCGCAGTACCCAGGTCAATGGACa                                                                >  2:706887/1‑150 (MQ=255)
                                                                                         gCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCa                                                           <  1:191787/149‑1 (MQ=255)
                                                                                         gCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCAt                                                          <  2:831600/150‑1 (MQ=255)
                                                                                            tCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTg                                                        <  1:104405/149‑1 (MQ=255)
                                                                                             cgcgTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGaa                                                     <  2:1205455/151‑1 (MQ=255)
                                                                                              gcgTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGaaa                                                    <  1:1177968/151‑1 (MQ=255)
                                                                                                   aTGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATg                                               <  1:869333/151‑1 (MQ=255)
                                                                                                   aTGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATg                                               <  1:706214/151‑1 (MQ=255)
                                                                                                   aTGACATACTGCATCTACGCTTTTAGGTGAACCGGCACTATCCTGTCTACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGACGGTATTCGCAGTCCCCAGGACAATGGACCAGTCATGGCAAACCATg                                               <  1:176221/151‑1 (MQ=255)
                                                                                                      aCCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCa                                            <  2:729228/151‑1 (MQ=255)
                                                                                                       ccTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCAc                                           >  2:808152/1‑151 (MQ=255)
                                                                                                       ccTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCAc                                           <  2:1144187/151‑1 (MQ=255)
                                                                                                       ccTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCAc                                           <  2:328379/151‑1 (MQ=255)
                                                                                                        cTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACg                                          >  2:177156/1‑151 (MQ=255)
                                                                                                          aCTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACg                                          <  1:117308/149‑1 (MQ=255)
                                                                                                          aCTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACg                                          <  1:134239/149‑1 (MQ=255)
                                                                                                             gCATCTACGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAAtt                                     <  1:579710/151‑1 (MQ=255)
                                                                                                                   aCGCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTc                                <  1:134558/150‑1 (MQ=255)
                                                                                                                     gCTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAAc                             >  1:655208/1‑151 (MQ=255)
                                                                                                                      cTTTTCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAAc                             >  2:205026/1‑150 (MQ=255)
                                                                                                                       ttttCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACAt                           >  2:611520/1‑151 (MQ=255)
                                                                                                                        tttCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACa                            <  1:1035184/149‑1 (MQ=255)
                                                                                                                        tttCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACAt                           >  2:291873/1‑150 (MQ=255)
                                                                                                                          tCGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATAc                         <  1:331429/150‑1 (MQ=255)
                                                                                                                           cGGTGAACCGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACt                        >  1:1281485/1‑150 (MQ=255)
                                                                                                                                  ccGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCATTACCCAGGTCAATGGACAAGTCATTGGAAACCATGCCACGAAATTTTTTCAACATACTAAGGGAt                 <  1:362225/150‑1 (MQ=255)
                                                                                                                                  ccGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATa                <  1:193467/151‑1 (MQ=255)
                                                                                                                                   cGGCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATaa               <  1:1138911/151‑1 (MQ=255)
                                                                                                                                    ggCACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAAt              <  2:625181/151‑1 (MQ=255)
                                                                                                                                      cACGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATc             <  1:369411/150‑1 (MQ=255)
                                                                                                                                        cGATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTg          <  2:476217/151‑1 (MQ=255)
                                                                                                                                         gATCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTg          >  1:556483/1‑150 (MQ=255)
                                                                                                                                           tCCTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaa        <  2:1208020/150‑1 (MQ=255)
                                                                                                                                             cTGACGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGc     >  1:1086564/1‑151 (MQ=255)
                                                                                                                                                aCGACTGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTgg  <  2:90830/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATTTTTTCAGTCACGAAGAAGTCGGCGATAACGCCGTCTTTCATTGGGCGAATGGCAGCAATATTGCCCGGCGTACGGCCCAGCATCTGCTTCGCGTCATGACCTACTGCAGCTACGCTTTTCGGTGAACCGGCACGATCCTGACGAATGGCCACCACGGAAGGCTCATTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGG  >  NC_000913/3400818‑3401112

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