Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,551,422 T→A 19.7% intergenic (‑126/‑214) gntP ← / → uxuA fructuronate transporter/mannonate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,551,4220TA19.7% 283.3 / 61.0 131intergenic (‑126/‑214)gntP/uxuAfructuronate transporter/mannonate hydrolase
Reads supporting (aligned to +/- strand):  ref base T (53/52);  new base A (11/15);  total (64/67)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

GTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGC  >  NC_000913/4551323‑4551519
                                                                                                   |                                                                                                 
gttgttATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAt                                                                                                   >  2:1244351/1‑100 (MQ=255)
    ttATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgc                                                                                               <  1:1975185/100‑1 (MQ=255)
    ttATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgc                                                                                               <  2:139201/100‑1 (MQ=255)
    ttATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACACCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgc                                                                                               <  1:2274063/100‑1 (MQ=255)
     tATTTTGGGCGATATGTTATGTAAATTGGTCAGCCATTGATGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgca                                                                                              <  1:1820972/100‑1 (MQ=255)
     tATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgca                                                                                              <  2:2576055/100‑1 (MQ=255)
      aTTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgcag                                                                                             >  1:381603/1‑100 (MQ=255)
      aTTTTGGGCGATATGTTATCTAAACTGGTCAACCCTGCTCGCGATGAATGTCCAATCCTCTGATCAATACCCATCCAGTACCCCTTCCTGCAATATgca                                                                                              <  1:1022604/99‑1 (MQ=255)
       tattGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgcagc                                                                                            >  2:1745422/3‑100 (MQ=255)
        tttGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgcagca                                                                                           <  2:1671655/100‑1 (MQ=255)
         ttGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGcaa                                                                                          >  1:992755/1‑100 (MQ=255)
          tGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATGACCATCGATTACCCTTTGCTGCAATTTGCAGcaac                                                                                         <  2:770749/100‑1 (MQ=255)
          tGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaa                                                                                         <  1:399293/100‑2 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCTTTGTTGCGATGAATGTCACATCCTCTGATCAATAACTATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                        <  1:1552678/100‑3 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                        <  1:1429728/100‑3 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                        <  2:2230491/100‑3 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                        <  2:1910706/100‑3 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                        <  1:2593986/100‑3 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                        <  2:2407095/100‑3 (MQ=255)
            ggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATGGATTACCCTTTGCTGCATATTGCAGcaaaag                                                                                       <  2:1742103/100‑4 (MQ=255)
            ggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                        <  1:1727163/99‑3 (MQ=255)
             gCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCaacaac                                                                                      >  2:37248/1‑100 (MQ=255)
             gCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCaacaac                                                                                      <  2:6330/100‑1 (MQ=255)
             gCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaaggg                                                                                      <  2:2383881/100‑5 (MQ=255)
             gCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaaggg                                                                                      <  2:2277503/100‑5 (MQ=255)
             gCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACAGCATCTGATCAATAACCAGCGATTACCCTTTGCTGCAAATTGCAGcaaaggg                                                                                      <  1:274989/100‑5 (MQ=255)
              cGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAAcc                                                                                     <  1:1005632/100‑1 (MQ=255)
               gATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCa                                                                                    <  1:1193681/100‑1 (MQ=255)
                 tatgttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATg                                                                                  <  2:187087/100‑1 (MQ=255)
                 natgttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATg                                                                                  >  2:789219/2‑100 (MQ=255)
                   tgttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATga                                                                                 >  1:2114165/1‑99 (MQ=255)
                    gttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATgaga                                                                               >  2:270440/1‑100 (MQ=255)
                    gttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATgaga                                                                               >  1:1500939/1‑100 (MQ=255)
                     ttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGagag                                                                              <  2:1309585/100‑1 (MQ=255)
                      tatgtAAATTTGTCAACCATTGTATCGATGAATGTAACATCCTCTGATCAAGAACCATCGATTACCTTTTGCTGCAATTTGCGGCAACAACCATGagag                                                                              <  1:52546/99‑1 (MQ=255)
                        tgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGa                                                                           >  1:652027/1‑100 (MQ=255)
                          tAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGaaa                                                                         >  2:1324542/1‑100 (MQ=255)
                           aaaTTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAAt                                                                        >  2:477553/1‑100 (MQ=255)
                            aaTTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAAtt                                                                       <  1:243380/100‑1 (MQ=255)
                               tGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCtt                                                                    >  2:1482011/1‑100 (MQ=255)
                                ggTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtg                                                                   >  1:135650/1‑100 (MQ=255)
                                ggTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtg                                                                   <  1:509309/100‑1 (MQ=255)
                                 gTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtgt                                                                  >  2:1076781/1‑100 (MQ=255)
                                 gTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtgt                                                                  <  1:1687966/100‑1 (MQ=255)
                                 gTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtgt                                                                  >  1:1501297/1‑100 (MQ=255)
                                   cAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGa                                                                <  2:1501306/100‑1 (MQ=255)
                                    aaCCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGAt                                                               >  1:1983743/1‑100 (MQ=255)
                                          tgttgCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGtt                                                         <  1:776789/100‑1 (MQ=255)
                                           gttgCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTa                                                        >  1:1360339/1‑100 (MQ=255)
                                           gttgCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTa                                                        >  2:1296630/1‑100 (MQ=255)
                                             tgCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCCGAAACAACCATAAGAGTGAAATTCTTGTGATGTGGTTAAc                                                      <  1:763601/100‑1 (MQ=255)
                                              gCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAAcc                                                     <  2:331522/100‑1 (MQ=255)
                                               cGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCa                                                    >  2:1833358/1‑100 (MQ=255)
                                                 aTGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAAt                                                  <  1:1838649/100‑1 (MQ=255)
                                                    aaTGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTc                                               >  2:1588692/1‑100 (MQ=255)
                                                     aTGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCa                                              >  1:2500378/1‑100 (MQ=255)
                                                            aTCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCg                                       <  1:1995480/100‑1 (MQ=255)
                                                              ccTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCggg                                     >  2:78850/1‑100 (MQ=255)
                                                              ccTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCggg                                     >  1:159987/1‑100 (MQ=255)
                                                                tctGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGtt                                   <  1:1891146/100‑1 (MQ=255)
                                                                tctGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGtt                                   >  1:1966873/1‑100 (MQ=255)
                                                                tctGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGtt                                   <  2:1747081/100‑1 (MQ=255)
                                                                 ctGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTg                                  <  2:2553465/100‑1 (MQ=255)
                                                                   gATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAGTTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGAc                                <  1:1866146/100‑1 (MQ=255)
                                                                   gATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGa                                 >  1:278524/1‑99 (MQ=255)
                                                                   gATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGAc                                <  1:1866424/100‑1 (MQ=255)
                                                                   gATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGAc                                >  2:1978889/1‑100 (MQ=255)
                                                                   gATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGAc                                >  2:1909047/1‑100 (MQ=255)
                                                                    aTCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACa                               <  1:1077486/100‑1 (MQ=255)
                                                                     tCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACAt                              >  2:39353/1‑100 (MQ=255)
                                                                     tCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACAt                              >  1:1107571/1‑100 (MQ=255)
                                                                      cAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATg                             >  1:1309785/1‑100 (MQ=255)
                                                                       aaTAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGt                            >  1:738133/1‑100 (MQ=255)
                                                                       aaTAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGt                            >  2:1923663/1‑100 (MQ=255)
                                                                        aTAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTc                           <  1:2446387/100‑1 (MQ=255)
                                                                          aaCCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCtt                         >  2:181138/1‑100 (MQ=255)
                                                                           aCCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCt                          >  1:1333235/1‑98 (MQ=255)
                                                                            ccATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGCTAACCAATTGCAGAATTCGGGTTGACATGTCTTAc                       >  1:2271190/1‑100 (MQ=255)
                                                                             cATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTAcc                      <  2:545622/100‑1 (MQ=255)
                                                                               tCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCaa                    >  2:115970/1‑100 (MQ=255)
                                                                                 gATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCaaaa                  >  1:370673/1‑100 (MQ=255)
                                                                                   ttACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAgg                <  1:1842780/100‑1 (MQ=255)
                                                                                   ttACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAgg                <  2:334532/100‑1 (MQ=255)
                                                                                   ttACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACAAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAgg                <  1:939125/100‑1 (MQ=255)
                                                                                     aCCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTa              >  2:586456/1‑100 (MQ=255)
                                                                                      gttgttGCTGCAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAg             >  2:2513655/5‑100 (MQ=255)
                                                                                      ggtgttGCTGCACATTGGAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGCTAACCAAGCTCAGAACTCGGGTTGACATGTCTTACCAAAAg                 >  1:2087140/5‑96 (MQ=255)
                                                                                      cccTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAg             >  2:1208572/1‑100 (MQ=255)
                                                                                        tgttGCTGCAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGaa           >  2:163121/3‑100 (MQ=255)
                                                                                        tgttGCTGCAAATGGCAGCAACAACCATGAGAGGGAAATTCTTTTGATGTGGTTAACAAATTTCAGAATTCTGGTTGACATGTCTTATCAAAAGGTAGaa           >  1:2196851/3‑100 (MQ=255)
                                                                                        cTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGaa           >  1:2095873/1‑100 (MQ=255)
                                                                                         gttGCTGCAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAAc          >  1:113005/2‑100 (MQ=255)
                                                                                         cttGCTGCAAATTGCAGCAACAACCATGAGATTAAAATTCCTGTGATGTGGTTAACCAATTTCAGCATTCGGTTTGACATGTCTTACCAAAAGGTAGAAc          >  1:915137/2‑100 (MQ=255)
                                                                                          ttGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACt         <  2:1937246/100‑1 (MQ=255)
                                                                                           tgctgcAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACtt        >  1:144884/1‑100 (MQ=255)
                                                                                            gctgcAATTTGCAGCAACAACCCTGAGAGTGAAATTCTTGTGATGTGGTGAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACtt        >  2:688362/1‑99 (MQ=255)
                                                                                             ctgcAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTtat      >  2:666816/1‑100 (MQ=255)
                                                                                               gcAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATAc    >  1:1183044/1‑100 (MQ=255)
                                                                                                 aaTTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGc  <  2:717072/100‑1 (MQ=255)
                                                                                                 aaTTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGc  <  2:717181/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGC  >  NC_000913/4551323‑4551519

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