Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 191,739 T→C 6.3% intergenic (+31/‑116) tsf → / → pyrH translation elongation factor EF‑Ts/uridylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913191,7390TC 6.3% 612.4 / 22.2 190intergenic (+31/‑116)tsf/pyrHtranslation elongation factor EF‑Ts/uridylate kinase
Reads supporting (aligned to +/- strand):  ref base T (96/82);  new base C (6/6);  total (102/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-02

Only 100 of 197 total aligned reads displayed.

TGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAAC  >  NC_000913/191593‑191885
                                                                                                                                                   |                                                                                                                                                  
 gAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGcttctt                                                                                                                                                 <  1:691686/149‑1 (MQ=255)
 gAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGcttctt                                                                                                                                                 >  3:345414/1‑149 (MQ=255)
  aaaGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCttt                                                                                                                                                <  4:675750/149‑1 (MQ=255)
  aaaGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCttt                                                                                                                                                >  1:594363/1‑149 (MQ=255)
  aaaGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCttt                                                                                                                                                <  1:387256/149‑1 (MQ=255)
  aaaGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCttt                                                                                                                                                >  3:167992/1‑149 (MQ=255)
    agagCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCttttt                                                                                                                                              <  2:366805/149‑1 (MQ=255)
    agagCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCttttt                                                                                                                                              <  2:1446049/149‑1 (MQ=255)
    agagCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCttttt                                                                                                                                              <  2:1100368/149‑1 (MQ=255)
    agagCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGAAGCCG‑CCCTCAGGCGGCTCCttttt                                                                                                                                              <  2:1358700/149‑1 (MQ=255)
        cATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGc                                                                                                                                          <  3:193565/149‑1 (MQ=255)
        cATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGAAGCCG‑CCCTCAGGCGGCTCCTTTTtgata                                                                                                                                          <  3:1332906/149‑4 (MQ=38)
         aTAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGcc                                                                                                                                         >  3:246191/1‑149 (MQ=255)
         aTAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGcc                                                                                                                                         >  2:1044596/1‑149 (MQ=255)
             cGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGcttct                                                                                                                                                  <  3:874578/136‑1 (MQ=255)
                 gAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTTTTTTTGTGCCCATCTTGt                                                                                                                                 <  1:1053642/149‑1 (MQ=255)
                  aaGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTa                                                                                                                                >  1:760658/1‑149 (MQ=255)
                      gACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAAtt                                                                                                                            >  1:1028214/1‑149 (MQ=255)
                       aCTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTc                                                                                                                           <  3:291749/149‑1 (MQ=255)
                       aCTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTc                                                                                                                           >  2:100671/1‑149 (MQ=255)
                         tGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCTATCTTGTAAATTCAg                                                                                                                         <  2:1327540/149‑1 (MQ=255)
                          ggCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGc                                                                                                                        >  2:698938/1‑149 (MQ=255)
                          ggCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGc                                                                                                                        >  2:983545/1‑149 (MQ=255)
                          ggCTTCATCCGCGTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGc                                                                                                                        <  2:1236677/149‑1 (MQ=255)
                              tCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAAc                                                                                                                    >  4:846145/1‑149 (MQ=255)
                              tCATCCGCTTCGAAGGGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCAAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAAc                                                                                                                    >  4:171339/1‑149 (MQ=255)
                                    gCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTtgt                                                                                                              <  4:688311/149‑1 (MQ=255)
                                    gCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTtgt                                                                                                              <  3:967625/149‑1 (MQ=255)
                                    gCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTtgt                                                                                                              <  2:707948/149‑1 (MQ=255)
                                     cTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTtgtg                                                                                                             >  1:1084875/1‑149 (MQ=255)
                                         gAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGc                                                                                                         >  4:536172/1‑149 (MQ=255)
                                               ggTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCt                                                                                                   >  4:390411/1‑149 (MQ=255)
                                                gTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTg                                                                                                  >  4:337642/1‑149 (MQ=255)
                                                 tGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGa                                                                                                 >  4:1223260/1‑149 (MQ=255)
                                                    aGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGaaaa                                                                                              <  1:36706/149‑1 (MQ=255)
                                                    aGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGaaaa                                                                                              >  1:281263/1‑149 (MQ=255)
                                                      gCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGc                                                                                            <  4:1350054/149‑1 (MQ=255)
                                                         tCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGAc                                                                                         <  2:726535/149‑1 (MQ=255)
                                                         tCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGAc                                                                                         >  3:368815/1‑149 (MQ=255)
                                                           gagaAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGt                                                                                       >  3:337171/1‑149 (MQ=255)
                                                           gagaAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGt                                                                                       >  4:498099/1‑149 (MQ=255)
                                                            agaAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTa                                                                                      >  3:271756/1‑149 (MQ=255)
                                                                 gTTGAGACTGACTTTGCAGCAGCAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATg                                                                                 >  4:1232387/1‑149 (MQ=255)
                                                                 gTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATg                                                                                 >  3:842931/1‑149 (MQ=255)
                                                                  ttGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGt                                                                                >  4:190865/1‑149 (MQ=255)
                                                                  ttGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGt                                                                                <  2:762647/149‑1 (MQ=255)
                                                                   tGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTc                                                                               <  4:1286404/149‑1 (MQ=255)
                                                                     agaCTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGc                                                                             <  1:1369492/149‑1 (MQ=255)
                                                                       actgactTTGCAGCAGAAGTTGCTGCGATGTCCGAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTa                                                                           <  1:872118/149‑1 (MQ=255)
                                                                        ctgactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAg                                                                          >  3:920930/1‑149 (MQ=255)
                                                                         tgactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACaa                                                                                   >  3:589695/1‑139 (MQ=255)
                                                                         tgactTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACaa                                                                                   <  4:589695/139‑1 (MQ=255)
                                                                            ctTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGGGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTAtt                                                                      <  4:604123/149‑1 (MQ=255)
                                                                              ttGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTaa                                                                    <  4:225235/149‑1 (MQ=255)
                                                                                         gTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAAt                                                         >  3:992285/1‑149 (MQ=255)
                                                                                             ctgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGtt                                                     <  3:158635/149‑1 (MQ=255)
                                                                                             ctgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGtt                                                     <  2:1084875/149‑1 (MQ=255)
                                                                                             ctgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGtt                                                     <  2:803537/149‑1 (MQ=255)
                                                                                              tgcGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGtt                                                     <  3:507337/148‑1 (MQ=255)
                                                                                                     tCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTTAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGtct                                             >  4:905063/1‑149 (MQ=255)
                                                                                                     tCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGtct                                             >  3:665259/1‑149 (MQ=255)
                                                                                                     tCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGtct                                             <  1:111953/149‑1 (MQ=255)
                                                                                                      ccAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTgc                                                                                                      <  4:92671/91‑1 (MQ=255)
                                                                                                      ccAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGtctc                                            <  2:1536046/149‑1 (MQ=255)
                                                                                                       cAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGtt                                                     >  2:47085/1‑139 (MQ=255)
                                                                                                       cAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCa                                           >  1:837591/1‑149 (MQ=255)
                                                                                                       cAAGCAGTCCTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCa                                           >  3:496485/1‑149 (MQ=255)
                                                                                                        aaGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTTGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAg                                          <  4:368815/149‑1 (MQ=255)
                                                                                                        aaGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGtctc                                            <  2:1207249/147‑1 (MQ=255)
                                                                                                            aGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaa                                      <  3:1415821/149‑1 (MQ=255)
                                                                                                            aGTCTTAATTATCAAAAAGGAGCAG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAACAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGCCTCAGgaaa                                      <  4:1423878/149‑1 (MQ=255)
                                                                                                               cTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaagaa                                   <  4:351622/149‑1 (MQ=255)
                                                                                                                ttAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaagaaa                                  >  1:1198200/1‑149 (MQ=255)
                                                                                                                   aTTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACAt                               >  2:1150407/1‑149 (MQ=255)
                                                                                                                   aTTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACAt                               >  3:601680/1‑149 (MQ=255)
                                                                                                                     tATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATgg                             <  3:1223260/149‑1 (MQ=255)
                                                                                                                      aTCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGc                            >  1:1554081/1‑149 (MQ=255)
                                                                                                                        cAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTa                          >  3:224699/1‑149 (MQ=255)
                                                                                                                        cAAAAAGAAGCCGCCCTCA‑GGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGAc                                                                                         <  1:557934/86‑1 (MQ=21)
                                                                                                                        cAAAAAGAAGCCGCCCTCA‑GGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGAc                                                                                         >  2:557934/1‑86 (MQ=21)
                                                                                                                         aaaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAc                         >  4:1494440/1‑149 (MQ=255)
                                                                                                                         aaaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAc                         >  2:145485/1‑149 (MQ=255)
                                                                                                                         aaaaaGAAGCCGCCCTCA‑GGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGc                                                                                            <  4:279504/82‑1 (MQ=21)
                                                                                                                         aaaaaGAAGCCG‑CCCTCAGGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAc                         >  4:1312496/1‑149 (MQ=255)
                                                                                                                          aaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTgg                                                                                                            >  1:1319969/1‑65 (MQ=255)
                                                                                                                          aaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAAt                                                         <  1:1224348/116‑1 (MQ=255)
                                                                                                                          aaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAcc                        >  4:310681/1‑149 (MQ=255)
                                                                                                                          aaaaGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAAAGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAcc                        >  3:1484947/1‑149 (MQ=255)
                                                                                                                          aaaaGAAGCCG‑CCCTCAGGCGGCTCCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTAcc                        >  1:649529/1‑149 (MQ=255)
                                                                                                                             aGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAAt                     <  3:317671/149‑1 (MQ=255)
                                                                                                                             aGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAAt                     <  1:166314/149‑1 (MQ=255)
                                                                                                                                  ccg‑ccTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCaaa                >  2:1491955/1‑149 (MQ=255)
                                                                                                                                  ccg‑ccTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCaaa                >  1:150794/1‑149 (MQ=255)
                                                                                                                                  ccg‑ccTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCaaa                >  2:1520591/1‑149 (MQ=255)
                                                                                                                                               ggcTTCTTTTTGTGCCCATCTTGTAAATTCAGCTGACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGa                                                   <  4:109164/102‑1 (MQ=255)
                                                                                                                                               ggcTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGa                                                   >  3:109164/1‑102 (MQ=255)
                                                                                                                                                gcTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaa                                      <  2:1294846/114‑1 (MQ=255)
                                                                                                                                                gcTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGgaaa                                      >  1:1294846/1‑114 (MQ=255)
                                                                                                                                                gcTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATaaa   >  1:710538/1‑149 (MQ=255)
                                                                                                                                                 cttcttTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAAc  >  1:1025065/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGAAAGAGCATAACGCTGAAGTGACTGGCTTCATCCGCTTCGAAGTGGGTGAAGGCATCGAGAAAGTTGAGACTGACTTTGCAGCAGAAGTTGCTGCGATGTCCAAGCAGTCTTAATTATCAAAAAGGAGCCG‑CCTGAGGGCGGCTTCTTTTTGTGCCCATCTTGTAAATTCAGCTAACCCTTGTGGGGCTGCGCTGAAAAGCGACGTACAATGTCGCTAGTATTAATTCATTTCAATCGTTGACAGTCTCAGGAAAGAAACATGGCTACCAATGCAAAACCCGTCTATAAAC  >  NC_000913/191593‑191885

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑