Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 87,881 A→G 6.8% intergenic (+33/‑147) 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,8810AG 6.8% 301.9 / 9.8 103intergenic (+33/‑147)ilvH/craacetolactate synthase 3, small subunit, valine‑sensitive/transcriptional repressor‑activator for carbon metabolism
Reads supporting (aligned to +/- strand):  ref base A (42/54);  new base G (4/3);  total (46/57)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

Only 100 of 107 total aligned reads displayed.

GCACCAGCGGTAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATTGT  >  NC_000913/87733‑88029
                                                                                                                                                    |                                                                                                                                                    
gCACCAGCGGTAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACa                                                                                                                                                      >  1:858372/1‑149 (MQ=255)
      gCGGTAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCAc                                                                                                                                                >  3:336368/1‑149 (MQ=255)
        ggTAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGt                                                                                                                                              >  1:594848/1‑149 (MQ=255)
          tAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTg                                                                                                                                            <  4:238534/149‑1 (MQ=255)
          tAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTg                                                                                                                                            <  1:639313/149‑1 (MQ=255)
           aaGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTgg                                                                                                                                           <  3:35859/149‑1 (MQ=255)
            aGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTggg                                                                                                                                          <  2:337821/149‑1 (MQ=255)
                   tGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACGTGACATGTTGGGCtgtaaa                                                                                                                                   <  1:296923/149‑6 (MQ=255)
                   tGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCtttttt                                                                                                                                   <  2:288838/149‑1 (MQ=255)
                   tGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCtttttt                                                                                                                                   >  4:465853/1‑149 (MQ=255)
                    gCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACGTGACATGTTGGGc                                                                                                                                         <  2:738093/142‑1 (MQ=255)
                    gCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACGTGACATGTTGGGc                                                                                                                                         >  1:738093/1‑142 (MQ=255)
                    gCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCttttttt                                                                                                                                  >  3:304326/1‑149 (MQ=255)
                    gCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCttttttt                                                                                                                                  <  2:704194/149‑1 (MQ=255)
                    gCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCttttttt                                                                                                                                  <  1:422069/149‑1 (MQ=255)
                      aTTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGc                                                                                                                                <  1:827819/149‑1 (MQ=255)
                      aTTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGc                                                                                                                                >  4:664711/1‑149 (MQ=255)
                      aTTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGc                                                                                                                                <  2:197224/149‑1 (MQ=255)
                               aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCg                                                                                                                               >  4:336775/1‑141 (MQ=255)
                               aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCg                                                                                                                               <  3:336775/141‑1 (MQ=255)
                               aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGt                                                                                                                       <  2:1023368/149‑1 (MQ=255)
                               aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGt                                                                                                                       <  4:336368/149‑1 (MQ=255)
                               aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGt                                                                                                                       <  3:145426/149‑1 (MQ=255)
                               aTCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGt                                                                                                                       <  3:354067/149‑1 (MQ=255)
                                 cGATTAGCGATGTGGCGAATATTGTGGAGATTACTCGCTCTGGTGTGGTCGGACTTTCGCGCGACGATAAAATAATGCCTTCAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCATGTTGCGCTTTTAGTACAAAACCAGTgg                                                                                                                     <  2:477431/149‑1 (MQ=255)
                                  gATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTggg                                                                                                                    <  2:96539/149‑1 (MQ=255)
                                    ttCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGCCTTTCGCGCGGCGATGAAATAATTCGTTGAGAATGATCTCAATGCGCTATTTACAGCCCAATATTTAACGTTGGGCTTTTTTTGCGGAATCAGTgggta                                                                                                                  >  2:1101040/1‑147 (MQ=255)
                                       gcgATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCt                                                                                                               <  1:690302/149‑1 (MQ=255)
                                          aTGTGGCGGAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGTACTTTCGCGCGGCGATAAAATGATGCGTTGGGAATGGTCTCAATGCGCAATTTCCAGCCCAACATGTCTCGTTGGGCTTTTTTTGCGACATCAGTGGGAACCTGGa                                                                                                            >  2:1088729/1‑149 (MQ=255)
                                          aTGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGa                                                                                                            <  3:372863/149‑1 (MQ=255)
                                           tgtgGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGaa                                                                                                           <  2:82543/149‑1 (MQ=255)
                                               gCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTg                                                                                                              >  4:514346/1‑142 (MQ=255)
                                               gCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTg                                                                                                              <  3:514346/142‑1 (MQ=255)
                                               gCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATaaa                                                                                                       <  4:666564/149‑1 (MQ=255)
                                                cGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCTCGTCGATAAAATAATGCGCTGAGAAGTATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAGCCTGGCATaaaa                                                                                                      <  2:1000651/149‑1 (MQ=255)
                                                cGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATaaaa                                                                                                      <  2:39769/149‑1 (MQ=255)
                                                cGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCCCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATaaaa                                                                                                      <  4:623781/149‑1 (MQ=255)
                                                 gAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAg                                                                                                     >  4:24669/1‑149 (MQ=255)
                                                 gAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAg                                                                                                     >  2:286426/1‑149 (MQ=255)
                                                 gAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAg                                                                                                     >  4:372575/1‑149 (MQ=255)
                                                 gAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAg                                                                                                     >  4:457932/1‑149 (MQ=255)
                                                  aaaaTTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGc                                                                                                    <  3:654494/149‑1 (MQ=255)
                                                      ttGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGtt                                                                                                <  4:525790/149‑1 (MQ=255)
                                                        gtgGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGc                                                                                              <  2:594848/149‑1 (MQ=255)
                                                            aGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCa                                                                                          >  3:518291/1‑149 (MQ=255)
                                                              gTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGt                                                                                        <  3:850909/149‑1 (MQ=255)
                                                              gTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGt                                                                                        >  3:370985/1‑149 (MQ=255)
                                                                 gctcgctcTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTaa                                                                                     <  4:892613/149‑1 (MQ=255)
                                                                       tctGGTGTGGTCGGACTTTCGTGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCt                                                                               <  4:224721/149‑1 (MQ=255)
                                                                       tctGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAGTAAAAGCAGTTGCCGCAGTTAATTTTcc                                                                               >  4:930615/1‑148 (MQ=255)
                                                                       tctGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCt                                                                               >  3:479969/1‑149 (MQ=255)
                                                                                 tCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGCACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGAt                                                                     >  4:590800/1‑149 (MQ=255)
                                                                                 tCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGAt                                                                     >  2:257973/1‑149 (MQ=255)
                                                                                 tCGGACTGTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGAt                                                                     >  4:347488/1‑149 (MQ=255)
                                                                                    gACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGtt                                                                  >  3:1017881/1‑149 (MQ=255)
                                                                                    gACTTTCGCGCGACAAAGAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCAGGTTGGGCTGTTTTTGCGAAATCAGTGGGAACCTGGAGTAACAGCAGTTGCCGCAGTTAATTTTCTGCGCTGAGATGtt                                                                  >  2:790862/1‑149 (MQ=255)
                                                                                      cTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTaa                                                                <  3:465853/149‑1 (MQ=255)
                                                                                         tCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGa                                                             <  1:961374/149‑1 (MQ=255)
                                                                                           gcgcgGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAAt                                                           >  3:744039/1‑149 (MQ=255)
                                                                                                      aaaTAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACACGACACGTTGGGCTTTTTTTGCGAAATCTGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTGGATGTTAATGAATTTAACCCATAc                                                <  3:913764/149‑1 (MQ=255)
                                                                                                       aataatGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATAcc                                               >  2:1097536/1‑149 (MQ=255)
                                                                                                               gTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAAt                                       <  1:20819/149‑1 (MQ=255)
                                                                                                                ttGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTACCCCATACCAGTACAATg                                      <  2:314179/149‑1 (MQ=255)
                                                                                                                ttGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATg                                      >  2:204790/1‑149 (MQ=255)
                                                                                                                ttGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTGCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATg                                      <  3:590090/149‑1 (MQ=255)
                                                                                                                  gagaATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGc                                    <  3:664711/149‑1 (MQ=255)
                                                                                                                    gaATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATAc                                                <  4:470958/135‑1 (MQ=255)
                                                                                                                    gaATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATAc                                                >  3:470958/1‑135 (MQ=255)
                                                                                                                    gaATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTa                                  >  3:1054859/1‑149 (MQ=255)
                                                                                                                     aaTGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCttt                                                                                                                                      <  2:458127/48‑1 (MQ=255)
                                                                                                                     aaTGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCttt                                                                                                                                      >  1:458127/1‑48 (MQ=255)
                                                                                                                            tcAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGttttt                          >  1:158507/1‑149 (MQ=255)
                                                                                                                            tcAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGttttt                          <  4:82833/149‑1 (MQ=255)
                                                                                                                               aTGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACa                       <  1:1097213/149‑1 (MQ=255)
                                                                                                                               aTGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACa                       <  4:1017881/149‑1 (MQ=255)
                                                                                                                                   gcAATTTACAGCCCAACATGTCACGGTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTAGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACAtttt                   >  4:908245/1‑149 (MQ=255)
                                                                                                                                    cAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTa                  <  1:1097536/149‑1 (MQ=255)
                                                                                                                                     aaTTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATTTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTAc                 <  4:916497/149‑1 (MQ=255)
                                                                                                                                     aaTTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAAt                                       >  3:260771/1‑127 (MQ=255)
                                                                                                                                     aaTTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAAt                                       <  4:260771/127‑1 (MQ=255)
                                                                                                                                     aaTTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTAc                 >  2:86424/1‑149 (MQ=255)
                                                                                                                                     aaTTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTAc                 >  1:300794/1‑149 (MQ=255)
                                                                                                                                     aaTTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTAc                 >  4:64183/1‑149 (MQ=255)
                                                                                                                                     aaTTTACAACCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAAACTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTAc                 >  4:820165/1‑149 (MQ=255)
                                                                                                                                      aaaaaaaGCCCAACGTGACATGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACg                >  1:651307/7‑149 (MQ=255)
                                                                                                                                      aaaaaaaGCCCAACGTGACATGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACg                >  2:359947/7‑149 (MQ=255)
                                                                                                                                          aCAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAg            <  1:204790/149‑1 (MQ=255)
                                                                                                                                            aGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAggg          <  3:152885/149‑1 (MQ=255)
                                                                                                                                             gCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAgggg         <  3:64183/149‑1 (MQ=255)
                                                                                                                                             gCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGAACAATGGCTATGGTTTTTACATTTTACGCAAgggg         >  2:511480/1‑149 (MQ=255)
                                                                                                                                              cccAACGTGACATGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAAt                                                           <  1:103119/98‑1 (MQ=255)
                                                                                                                                              cccAACGTGACATGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAAt                                                           >  2:103119/1‑98 (MQ=255)
                                                                                                                                              cccAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGc        >  4:86822/1‑149 (MQ=255)
                                                                                                                                               ccAACATGTCACGTTGGGCATTTTTTGCGTAATCAGTGGGTACCTAGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGTCTTTTACATTTTACGCAAGGGGCa       >  4:176549/1‑149 (MQ=255)
                                                                                                                                                cAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCaa      >  2:799511/1‑149 (MQ=255)
                                                                                                                                                   cATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGttt                            <  3:292453/124‑1 (MQ=255)
                                                                                                                                                   cATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGttt                            >  4:292453/1‑124 (MQ=255)
                                                                                                                                                    aTGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATtgt  >  3:39471/1‑149 (MQ=255)
                                                                                                                                                    aTGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATtgt  <  3:820165/149‑1 (MQ=255)
                                                                                                                                                    aTGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATtgt  <  4:372732/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCACCAGCGGTAAGCTTGATGCATTTTTAGCATCGATTCGCGATGTGGCGAAAATTGTGGAGGTTGCTCGCTCTGGTGTGGTCGGACTTTCGCGCGGCGATAAAATAATGCGTTGAGAATGATCTCAATGCGCAATTTACAGCCCAACATGTCACGTTGGGCTTTTTTTGCGAAATCAGTGGGAACCTGGAATAAAAGCAGTTGCCGCAGTTAATTTTCTGCGCTTAGATGTTAATGAATTTAACCCATACCAGTACAATGGCTATGGTTTTTACATTTTACGCAAGGGGCAATTGT  >  NC_000913/87733‑88029

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