Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,786,899 T→A 6.7% L21H (CTT→CAT) ‡ envC → activator of AmiB,C murein hydrolases, septal ring factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,786,8990TA 6.7% 279.9 / 7.3 101L21H (CTT→CAT) ‡envCactivator of AmiB,C murein hydrolases, septal ring factor
Reads supporting (aligned to +/- strand):  ref base T (46/48);  new base A (5/2);  total (51/50)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

Only 100 of 103 total aligned reads displayed.

GGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAA  >  NC_000913/3786751‑3787045
                                                                                                                                                    |                                                                                                                                                  
ggTATGGAACTCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCtt                                                                                                                                                   <  1:1059015/150‑1 (MQ=255)
ggTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCtt                                                                                                                                                   <  2:1455114/150‑1 (MQ=255)
    tGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAgcg                                                                                                                                               >  2:214660/1‑150 (MQ=255)
    tGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAgcg                                                                                                                                               >  2:583281/1‑150 (MQ=255)
     ggAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCAAGAGGGGAATGGCGATTAATACAAAGACACGGGACGTGAAACCGCGCAGCTTTGCAATCAGGCCCCTCATCTACGCCCGCCTTCTTAgcgg                                                                                                                                              >  1:228699/1‑149 (MQ=255)
      gAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTAATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGAAGGTTTGAAATCAGGCCCATCATCTACGCCAGCTTTCTTATCGCt                                                                                                                                             >  1:1363601/1‑150 (MQ=255)
       aaaTCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTg                                                                                                                                            >  2:1527438/1‑150 (MQ=255)
       aaaTCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTg                                                                                                                                            >  1:1602122/1‑150 (MQ=255)
       aaaTCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTg                                                                                                                                            >  2:380508/1‑150 (MQ=255)
          tCCAGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACCCGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCg                                                                                                                                         <  1:138563/150‑1 (MQ=255)
           cccGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGt                                                                                                                                        <  1:1432879/150‑1 (MQ=255)
            ccGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTa                                                                                                                                       <  1:996801/150‑1 (MQ=255)
             cGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGCTAAGAACGCTGGCGTag                                                                                                                                      <  2:733069/150‑2 (MQ=255)
              gCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCCTGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAACCCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTAtt                                                                                                                                     >  2:260026/1‑150 (MQ=255)
                  gagaTGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTAttgttg                                                                                                                                 >  1:419843/1‑150 (MQ=255)
                         cTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGAGCCt                                                                                                                          >  1:1043883/1‑150 (MQ=255)
                              aaGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGGTTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTcg                                                                                                                     >  1:991200/1‑149 (MQ=255)
                              aaGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCCGCGTTGTTAGCGCTGGCGTATTGTTGTGCGCCTTTTcc                                                                                                                     >  2:538836/1‑150 (MQ=255)
                               aGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCCTCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCg                                                                                                                    <  2:1078025/150‑1 (MQ=255)
                                gCCGCTGTTCATCGGGGAATCAGCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGc                                                                                                                   <  1:9778/150‑1 (MQ=255)
                                 ccGCTGTTCCTCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCCATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGcc                                                                                                                  >  1:854620/1‑150 (MQ=255)
                                          aTCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGAt                                                                                                         <  2:805394/150‑1 (MQ=255)
                                          aTCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGAt                                                                                                         <  1:1538745/150‑1 (MQ=255)
                                           tCGTGGAATACTCCCTCCCCCTGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGTAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATg                                                                                                        <  1:1666322/150‑1 (MQ=255)
                                           tCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATg                                                                                                        <  1:1623141/150‑1 (MQ=255)
                                             gTGGAATACTCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGCCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAg                                                                                                      <  1:486891/150‑1 (MQ=255)
                                             gTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAg                                                                                                      >  1:189661/1‑150 (MQ=255)
                                               ggAATAATCCCTCCCCCTGAGGGGAAAGGCGAATAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCg                                                                                                    <  1:1563366/150‑1 (MQ=255)
                                               ggAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGCTAAGAACGCTGGCGTagatgat                                                                                                                                 <  2:1529595/121‑7 (MQ=38)
                                               ggAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGCTAAGAACGCTGGCGTagatgat                                                                                                                                 >  1:1529595/1‑115 (MQ=38)
                                                 aataatCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCTATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTg                                                                                                  <  1:1231714/150‑1 (MQ=255)
                                                   taatCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGAc                                                                                                >  2:352910/1‑150 (MQ=255)
                                                    aatCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGAcc                                                                                               <  2:885027/150‑1 (MQ=255)
                                                     atCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCa                                                                                              >  2:894991/1‑150 (MQ=255)
                                                     atCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCa                                                                                              <  1:1302818/150‑1 (MQ=255)
                                                     atCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGGGAAACCGCGCAGGGTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCa                                                                                              <  1:791658/150‑1 (MQ=255)
                                                      tccctcccCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCaa                                                                                             >  2:288915/1‑150 (MQ=255)
                                                      tccctcccCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCaa                                                                                             >  2:1180122/1‑150 (MQ=255)
                                                      tccctcccCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCaa                                                                                             >  2:1522040/1‑150 (MQ=255)
                                                         ctcccCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTc                                                                                          <  1:288915/150‑1 (MQ=255)
                                                         ctcccCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTc                                                                                          <  1:1522040/150‑1 (MQ=255)
                                                         ctcccCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTc                                                                                          <  1:352910/150‑1 (MQ=255)
                                                            cccATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCaaa                                                                                       >  1:1408144/1‑150 (MQ=255)
                                                            cccATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCaaa                                                                                       <  1:858278/150‑1 (MQ=255)
                                                            cccATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCaaa                                                                                       <  2:1626588/150‑1 (MQ=255)
                                                             ccATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTACTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCaaaa                                                                                      >  1:358431/1‑149 (MQ=255)
                                                              cctGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATc                                                                                     <  2:1082118/148‑1 (MQ=255)
                                                              cATGAGGGGAAAGGCGATTAAAACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATc                                                                                     <  1:1542724/150‑1 (MQ=255)
                                                                tGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTa                                                                                   <  1:1090204/150‑1 (MQ=255)
                                                                   ggggAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTc                                                                                >  2:1587012/1‑150 (MQ=255)
                                                                   ggggAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTGCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTc                                                                                <  1:287433/150‑1 (MQ=255)
                                                                     ggAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGCCCACCTCAAATCTATTCAg                                                                              >  1:1477323/1‑150 (MQ=255)
                                                                     ggAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAg                                                                              <  2:834511/150‑1 (MQ=255)
                                                                      gAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAgg                                                                             <  2:750148/150‑1 (MQ=255)
                                                                      gAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAgg                                                                             <  1:1587012/150‑1 (MQ=255)
                                                                          ggCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGa                                                                         <  2:1671569/150‑1 (MQ=255)
                                                                            cGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGCCCAACTCAAATCTATTCAGGCTGACa                                                                       >  1:1565493/1‑150 (MQ=255)
                                                                              aTTAATTCCATGACAAGGGCCGTGAAACCGCGCAGGTTTGCAATCAGTCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTTCGCCTTTTCCGCCCACGCGGATGAGCGTGACAAACTCAAATCTATACAGGCTGACaaa                                                                     >  1:312691/1‑148 (MQ=255)
                                                                              aTTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCGGACATc                                                                     <  1:1129028/150‑1 (MQ=255)
                                                                               ttAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGTCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATcg                                                                    <  1:1024229/150‑1 (MQ=255)
                                                                                     ccATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCaa                                                              >  2:1508988/1‑150 (MQ=255)
                                                                                      cATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCAATCATCTACGCCTGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGCGGGGGCCCCACTGAAATATGTTACGGCTGACTTCGCCGCaaa                                                             >  1:1088871/1‑150 (MQ=255)
                                                                                        tGACACTGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATAATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACAACGCCGCAAAag                                                           <  2:921354/150‑1 (MQ=255)
                                                                                        tGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAag                                                           <  1:1731844/150‑1 (MQ=255)
                                                                                        tGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAag                                                           <  1:67417/150‑1 (MQ=255)
                                                                                         gACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCCACTCAAATCTATTCAGGCTGACATCGCCGCAAAaga                                                          >  1:1269826/1‑150 (MQ=255)
                                                                                         gACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAaga                                                          >  1:737029/1‑150 (MQ=255)
                                                                                          acacGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAagag                                                         <  2:1168887/150‑1 (MQ=255)
                                                                                          acacGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAagag                                                         <  1:22016/150‑1 (MQ=255)
                                                                                            acGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCGTCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGAGTGACAAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAgcg                                                       <  2:855089/150‑1 (MQ=255)
                                                                                             cGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGCCCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAgcgc                                                      >  2:721324/1‑150 (MQ=255)
                                                                                               ggCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAgcgcgc                                                    <  1:244810/150‑1 (MQ=255)
                                                                                                  cGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCggg                                                 >  1:1315790/1‑149 (MQ=255)
                                                                                                  cGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGt                                                 >  1:946078/1‑150 (MQ=255)
                                                                                                  cGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGt                                                 >  1:627031/1‑150 (MQ=255)
                                                                                                  cGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGt                                                 >  1:643099/1‑150 (MQ=255)
                                                                                                    tGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTAc                                               >  2:764388/1‑150 (MQ=255)
                                                                                                     gAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACg                                              <  1:380508/150‑1 (MQ=255)
                                                                                                      aaaCCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGc                                             <  1:764026/150‑1 (MQ=255)
                                                                                                        aCCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCa                                           <  1:1194910/150‑1 (MQ=255)
                                                                                                          cgcgCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCACACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCTCAAAAGAGCGCGCGGTACGCCaaa                                         <  1:931766/150‑1 (MQ=255)
                                                                                                             gcAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGc                                      <  1:1527438/150‑1 (MQ=255)
                                                                                                              cAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGca                                     >  1:1306518/1‑150 (MQ=255)
                                                                                                                ggTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGcaac                                   >  1:1413748/1‑150 (MQ=255)
                                                                                                                  tttGCAATCAGGCCCATCATCTACGCCAGCGGTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTCCGCCTAAAGcaacaa                                 >  1:919008/1‑150 (MQ=255)
                                                                                                                   ttGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGcaacaac                                <  1:1503527/150‑1 (MQ=255)
                                                                                                                      cAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCaacaacaac                             >  1:583073/1‑150 (MQ=255)
                                                                                                                        aTCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAAcgc                           <  1:214660/150‑1 (MQ=255)
                                                                                                                              cccATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGCGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGc                     >  1:1103221/1‑150 (MQ=255)
                                                                                                                                    caatACGCCAGCGCTAAGAACGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATACAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTc               >  2:1500425/4‑150 (MQ=255)
                                                                                                                                     tcTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCg              <  2:643099/150‑1 (MQ=255)
                                                                                                                                     tcTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCg              <  1:1509548/150‑1 (MQ=255)
                                                                                                                                     aatACGCCAGCGCTAAGAACGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGCGCGCGCGGTACGCCAAAAGCGACAACAACGCGCAAGCGTGCTCg              >  1:307529/3‑150 (MQ=255)
                                                                                                                                      atACGCCAGCGCTAAGACCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTTGc             >  2:335622/2‑150 (MQ=255)
                                                                                                                                      atACGCCAGCGCTAAGAACGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGc             >  2:217626/2‑150 (MQ=255)
                                                                                                                                           ccAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAAt        >  1:1601060/1‑150 (MQ=255)
                                                                                                                                             aGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGATCGCACAATTg      <  2:1613545/150‑1 (MQ=255)
                                                                                                                                              gCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGa     >  1:988166/1‑150 (MQ=255)
                                                                                                                                                 ttCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGaaaa  <  1:1418739/150‑1 (MQ=255)
                                                                                                                                                 ttCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGAACGCAAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTCaaaa  >  2:1715520/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGTATGGAAATCCCGCAAGAGATGACTGGTAAGCCGCTGTTCATCGTGGAATAATCCCTCCCCATGAGGGGAAAGGCGATTAATACCATGACACGGGCCGTGAAACCGCGCAGGTTTGCAATCAGGCCCATCATCTACGCCAGCGTTCTTAGCGCTGGCGTATTGTTGTGCGCCTTTTCCGCCCACGCGGATGAGCGTGACCAACTCAAATCTATTCAGGCTGACATCGCCGCAAAAGAGCGCGCGGTACGCCAAAAGCAACAACAACGCGCAAGCCTGCTCGCACAATTGAAAA  >  NC_000913/3786751‑3787045

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