Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 87,884 T→A 5.0% intergenic (+36/‑144) ilvH → / → cra acetolactate synthase 3, small subunit, valine‑sensitive/transcriptional repressor‑activator for carbon metabolism

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091387,8840TA 5.0% 345.4 / 6.8 99intergenic (+36/‑144)ilvH/craacetolactate synthase 3, small subunit, valine‑sensitive/transcriptional repressor‑activator for carbon metabolism
Reads supporting (aligned to +/- strand):  ref base T (44/50);  new base A (3/2);  total (47/52)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

CAGCGGTAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATTGT  >  NC_000913/87737‑88029
                                                                                                                                                   |                                                                                                                                                 
cAGCGGTAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTc                                                                                                                                                  <  1:543322/149‑1 (MQ=255)
    ggTAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGTGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGt                                                                                                                                              >  4:400714/1‑149 (MQ=255)
     gTAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGtt                                                                                                                                             >  1:872978/1‑149 (MQ=255)
      tAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTg                                                                                                                                            <  2:805983/149‑1 (MQ=255)
         gCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGc                                                                                                                                         >  1:257662/1‑149 (MQ=255)
           ttGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCtt                                                                                                                                       >  1:347182/1‑149 (MQ=255)
            tGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCttt                                                                                                                                      >  4:701137/1‑149 (MQ=255)
             gATGCAGTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCtttt                                                                                                                                     >  2:621873/1‑149 (MQ=255)
                gCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCttttttt                                                                                                                                  >  3:859194/1‑149 (MQ=255)
                gCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCttttttt                                                                                                                                  <  2:570349/149‑1 (MQ=255)
                gCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCttttttt                                                                                                                                  >  3:169851/1‑149 (MQ=255)
                  aTTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGc                                                                                                                                <  2:650280/149‑1 (MQ=255)
                  aTTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGc                                                                                                                                <  2:590041/149‑1 (MQ=255)
                       tAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAAt                                                                                                                           <  4:659063/149‑1 (MQ=255)
                         gCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCa                                                                                                                         >  2:564151/1‑149 (MQ=255)
                           aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGt                                                                                                                       >  2:27086/1‑149 (MQ=255)
                           aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGt                                                                                                                       <  2:47564/149‑1 (MQ=255)
                           aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGt                                                                                                                       >  1:47564/1‑149 (MQ=255)
                           aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGt                                                                                                                       >  3:807075/1‑149 (MQ=255)
                           aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGt                                                                                                                       >  2:596565/1‑149 (MQ=255)
                           aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGt                                                                                                                       <  2:591036/149‑1 (MQ=255)
                            tCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAGATAATGCGTTGAGAATGATGTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTg                                                                                                                      <  2:718138/149‑1 (MQ=255)
                            tCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTg                                                                                                                      >  4:350427/1‑149 (MQ=255)
                            tCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTg                                                                                                                      >  3:564130/1‑149 (MQ=255)
                                ttCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTgggga                                                                                                                  >  1:686546/1‑147 (MQ=255)
                                    cgATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTg                                                                                                              <  1:729448/149‑1 (MQ=255)
                                    cgATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTg                                                                                                              <  1:564151/149‑1 (MQ=255)
                                     gATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGACGATAAAATAATGCGTCGAGACTCATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTgg                                                                                                             <  2:122039/149‑1 (MQ=255)
                                      aTGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGa                                                                                                            >  4:597610/1‑149 (MQ=255)
                                           gCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATaaa                                                                                                       >  1:743541/1‑149 (MQ=255)
                                           gCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATaaa                                                                                                       >  1:454791/1‑149 (MQ=255)
                                            cGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATaaaa                                                                                                      <  4:862377/149‑1 (MQ=255)
                                            cGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATaaaa                                                                                                      <  2:64764/149‑1 (MQ=255)
                                             gAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAg                                                                                                     >  1:603683/1‑149 (MQ=255)
                                             gAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAg                                                                                                     >  3:371912/1‑149 (MQ=255)
                                             gAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAg                                                                                                     >  2:490299/1‑149 (MQ=255)
                                             gAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAg                                                                                                     <  2:603683/149‑1 (MQ=255)
                                             gAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAg                                                                                                     <  3:92933/149‑1 (MQ=255)
                                                         ggTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAg                                                                                         <  2:405137/149‑1 (MQ=255)
                                                          gTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGt                                                                                        <  3:247058/149‑1 (MQ=255)
                                                               tcgctcTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAAtt                                                                                   >  3:810426/1‑149 (MQ=255)
                                                               tcgctcTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAAtt                                                                                   <  2:601623/149‑1 (MQ=255)
                                                               tcgctcTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAAtt                                                                                   >  1:308147/1‑149 (MQ=255)
                                                                  ctctGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTc                                                                                <  1:596565/149‑1 (MQ=255)
                                                                   tctGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCt                                                                               >  2:36881/1‑149 (MQ=255)
                                                                       gtgtGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTgcgc                                                                           <  1:396938/149‑1 (MQ=255)
                                                                       gtgtGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTgcgc                                                                           >  2:396938/1‑149 (MQ=255)
                                                                         gtgGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCtt                                                                         <  4:810426/149‑1 (MQ=255)
                                                                          tgGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTa                                                                        <  2:144045/149‑1 (MQ=255)
                                                                             tCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGAt                                                                     >  4:727339/1‑149 (MQ=255)
                                                                             tCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGAt                                                                     <  2:488133/149‑1 (MQ=255)
                                                                             tCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGAt                                                                     >  2:879488/1‑149 (MQ=255)
                                                                                  cTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTaa                                                                <  4:13141/149‑1 (MQ=255)
                                                                                  cTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTaa                                                                <  1:36881/149‑1 (MQ=255)
                                                                                   tttCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAAt                                                               <  4:807075/149‑1 (MQ=255)
                                                                                   tttCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAAt                                                               <  4:802512/149‑1 (MQ=255)
                                                                                   tttCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAAt                                                               >  4:787395/1‑149 (MQ=255)
                                                                                   tttCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAAt                                                               >  1:524164/1‑149 (MQ=255)
                                                                                      cgcgcgGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGaa                                                            <  4:371912/149‑1 (MQ=255)
                                                                                         gcggcgATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAAttt                                                         <  4:658677/149‑1 (MQ=255)
                                                                                         gcggcgATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAAttt                                                         >  3:488137/1‑149 (MQ=255)
                                                                                            gcgATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAAc                                                      >  1:772639/1‑149 (MQ=255)
                                                                                            gcgATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAAc                                                      <  1:27086/149‑1 (MQ=255)
                                                                                              gATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAAccc                                                    >  1:346315/1‑149 (MQ=255)
                                                                                                tAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCttttttg                                                                                                                                  >  4:642566/1‑68 (MQ=255)
                                                                                                tAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCttttttg                                                                                                                                  <  3:642566/69‑2 (MQ=255)
                                                                                                tAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCAt                                                  <  1:209625/149‑1 (MQ=255)
                                                                                                   aataatGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAAt                                                               <  4:616851/133‑1 (MQ=255)
                                                                                                   aataatGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAAt                                                               >  3:616851/1‑133 (MQ=255)
                                                                                                      aaggCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAATATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGTAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGt                                            >  4:85702/4‑149 (MQ=255)
                                                                                                          cGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACaa                                        <  2:772639/149‑1 (MQ=255)
                                                                                                           tgcctgaATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAAt                                       <  4:29021/144‑1 (MQ=255)
                                                                                                            ttGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATg                                      >  1:155089/1‑149 (MQ=255)
                                                                                                            ttGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATg                                      <  3:577298/149‑1 (MQ=255)
                                                                                                            ttGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATg                                      <  3:171815/149‑1 (MQ=255)
                                                                                                               agaATGATCTCAATGCGCAATTTACAGCTCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCt                                   >  3:58362/1‑149 (MQ=255)
                                                                                                               agaATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCt                                   >  2:793660/1‑149 (MQ=255)
                                                                                                                 aaagATCTCAATGCGCAATTTACAGCCCAACATGACACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAACGGCTAt                                 >  3:576420/4‑149 (MQ=255)
                                                                                                                 aaTGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTgg                                                                                                                     <  1:398295/65‑1 (MQ=255)
                                                                                                                 aaTGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTgg                                                                                                                     >  2:398295/1‑65 (MQ=255)
                                                                                                                 aaTGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAg                                                                                                     <  4:851943/81‑1 (MQ=255)
                                                                                                                 aaTGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAg                                                                                                     >  3:851943/1‑81 (MQ=255)
                                                                                                                    gATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGt                              >  1:835323/1‑149 (MQ=255)
                                                                                                                      tctcAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGAt                                                                     <  3:449735/108‑1 (MQ=255)
                                                                                                                      tctcAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGAt                                                                     >  4:449735/1‑108 (MQ=255)
                                                                                                                      tctcAATGCGCAATTTACAACCCAATATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGTCGCAGTTAATTGTCTTTGCTTATATGTTAATGAATTTAACACATACCAGTACAATGGCTATGGttt                            <  2:113003/149‑1 (MQ=255)
                                                                                                                       ctcAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGtttt                           <  3:128226/149‑1 (MQ=255)
                                                                                                                              cgcAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACAttt                    <  2:508404/149‑1 (MQ=255)
                                                                                                                                cAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTa                  <  2:347182/149‑1 (MQ=255)
                                                                                                                                     aaaaGCCCGACGTGACATGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACAt                      >  3:799057/4‑140 (MQ=255)
                                                                                                                                     aaaaGCCCAACGTGACATGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACAt                      <  4:799057/137‑1 (MQ=255)
                                                                                                                                      aaaGCCCAACGTGACATGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATg                                >  3:829862/3‑129 (MQ=255)
                                                                                                                                      aaaGCCCAACGTGACATGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATg                                <  4:829862/127‑1 (MQ=255)
                                                                                                                                      aCAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAg            <  2:155089/149‑1 (MQ=255)
                                                                                                                                         gCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAgggg         >  1:110993/1‑149 (MQ=255)
                                                                                                                                         gCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAgggg         <  1:438021/149‑1 (MQ=255)
                                                                                                                                             aaCATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAAt     <  1:621873/149‑1 (MQ=255)
                                                                                                                                              aCATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAAtt    <  1:446264/149‑1 (MQ=255)
                                                                                                                                                aTGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATtgt  <  1:879488/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAGCGGTAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATTGT  >  NC_000913/87737‑88029

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