Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,154,926 A→T 5.8% N272I (AAC→ATC)  yceG → septation protein, ampicillin sensitivity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,154,9260AT 5.8% 307.8 / 3.7 95N272I (AAC→ATC) yceGseptation protein, ampicillin sensitivity
Reads supporting (aligned to +/- strand):  ref base A (48/41);  new base T (3/3);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAGGTGGTCACACGTTTAATACCAATCTTG  >  NC_000913/1154780‑1155072
                                                                                                                                                  |                                                                                                                                                  
tGAACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACATCGTATAaca                                                                                                                                                  >  5:22950/1‑149 (MQ=255)
tGAACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAaca                                                                                                                                                  <  8:161199/149‑1 (MQ=255)
tGAACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGACTGCAGACCGACCAGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAaca                                                                                                                                                  >  3:411192/1‑149 (MQ=255)
  aaCGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAacac                                                                                                                                                 <  3:248604/148‑1 (MQ=255)
  aaCGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACAcc                                                                                                                                                <  7:191188/149‑1 (MQ=255)
  aaCGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACAcc                                                                                                                                                >  3:396953/1‑149 (MQ=255)
       gATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATAc                                                                                                                                           <  3:87876/149‑1 (MQ=255)
        aTAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGTCAAACTTTCTCGTGCAGACCTGGAAAAGCAGACAGCGTATAACACCTATAc                                                                                                                                           >  5:99889/1‑148 (MQ=255)
                 ccTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGt                                                                                                                                 <  3:405162/149‑1 (MQ=255)
                  cTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTc                                                                                                                                >  3:123130/1‑149 (MQ=255)
                   tCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCt                                                                                                                               >  2:523331/1‑149 (MQ=255)
                   tCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAACGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCt                                                                                                                               >  7:71193/1‑149 (MQ=255)
                      gTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACTGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACAATTACCGTTCTgcc                                                                                                                            >  2:183405/1‑149 (MQ=255)
                          ttATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTgccgc                                                                                                                          >  7:164503/1‑147 (MQ=255)
                           tATGAACCGTTTACGTATTGGTATGCGGCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCa                                                                                                                        <  6:22950/148‑1 (MQ=255)
                            aTCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAg                                                                                                                       <  8:164503/148‑1 (MQ=255)
                            aTCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAg                                                                                                                       <  3:20343/148‑1 (MQ=255)
                             tCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCtata                                                                                                                                            <  2:559128/126‑1 (MQ=255)
                             tCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGt                                                                                                                     >  4:58697/1‑149 (MQ=255)
                               caCCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCt                                                                                                                                               >  7:398409/2‑121 (MQ=255)
                               aaCCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATATCACCTATACCATTACCGGTCTGCCGCCAGGTGc                                                                                                                   >  3:344635/1‑149 (MQ=255)
                               aaCCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCt                                                                                                                                               <  8:398409/121‑1 (MQ=255)
                               aaCCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCTAGGTGc                                                                                                                   >  7:112716/1‑149 (MQ=255)
                               aaCCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTAGATCACCTATAACATTACCGGTCTGCCGCCAGGTGc                                                                                                                   >  3:344638/1‑149 (MQ=255)
                               aaCCGTTTACGCATTGGTATGCGCCTGCAGACAGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGTCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATATCACCTATACCATTACCGGTCTGCCGCCAGGTGc                                                                                                                   >  7:501170/1‑149 (MQ=255)
                                  cGTTTACGCATTGGTATGCGCCTGCAGCCCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGa                                                                                                                 >  4:397357/1‑148 (MQ=255)
                                      tACGGATTGGTATGGGCCTGCAGACGGACGCGACCGTGATTTACGGGATGGGATAGCGTTATATTGGCATACTTTCTCGTGCAGACCTGGAAACGCGGACAGCGTATATCACCTATACCATTTCCGGTCTGCCGCGAGGTGCGATTGc                                                                                                             <  4:396953/148‑1 (MQ=255)
                                       aCGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGa                                                                                                           >  3:229208/1‑149 (MQ=255)
                                           aTTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTgc                                                                                                                             <  7:460359/127‑1 (MQ=255)
                                           aTTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTgc                                                                                                                             >  8:460359/1‑127 (MQ=255)
                                           aTTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGcc                                                                                                       >  8:16680/1‑149 (MQ=255)
                                            ttGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATCCCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCg                                                                                                      <  1:523331/149‑1 (MQ=255)
                                            ttGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCg                                                                                                      <  1:2104/149‑1 (MQ=255)
                                                    cgcCTGCAGACCGACCCGACCGTGATTTGCGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCtata                                                                                                                                            >  4:328325/1‑103 (MQ=255)
                                                    cgcCTGCAGACCGACCCGACCGTGATTTGCGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCtata                                                                                                                                            <  3:328325/103‑1 (MQ=255)
                                                           aGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCt                                                                                        >  8:290839/1‑148 (MQ=255)
                                                            gaccgaccCGCCCGAGAGTTACGTGAGGCGAGAGCGTTATAATGGCAAACTATCTCGTGCAGACCTCGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATACCGACGCCGGGGGCGGATTCGCTGa                                                                                      <  1:183405/149‑1 (MQ=255)
                                                            gaccgaccCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTg                                                                                       >  7:139223/1‑148 (MQ=255)
                                                             accgaccCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGaa                                                                                     <  7:290839/149‑1 (MQ=255)
                                                              ccgacccgacCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAg                                                                                    <  4:229208/149‑1 (MQ=255)
                                                                  cccgaccGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAg                                                                                    <  7:276863/145‑1 (MQ=255)
                                                                  cccgaccGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAg                                                                                    >  8:276863/1‑145 (MQ=255)
                                                                    cgaccGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTgc                                                                                                                             <  4:233465/102‑1 (MQ=255)
                                                                    cgaccGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTgc                                                                                                                             >  3:233465/1‑102 (MQ=255)
                                                                      accGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGa                                                                                      >  5:478429/1‑139 (MQ=255)
                                                                      accGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGa                                                                                      <  6:478429/139‑1 (MQ=255)
                                                                        cGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTgcg                                                                           >  1:24558/1‑148 (MQ=255)
                                                                        cGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTgcg                                                                           >  1:214610/1‑148 (MQ=255)
                                                                              ttACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCgg                                                                    >  3:75051/1‑149 (MQ=255)
                                                                                aCGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGc                                                                   <  5:201669/148‑1 (MQ=255)
                                                                                    gATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCaaaaa                                                              >  7:130432/1‑149 (MQ=255)
                                                                                    gATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCaaaaa                                                              <  2:24558/149‑1 (MQ=255)
                                                                                    gATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCaaaaa                                                              >  3:50256/1‑149 (MQ=255)
                                                                                      tGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCAGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACg                                                            >  3:195354/1‑149 (MQ=255)
                                                                                                ttATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTAtctc                                                   >  3:135394/1‑148 (MQ=255)
                                                                                                    aaTGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTAttt                                              >  7:349083/1‑149 (MQ=255)
                                                                                                     aTGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTAtttt                                             <  4:50256/149‑1 (MQ=255)
                                                                                                     aTGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGGCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTAtttt                                             <  4:411192/149‑1 (MQ=255)
                                                                                                      tGGCAAACTTTCTCGTGCAGACCTGGAAACTCCGACAGCGTATAACACCTATGCCATTACCGGTCTGCCGCCAGGTGAGATAGCGACGCCGGTGGCGGATTCGCTGAATGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATATtg                                            >  2:302664/1‑149 (MQ=255)
                                                                                                        gCAAACTTTCTCGTGCAGAACTGGAAACGCCGACAGCGTATAAAACCTATAACATTTCCGTTCTGCCGCCAGGTGCGATAGCGAAGCCGGGGGCGGATTCGCTGAAGGCTGATGCGCATCCGGCAAAAACGCCGTATCTCTATTTtgtg                                          <  8:501170/149‑1 (MQ=255)
                                                                                                         cAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTg                                                                                       >  8:564500/1‑103 (MQ=255)
                                                                                                         cAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTg                                                                                       <  7:564500/103‑1 (MQ=255)
                                                                                                         cAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTgg                                         <  7:16680/149‑1 (MQ=255)
                                                                                                         cAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTgg                                         >  5:140812/1‑149 (MQ=255)
                                                                                                          aaaCTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTtgt                                           >  6:151813/1‑146 (MQ=255)
                                                                                                          aaaCTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTtgt                                           <  5:151813/146‑1 (MQ=255)
                                                                                                          aaaCTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGc                                        >  1:235928/1‑149 (MQ=255)
                                                                                                           aaCTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGa                                                                                      <  3:320503/102‑1 (MQ=255)
                                                                                                           aaCTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGa                                                                                      >  4:320503/1‑102 (MQ=255)
                                                                                                            aCTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGgctgct                                                                              >  3:300709/1‑109 (MQ=255)
                                                                                                            aCTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGgctgct                                                                              <  4:300709/109‑1 (MQ=255)
                                                                                                             cTTTCTCGTGCAGACCTGGACACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGa                                     <  7:196944/149‑1 (MQ=255)
                                                                                                             cTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTa                                                 >  2:176024/1‑137 (MQ=255)
                                                                                                             cTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTa                                                 <  1:176024/137‑1 (MQ=255)
                                                                                                              tttCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGAt                                    >  1:292101/1‑149 (MQ=255)
                                                                                                                tctcGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATg                                   <  8:139223/148‑1 (MQ=255)
                                                                                                                tctcGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATg                                   <  4:135394/148‑1 (MQ=255)
                                                                                                                    gTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGt                                 <  3:58697/146‑1 (MQ=255)
                                                                                                                     tGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAg                             >  6:215686/1‑149 (MQ=255)
                                                                                                                     tGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAg                             <  8:408608/149‑1 (MQ=255)
                                                                                                                      gCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAgg                            >  1:220972/1‑149 (MQ=255)
                                                                                                                         gACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAggtgg                         <  2:214610/149‑1 (MQ=255)
                                                                                                                          aCCTGGAAACGCCGACAGCGTATATCACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAggtggt                        <  4:344635/149‑1 (MQ=255)
                                                                                                                          aCCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAggtggt                        >  5:322601/1‑149 (MQ=255)
                                                                                                                          aCCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAggtggt                        >  2:45832/1‑149 (MQ=255)
                                                                                                                           ccTGGAAACGCCGACAGCGTATATCACCTATACCATTACCGGGCTGCCGCCAGTTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAggtggt                        <  4:344638/148‑1 (MQ=255)
                                                                                                                            cTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAGGTGGTca                      >  6:498672/1‑149 (MQ=255)
                                                                                                                              ggAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTAt                                                      <  5:240257/115‑1 (MQ=255)
                                                                                                                              ggAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTAt                                                      >  6:240257/1‑115 (MQ=255)
                                                                                                                                     ccGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAGGTGGTCCCACGTTTaa             >  5:98435/1‑149 (MQ=255)
                                                                                                                                      cGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAGGTGGTCACACGTTTAAt            <  1:45832/149‑1 (MQ=255)
                                                                                                                                        aCAGCGTATAACACCTATACCATTACCGGTCTGCCGCAAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCAGATGGTAAAGGTCGTCACACGTTTAATAc          >  7:238844/1‑149 (MQ=255)
                                                                                                                                         cAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTtgtg                                          <  6:474202/116‑1 (MQ=255)
                                                                                                                                         cAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTtgtg                                          >  5:474202/1‑116 (MQ=255)
                                                                                                                                          aGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAGGTGGTCACACGTTTAATACCa        <  2:220972/149‑1 (MQ=255)
                                                                                                                                             gTATAACACCTATACCATTACCGGTCTGCCGCCTGGTGCGTTAGCGACGCGGGGGGCGGATTCGCTGAAGGGTGCTGCGCATCCGGCAATAACGCCGTATCTCTATTTTGTTGCCGATTGTAAAGGTGGTCACACGTTTAATACCATTc     <  8:238844/149‑1 (MQ=255)
                                                                                                                                              tataACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAGCTGGTCACACGTTTAATACCAATCt    >  5:227194/1‑149 (MQ=255)
                                                                                                                                                taACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAGGTGGTCACACGTTTAATACCAATCtt   <  5:498672/148‑1 (MQ=255)
                                                                                                                                                taACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAGGTGGTCACACGTTTAATACCAATCTTg  >  1:390233/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGAACGCGATAAGGTTGCCTCAGTATTTATCAACCGTTTACGCATTGGTATGCGCCTGCAGACCGACCCGACCGTGATTTACGGGATGGGAGAGCGTTATAATGGCAAACTTTCTCGTGCAGACCTGGAAACGCCGACAGCGTATAACACCTATACCATTACCGGTCTGCCGCCAGGTGCGATAGCGACGCCGGGGGCGGATTCGCTGAAGGCTGCTGCGCATCCGGCAAAAACGCCGTATCTCTATTTTGTGGCCGATGGTAAAGGTGGTCACACGTTTAATACCAATCTTG  >  NC_000913/1154780‑1155072

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