Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,653,367 G→A 21.3% intergenic (+28/‑488) sseA → / → ryfA 3‑mercaptopyruvate sulfurtransferase/novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,653,3670GA21.3% 133.8 / 33.5 66intergenic (+28/‑488)sseA/ryfA3‑mercaptopyruvate sulfurtransferase/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base G (22/30);  new base A (7/7);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-02

GTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATT  >  NC_000913/2653102‑2653577
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                  
gTGGACGGAACTGGTGCGCGCAGGCGAACTAAAAACGACAGATGAACTGGCTGCGATATTTTGTCGTCGCGGCGTCATCTACGACAAACCACTTATCATCAGCTGCGGCGCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAAGGTGAAACTGTACGACGGCGCATGGA‑‑GTGCATGGGGGGCGCGGGCGGTATTACCGGTTGAGCCCGTGAAATAGGTATGTTACGGGAAATAAAAAACCGCCGCATTAGGct                                                                                                                                                                                                           >  4:368306/1‑274 (MQ=255)
      ggAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGGTTTACCGGTTGAGCCGGTGAAATCAGTATTTTACAGGCAATAAAAAAACGCCGAATTTGGCGGtttt                                                                                                                                                                                                      >  2:102896/1‑274 (MQ=255)
      ggAAATGGTGCAATAAGCGGAAATAAAAACCGCCGATGAACTGGGTGCGATATTTTTGTGTCTTGGGGTCAGCTACGACAAACCAATTATCGTGCACTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGGGACGCTGTATGTGCCAAACGTGATACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAACCAGTGAAATAAGTGTTTTACAGGCAATAAAAGACCGCCGAATTTGGCGGttttt                                                                                                                                                                                                     <  4:102895/275‑1 (MQ=255)
           tggtgacgtaaGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTg                                                                                                                                                                                                <  2:430478/266‑1 (MQ=255)
             gTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATACGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGc                                                                                                                                                                                               >  1:40102/1‑274 (MQ=255)
                  cgAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGCTTTCACGGTTGAGCAGCGGAAATAAGGTTTTTAACGGCATTAAAAAACCGCCGTATTTGGCGGGTTTTTAtttaatgtg                                                                                                                                                                                          >  4:455299/1‑267 (MQ=255)
                                      ccGTTGAACTGGATGCGTTATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTTCGGCTCTGGTGTAACGGCATCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCACATGGA‑‑ATGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCcagc                                                                                                                                                                      <  2:163059/274‑1 (MQ=255)
                                            aaCTGGATGCATTATTTTTTGGTCGCTGCGTCAGCTACGACAAACCAATTATCGTCAGTTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGa                                                                                                                                                               <  2:464018/275‑1 (MQ=255)
                                             actgggtgtgtgatTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTTGCGTCAGCTGCGGCTATGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCGGGATGTGCCAAACGTGAAACTGTACGACTGCGCATGGA‑‑GTGAAGGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGAc                                                                                                                                                              <  2:271218/263‑1 (MQ=255)
                                                 gATGCGATATTTTTTGGTCGCGGCGTCAGCTACGAAAAACCAATGATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTtgt                                                                                                                                                          <  2:178477/275‑1 (MQ=255)
                                                                                       aCCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACa                                                                                                                      >  2:164593/1‑273 (MQ=255)
                                                                                           aTTAACGTCAGTTGCGGTTCTGGGTTAACGGCAGGCGTGGTTGTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGAAGGGG‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGAAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAc                                                                                                                  <  4:164591/273‑1 (MQ=255)
                                                                                                    tataagggacagggTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTaa                                                                                                        <  2:12887/262‑1 (MQ=255)
                                                                                                                 gtgtAACGGCAGCCGTGGTTTTGTAAGGACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTc                                                                                           <  2:249851/274‑1 (MQ=255)
                                                                                                                     aaCGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGAGTGTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAg                                                                                        <  2:222593/275‑1 (MQ=255)
                                                                                                                          cAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAAttc                                                                                                                                                                                                               <  8:407871/149‑2 (MQ=255)
                                                                                                                           aGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCCTGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTTTGTTACCTGAGCACAGCACGTGTTCTCCAGCAAAAAAATCCCCACCATTCACCCTCtttttctgtgtcctttta                                                                                 >  4:12887/1‑261 (MQ=255)
                                                                                                                           aGCCGGTGTTTTGTTAGCAGGCGCGACTCTGGATGGGCCAAACGTAAAACTGTACGACGGGGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAagag                                                                                <  1:37908/275‑1 (MQ=255)
                                                                                                                               gTGGGTTTGTTAGCACTCGCGACGCTGGAGGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCa                                                                                                                           <  2:172111/228‑1 (MQ=255)
                                                                                                                                tGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGt                                                                                                                                                                                                         <  5:8000/149‑1 (MQ=255)
                                                                                                                                   ttttGTTAGCACTCTCGACGCTGGATGTGCCAAACGTGAAACTTTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAatt                                                                         <  2:241552/274‑1 (MQ=255)
                                                                                                                                   ttttGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGtttt                                                                                                                                                                                                      <  8:260119/149‑1 (MQ=255)
                                                                                                                                       gTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTAt                                                                                                                                                                                                  <  6:65353/149‑1 (MQ=255)
                                                                                                                                                cgcgACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAgg                                                                                                                                                                   <  2:114372/171‑1 (MQ=255)
                                                                                                                                                           aTGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGtt                                                                                                                                                                 >  2:171033/1‑162 (MQ=255)
                                                                                                                                                               gCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACt                                                                                                                                                             <  4:171031/162‑1 (MQ=255)
                                                                                                                                                               gCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCg                                            <  2:455304/275‑1 (MQ=255)
                                                                                                                                                                 cAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCca                                                                                                                                                                        >  7:262928/1‑149 (MQ=255)
                                                                                                                                                                  aaaCGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGttt                                                                                                                                                                                                       <  6:439033/117‑1 (MQ=255)
                                                                                                                                                                  aaaCGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCAAATTCGGCGGttt                                                                                                                                                                                                       >  5:439033/1‑117 (MQ=255)
                                                                                                                                                                         aaaCTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGGGGCCTTTCCAGCAGGTTGACTTGTGTTACATTAGCAACGGAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGTAGTTAATTATTGtt                                                                   >  4:446738/1‑242 (MQ=255)
                                                                                                                                                                          aaCTGTACGACGGGGCATGGG‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTccg                                  <  2:368311/274‑1 (MQ=255)
                                                                                                                                                                             tGTACGACGGGGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGtttatt                                                               <  2:446743/242‑1 (MQ=255)
                                                                                                                                                                               tACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCtttt                                                                                               >  2:24670/1‑208 (MQ=255)
                                                                                                                                                                                  gacgGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCt                                                                                                                                                                             >  7:190299/1‑127 (MQ=255)
                                                                                                                                                                                  gacgGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCt                                                                                                                                                                             <  8:190299/127‑1 (MQ=255)
                                                                                                                                                                                       cgcATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGc                                                                                       >  2:513142/1‑208 (MQ=255)
                                                                                                                                                                                       cgcATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAatta                                                                        >  2:389212/1‑223 (MQ=255)
                                                                                                                                                                                         cATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGc                                                                                       >  2:397890/1‑206 (MQ=255)
                                                                                                                                                                                         cATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTAACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGtat                                                       >  2:287204/1‑238 (MQ=255)
                                                                                                                                                                                            aggaatgtaagggggggggCGGACAAATTTACCGGTGGGTCTAGTGAAATAAGTATTTTACAGGCAAAAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGGCTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTaataaa                 <  2:483795/264‑1 (MQ=38)
                                                                                                                                                                                             gA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCaaaa                                                                                   <  4:397885/206‑1 (MQ=255)
                                                                                                                                                                                             gA‑‑GNNAATGGGGCGAGGGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATNTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCGGGTTCGCGGCCTTTCCAGCAGGTTGACGTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAGCAGTAACTCTCTTTTTGTCAAGCaaaa                                                                                   <  4:513137/206‑1 (MQ=255)
                                                                                                                                                                                                         ggCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGc                                                                                                                                                                                               >  7:413250/1‑88 (MQ=255)
                                                                                                                                                                                                         ggCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGc                                                                                                                                                                                               <  8:413250/88‑1 (MQ=255)
                                                                                                                                                                                                                        ttACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATa                                                                                                                   <  8:262928/149‑1 (MQ=255)
                                                                                                                                                                                                                            cGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCtt                                      >  2:380368/1‑222 (MQ=255)
                                                                                                                                                                                                                              gTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcaca                                                                                                                              >  5:46596/1‑132 (MQ=255)
                                                                                                                                                                                                                              gTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTcaca                                                                                                                              <  6:46596/132‑1 (MQ=255)
                                                                                                                                                                                                                                tGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACa                                                                                                                      >  4:391852/1‑138 (MQ=255)
                                                                                                                                                                                                                                    ccAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATAc                                                                                                                  <  2:391857/138‑1 (MQ=255)
                                                                                                                                                                                                                                               aGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGt                                                                                                                                      >  7:6785/1‑107 (MQ=255)
                                                                                                                                                                                                                                               aGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGt                                                                                                                                      <  8:6785/107‑1 (MQ=255)
                                                                                                                                                                                                                                                           gCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAagag                                                                                >  6:71056/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           gCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCt      >  4:384026/1‑223 (MQ=255)
                                                                                                                                                                                                                                                            cAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAagaga                                                                               <  5:71056/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               tAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCttatt  <  2:384031/223‑1 (MQ=255)
                                                                                                                                                                                                                                                               tAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAc                                                                                                              <  6:124104/115‑1 (MQ=255)
                                                                                                                                                                                                                                                               tAAAAAACCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCAc                                                                                                              >  5:124104/1‑115 (MQ=255)
                                                                                                                                                                                                                                                                   aaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCttttt                                                                                              <  6:182440/127‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCttttt                                                                                              >  5:182440/1‑127 (MQ=255)
                                                                                                                                                                                                                                                                   aaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAatta                                                                        >  6:96547/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAatta                                                                        >  6:446758/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aaaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAatta                                                                        >  8:376709/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCttttt                                                                                              <  6:464247/126‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaCCGCCAAATTCGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCttttt                                                                                              >  5:464247/1‑126 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                  
GTGGACGGAACTGGTGCGCGAAGGCGAACTAAAAACGACCGATGAACTGGATGCGATATTTTTTGGTCGCGGCGTCAGCTACGACAAACCAATTATCGTCAGCTGCGGCTCTGGTGTAACGGCAGCCGTGGTTTTGTTAGCACTCGCGACGCTGGATGTGCCAAACGTGAAACTGTACGACGGCGCATGGA‑‑GTGAATGGGGCGCGCGGGCAGATTTACCGGTTGAGCCAGTGAAATAAGTATTTTACAGGCAATAAAAAACCGCCGAATTTGGCGGTTTTTTATTGCTAGTCTGGTTCGCGGCCTTTCCAGCAGGTTGACTTGTGTTACATGAGCAACGCAGGTGCTTCACAGCAAAACAATACTCACCAGTAACTCTCTTTTTGTCAAGCAAAAGAGAGTAATTATTGTTTATTTAGCGTATTATCGACACCGGCCCTTTCCGCCGTGTTCGGTAATAAAATAACCTGGCTTATT  >  NC_000913/2653102‑2653577

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