Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,020,958 A→T 5.6% intergenic (‑39/‑180) sulA ← / → sxy SOS cell division inhibitor/CRP‑S‑dependent promoter expression factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,020,9580AT 5.6% 397.7 / 10.5 107intergenic (‑39/‑180)sulA/sxySOS cell division inhibitor/CRP‑S‑dependent promoter expression factor
Reads supporting (aligned to +/- strand):  ref base A (44/57);  new base T (3/3);  total (47/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

Only 100 of 116 total aligned reads displayed.

TGAACCCATTCCCGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTT  >  NC_000913/1020684‑1021110
                                                                                                                                                                                                                                                                                  |                                                                                                                                                        
tGAACCCATTCCCGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACCCCAAATTTTACTTGCTGCGGAAGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCCATCCAGCCCCTTTGATTTACTATATGGATgtaca                                                                                                                                                          >  4:284367/1‑275 (MQ=255)
           ccGTCTCAGTGTGTGTGGGGGTGTTAACCAGAGTGGCCATCGCGGTTGCTGACGGAGTTGCTGTAAAAATGGGAACAGAATATGGTTCGTAAACATGGGATGATCTTCGCGAAAGACAACTGCACTGATAAGCCCGGCTGTAGTGTTTTCAGAAGAGACACGGGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCGGTGAGGTACTGTATGGATGTACAGTACATCCAGt                                                                                                                                               <  4:460073/275‑1 (MQ=255)
           ccGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGt                                                                                                                                               <  1:9947/275‑1 (MQ=255)
                                               ccAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCt                                                                                                           >  2:378908/1‑275 (MQ=255)
                                                                                    aaCAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAtt                                                                      >  4:221682/1‑275 (MQ=255)
                                                                                     aCAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAtt                                                                      <  2:246215/274‑1 (MQ=255)
                                                                                         tAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCt                                                                 <  2:181850/275‑1 (MQ=255)
                                                                                                                       gggaaaGACAACTTCACTGTTAAGCCCGGCGGGAGTGTTTTCCGTAGAGACACGGGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCGGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCg                                   <  4:174793/270‑1 (MQ=255)
                                                                                                                              aCAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGTCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTaa                            <  2:505718/275‑1 (MQ=255)
                                                                                                                                 aCTTCACTGTTAAGCCCGGCTGTAGGGTTTTCCGTATAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACgg                         <  4:361896/275‑1 (MQ=255)
                                                                                                                                 aCTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACgg                         <  2:221683/275‑1 (MQ=255)
                                                                                                                                          aTAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAg                                                                                                                                                <  7:789693/147‑1 (MQ=255)
                                                                                                                                           tAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTg                                                                                                                                              >  8:371394/1‑148 (MQ=255)
                                                                                                                                           tAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTg                                                                                                                                              <  7:371394/148‑1 (MQ=255)
                                                                                                                                              gCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAAtgtg            <  2:284368/275‑1 (MQ=255)
                                                                                                                                                ccGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATc                                                                                                                                                   <  8:1125713/138‑1 (MQ=255)
                                                                                                                                                ccGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATc                                                                                                                                                   >  7:1125713/1‑138 (MQ=255)
                                                                                                                                                 cGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCat                                                                                                                                                >  6:402079/1‑139 (MQ=255)
                                                                                                                                                   gCTGTAGTGTTTTCCGTAGAGACACGCGCGATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGCGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCa       <  2:139666/275‑1 (MQ=255)
                                                                                                                                                   gCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGa                                                                                                                                             >  6:166740/1‑141 (MQ=255)
                                                                                                                                                   gCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGa                                                                                                                                             <  5:166740/141‑1 (MQ=255)
                                                                                                                                                     tGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGa                                                                                                                                    <  5:552455/148‑1 (MQ=255)
                                                                                                                                                        aGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagtaca                                                                                                                                                     <  7:676872/128‑1 (MQ=255)
                                                                                                                                                        aGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagtaca                                                                                                                                                     >  8:676872/1‑128 (MQ=255)
                                                                                                                                                         gtgtTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCtt  <  4:130917/274‑1 (MQ=255)
                                                                                                                                                           gtTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGaca                                                                                                                                           >  8:117773/1‑135 (MQ=255)
                                                                                                                                                           gtTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGaca                                                                                                                                           <  7:117773/135‑1 (MQ=255)
                                                                                                                                                                 cGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCa                                                                                                                                                 <  8:462656/123‑1 (MQ=255)
                                                                                                                                                                     gagaCACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCgg                                                                                                                    >  6:544887/1‑148 (MQ=255)
                                                                                                                                                                           cgcgCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCata                                                                                                                                               >  8:172498/1‑113 (MQ=255)
                                                                                                                                                                           cgcgCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAagag                                                                                                              >  7:763115/1‑148 (MQ=255)
                                                                                                                                                                              gcAATTTTACTTGTTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATc                                                                                                                                                   >  6:405574/1‑108 (MQ=255)
                                                                                                                                                                              gcAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATc                                                                                                                                                   <  5:405574/108‑1 (MQ=255)
                                                                                                                                                                               cAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGaa                                                                                                                  >  7:273000/1‑140 (MQ=255)
                                                                                                                                                                               cAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGaa                                                                                                                  <  8:273000/140‑1 (MQ=255)
                                                                                                                                                                                          gctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGt                                                                                               <  5:544887/148‑1 (MQ=255)
                                                                                                                                                                                                  tGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCg                                                                                                                     <  5:617724/118‑1 (MQ=255)
                                                                                                                                                                                                  tGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCg                                                                                                                     >  6:617724/1‑118 (MQ=255)
                                                                                                                                                                                                      aaCGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGAc                                                                                     >  7:887850/1‑146 (MQ=255)
                                                                                                                                                                                                      aaCGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGAc                                                                                     <  8:887850/146‑1 (MQ=255)
                                                                                                                                                                                                       acgacgAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCa                                                                                                                                                 <  7:89057/85‑1 (MQ=255)
                                                                                                                                                                                                       acgacgAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCa                                                                                                                                                 >  8:89057/1‑85 (MQ=255)
                                                                                                                                                                                                           cgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTa                             >  4:570440/1‑197 (MQ=255)
                                                                                                                                                                                                               gaaCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACgg                         <  2:570441/197‑1 (MQ=255)
                                                                                                                                                                                                                aaCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATaa                                                                         >  8:905636/1‑148 (MQ=255)
                                                                                                                                                                                                                    aTGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCATTGACAACAAAGATCAACCCTATTTTCGGAAag                                                                                                                >  7:687513/1‑105 (MQ=255)
                                                                                                                                                                                                                    aTGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAag                                                                                                                <  8:687513/105‑1 (MQ=255)
                                                                                                                                                                                                                    aTGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTa                                                                     >  7:315721/1‑148 (MQ=255)
                                                                                                                                                                                                                     tgtgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTaa                                                                    <  8:763115/148‑1 (MQ=255)
                                                                                                                                                                                                                       tgCATATCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATGTTCGTAGTTAACGGATCNGGTAATGTGa           >  4:558008/1‑203 (MQ=255)
                                                                                                                                                                                                                          aTAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCa                                                                                                                                                 <  6:398644/66‑1 (MQ=255)
                                                                                                                                                                                                                          aTAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCa                                                                                                                                                 >  5:398644/1‑66 (MQ=255)
                                                                                                                                                                                                                          aTAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTg                                                               <  8:315721/148‑1 (MQ=255)
                                                                                                                                                                                                                           tatccTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCa       <  2:558009/200‑1 (MQ=255)
                                                                                                                                                                                                                           tAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAt                                                                                                    >  8:300787/1‑110 (MQ=255)
                                                                                                                                                                                                                           tAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAt                                                                                                    <  7:300787/110‑1 (MQ=255)
                                                                                                                                                                                                                           tAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGtt                                                                                           >  8:626811/1‑119 (MQ=255)
                                                                                                                                                                                                                           tAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGtt                                                                                           <  7:626811/119‑1 (MQ=255)
                                                                                                                                                                                                                             gCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCa                                                                                                                                                 <  7:912844/63‑1 (MQ=255)
                                                                                                                                                                                                                             gCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCa                                                                                                                                                 >  8:912844/1‑63 (MQ=255)
                                                                                                                                                                                                                               cTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTg                                                               >  7:919126/1‑143 (MQ=255)
                                                                                                                                                                                                                                tGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCtt                                                         >  8:756236/1‑148 (MQ=255)
                                                                                                                                                                                                                                      gtACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAttt                                                                                                  <  8:850986/101‑1 (MQ=255)
                                                                                                                                                                                                                                      gtACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAttt                                                                                                  >  7:850986/1‑101 (MQ=255)
                                                                                                                                                                                                                                       tACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGa           >  2:44857/1‑187 (MQ=255)
                                                                                                                                                                                                                                           tAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCa       <  4:44856/187‑1 (MQ=255)
                                                                                                                                                                                                                                            aatcaatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACGTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTa                                                        <  5:195182/137‑1 (MQ=255)
                                                                                                                                                                                                                                              tcaatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCt                                                                 <  5:451818/126‑1 (MQ=255)
                                                                                                                                                                                                                                              tcaatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCt                                                                 >  6:451818/1‑126 (MQ=255)
                                                                                                                                                                                                                                                   ccAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAGGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTAttt                                      >  8:7937/1‑148 (MQ=255)
                                                                                                                                                                                                                                                      gCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAt                                                                       >  8:463548/1‑112 (MQ=255)
                                                                                                                                                                                                                                                      gCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAt                                                                       <  7:463548/112‑1 (MQ=255)
                                                                                                                                                                                                                                                       ccccTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTg                                                                                       <  7:258155/95‑1 (MQ=255)
                                                                                                                                                                                                                                                         ccTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAag                                                                                                                <  5:138372/68‑1 (MQ=255)
                                                                                                                                                                                                                                                         ccTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAag                                                                                                                >  6:138372/1‑68 (MQ=255)
                                                                                                                                                                                                                                                         ccTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTaa                                                                    >  6:785421/1‑112 (MQ=255)
                                                                                                                                                                                                                                                         ccTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTaa                                                                    <  5:785421/112‑1 (MQ=255)
                                                                                                                                                                                                                                                            gtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAGGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTa                             <  7:7937/148‑1 (MQ=255)
                                                                                                                                                                                                                                                            gtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                                                                                >  7:55272/1‑81 (MQ=255)
                                                                                                                                                                                                                                                            gtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                                                                                <  8:55272/81‑1 (MQ=255)
                                                                                                                                                                                                                                                            gtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTa                             >  2:260411/1‑148 (MQ=255)
                                                                                                                                                                                                                                                            gtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTa                             >  7:248920/1‑148 (MQ=255)
                                                                                                                                                                                                                                                            gtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGa           >  2:383854/1‑166 (MQ=255)
                                                                                                                                                                                                                                                                gTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACgg                         <  4:260410/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                gTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCa       <  4:383853/166‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAAc                           >  6:163957/1‑143 (MQ=255)
                                                                                                                                                                                                                                                                  tACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAAc                           <  5:163957/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGt                                               <  8:1124157/122‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGt                                               >  7:1124157/1‑122 (MQ=255)
                                                                                                                                                                                                                                                                        aTGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGt                                                                                               <  5:528927/70‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aTGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGaa                                                                                >  6:568747/1‑85 (MQ=255)
                                                                                                                                                                                                                                                                          ggATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAAt               <  8:248920/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGaaa                                                                               >  5:397713/1‑83 (MQ=255)
                                                                                                                                                                                                                                                                           gATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGaaa                                                                               <  6:397713/83‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATa                                                  <  5:116552/112‑1 (MQ=255)
                                                                                                                                                                                                                                                                           gATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATa                                                  >  6:116552/1‑112 (MQ=255)
                                                                                                                                                                                                                                                                              gtacagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAt                                                                                                    <  6:515823/59‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gtacagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAt                                                                                                    >  5:515823/1‑59 (MQ=255)
                                                                                                                                                                                                                                                                               tacagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTAtt                                       >  4:565995/1‑119 (MQ=255)
                                                                                                                                                                                                                                                                                acagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTa                                                        >  7:50037/1‑101 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                        
TGAACCCATTCCCGACTCAGTTTTTGTTGCGGTGTTAACCAGAGTTGCCAGCGCGATTGCTGACCGAGTTGCTGTAACAATGGCAACAGTAGAAGTTGCGTCATCATGGGCTGATCTTCGCGATAGACAACTTCACTGATAAGCCCGGCTGTAGTGTTTTCCGTAGAGACACGCGCAATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTTAAGAATAAGTTGCCTATTTTCGTAGTTAACGGATCCGTTAATGTGAATCATTCTT  >  NC_000913/1020684‑1021110

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