Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,551,422 T→A 17.2% 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,4220TA17.2% 357.0 / 65.6 152intergenic (‑126/‑214)gntP/uxuAfructuronate transporter/mannonate hydrolase
Reads supporting (aligned to +/- strand):  ref base T (62/63);  new base A (14/12);  total (77/75)
Fisher's exact test for biased strand distribution p-value = 8.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.11e-01

Only 100 of 155 total aligned reads displayed.

GTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCC  >  NC_000913/4551323‑4551520
                                                                                                    |                                                                                                  
gttgttATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAt                                                                                                    <  1:1371022/100‑1 (MQ=255)
gttgttAGTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAt                                                                                                    <  1:1216930/100‑1 (MQ=255)
   gttATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTg                                                                                                 >  2:1913006/1‑100 (MQ=255)
   gttATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTg                                                                                                 >  2:1063893/1‑100 (MQ=255)
      aTTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgcag                                                                                              <  2:2847725/100‑1 (MQ=255)
       ttttGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgcagc                                                                                             <  1:3325606/100‑1 (MQ=255)
       ttttGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgcagc                                                                                             >  1:283253/1‑100 (MQ=255)
       ttttGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTgcagc                                                                                             <  1:2059198/100‑1 (MQ=255)
          taggagaaatgttatgtAAATGGGTCAACCATGGTTGCGATGAATGTTACATCCGCTGATCAATAACCATCGATTACACTTTGCTGCAAATTgaagc                                                                                             <  2:3098161/89‑1 (MQ=255)
          tGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaa                                                                                          <  2:1810318/100‑2 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                         <  1:1149669/100‑3 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                         <  2:2858069/100‑3 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                         <  2:2427940/100‑3 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                         <  1:3494020/100‑3 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                         <  1:2667380/100‑3 (MQ=255)
            ggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGcaacaa                                                                                        >  1:1767608/1‑100 (MQ=255)
            ggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaagg                                                                                        <  1:579079/100‑4 (MQ=255)
             gCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaaggg                                                                                       <  1:1634675/100‑5 (MQ=255)
              cGATATGTTATGTAAATTGGTCAACAATTGTTGCGATGAATGTCACATCCTCTGTTCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGcaacaa                                                                                        >  2:2310575/1‑98 (MQ=255)
               gATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAAcc                                                                                      <  2:1430778/99‑1 (MQ=255)
                atatGTTCTGTAGATCGGTCAACCCTTGTTGCGATGCATGTCCCATCCTCTGAGCAAGAACCCTCGATTCCCCTTTGCTGCAATTTGCATCAACAAcc                                                                                      <  2:3547862/98‑1 (MQ=255)
                atatGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCAt                                                                                    <  2:1409606/100‑1 (MQ=255)
                 tatgttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATg                                                                                   >  1:940908/1‑100 (MQ=255)
                 tatgttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAAATTGCAGcaaagggtaat                                                                                   <  1:3610911/100‑9 (MQ=39)
                    gttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATgaga                                                                                >  2:1712478/1‑100 (MQ=255)
                    gttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATgaga                                                                                >  2:1134871/1‑100 (MQ=255)
                     ttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGagag                                                                               <  1:810817/100‑1 (MQ=255)
                     ttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGagag                                                                               <  1:189646/100‑1 (MQ=255)
                      tatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGt                                                                              >  2:361459/1‑100 (MQ=255)
                        tgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGa                                                                            >  1:3288665/1‑100 (MQ=255)
                         gtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCTTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGaa                                                                           >  2:1533694/1‑100 (MQ=255)
                          tAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGaaa                                                                          >  1:1907014/1‑100 (MQ=255)
                            aaTTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAAtt                                                                        >  1:1969858/1‑100 (MQ=255)
                                ggTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtg                                                                    <  2:2639458/100‑1 (MQ=255)
                                 gTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTTCCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtg                                                                    <  2:2020616/99‑1 (MQ=255)
                                 gTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtgt                                                                   >  2:1166318/1‑100 (MQ=255)
                                 gTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtgt                                                                   >  1:608060/1‑100 (MQ=255)
                                   cAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGa                                                                 >  2:1598416/1‑100 (MQ=255)
                                   cAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGa                                                                 >  1:2692034/1‑100 (MQ=255)
                                   cAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGa                                                                 >  1:2357151/1‑100 (MQ=255)
                                        aTTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTgg                                                            >  1:3566358/1‑100 (MQ=255)
                                         ttgttgCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGt                                                           <  2:2745775/100‑1 (MQ=255)
                                          tgttgCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGtt                                                          >  1:2246149/1‑100 (MQ=255)
                                           gttgCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTa                                                         <  1:1134880/100‑1 (MQ=255)
                                           gttgCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTa                                                         <  2:831101/100‑1 (MQ=255)
                                             tgCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAAc                                                       <  2:2859924/100‑1 (MQ=255)
                                               cGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAAc                                                       >  2:1398125/1‑98 (MQ=255)
                                               cGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCa                                                     >  2:2379871/1‑100 (MQ=255)
                                                gATGAATGCCCCATACTCTGGACAATAAAATCTGAGTACCCTTTGCTACAATTTGCAGCAACAACCAAGAGAGTAAAATTAACGTGATGTGGAAAACCa                                                     >  2:3630767/1‑99 (MQ=255)
                                                    aaTGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTGACCAATTTc                                                >  2:749920/1‑100 (MQ=255)
                                                      tGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAg                                              <  1:540405/100‑1 (MQ=255)
                                                      tGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAg                                              <  1:1718213/100‑1 (MQ=255)
                                                      tGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAg                                              >  1:1516373/1‑100 (MQ=255)
                                                      tGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAg                                              <  2:1016145/100‑1 (MQ=255)
                                                       gTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGa                                             >  1:766366/1‑100 (MQ=255)
                                                         cacaTCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAg                                              <  2:1350976/97‑1 (MQ=255)
                                                         cacaTCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAAt                                           >  2:3504192/1‑100 (MQ=255)
                                                             tCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCgg                                       >  2:2145676/1‑100 (MQ=255)
                                                              ccTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCggg                                      <  1:2219921/100‑1 (MQ=255)
                                                               ctctGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGt                                     <  2:1019710/100‑1 (MQ=255)
                                                               ctctGATCAATAACCAGCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATCCTTGTGATGTGGTTAACCAATTTCAGAATTCggg                                      <  2:1282190/99‑1 (MQ=255)
                                                                nctGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGtt                                    >  1:3623016/2‑100 (MQ=255)
                                                                 ctGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTg                                   <  2:3013047/100‑1 (MQ=255)
                                                                 ctGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTg                                   <  1:1198741/100‑1 (MQ=255)
                                                                   gATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGAc                                 <  2:638572/100‑1 (MQ=255)
                                                                   gATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGAc                                 >  2:635608/1‑100 (MQ=255)
                                                                       aaTAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGt                             >  1:2042500/1‑100 (MQ=255)
                                                                        aTAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTc                            <  1:3176001/100‑1 (MQ=255)
                                                                         tAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATCCTTGTGATGTGGTTAACAAATTTCAGAATACGGGTTGa                                  >  2:3328352/1‑93 (MQ=255)
                                                                          aaCCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCtt                          <  1:2520154/100‑1 (MQ=255)
                                                                          aaCCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCt                           <  1:1576330/99‑1 (MQ=255)
                                                                           aCCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTa                         >  2:1008153/1‑100 (MQ=255)
                                                                            ccATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTAc                        >  2:1928693/1‑100 (MQ=255)
                                                                            ccATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTAc                        >  1:2383508/1‑100 (MQ=255)
                                                                             cATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTAcc                       <  2:50071/100‑1 (MQ=255)
                                                                                   aTTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAg                  <  1:2169452/100‑1 (MQ=255)
                                                                                   aTTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAg                  <  2:2160778/100‑1 (MQ=255)
                                                                                    ttACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAgg                 <  2:3594064/100‑1 (MQ=255)
                                                                                    ttACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAgg                 <  1:301902/100‑1 (MQ=255)
                                                                                    atggttgttGCTGCAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCaaaa                   >  2:1826777/8‑98 (MQ=39)
                                                                                     tACCCTTTGCTGCAATTTGCATCAACAACCATGATAGTGAAATTCTTGTGATGTTGTTAACCAATTTCATAATTCGGGTTGACATGTCTTACCAAAAGGt                <  2:151081/100‑1 (MQ=255)
                                                                                      ggttgttGCTGCAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTa               >  1:3044300/6‑100 (MQ=255)
                                                                                      ggttgttGCTGCAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTa               >  2:1172715/6‑100 (MQ=255)
                                                                                      aCCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTa               >  2:547578/1‑100 (MQ=255)
                                                                                       gttgttGCTGCAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCaaaaag                 >  2:297768/5‑95 (MQ=255)
                                                                                       gttgttGCTGCAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAg              >  1:2015031/5‑100 (MQ=255)
                                                                                       gttgttGCTGCAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAg              >  1:3370248/5‑100 (MQ=255)
                                                                                       gccTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAg              >  2:2153610/2‑100 (MQ=255)
                                                                                         tgttGCTGCAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGaa            >  1:2603760/3‑100 (MQ=255)
                                                                                         cTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGaa            <  2:3288785/100‑1 (MQ=255)
                                                                                          tttGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAAc           >  1:290791/1‑100 (MQ=255)
                                                                                          tttGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAAc           <  2:934736/100‑1 (MQ=255)
                                                                                          gttGCTGCAAATTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAAc           >  1:3005575/2‑100 (MQ=255)
                                                                                           ttGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACt          <  1:129619/100‑1 (MQ=255)
                                                                                             gctgcAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTta        >  2:1283548/1‑100 (MQ=255)
                                                                                              ctgcAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTtat       >  1:1428579/1‑100 (MQ=255)
                                                                                               tgcAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTtata      >  2:3145702/1‑100 (MQ=255)
                                                                                                gcAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATAc     <  2:1775845/100‑1 (MQ=255)
                                                                                                  aaTTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGc   <  1:1455812/100‑1 (MQ=255)
                                                                                                   aTTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGcc  <  2:1172413/100‑1 (MQ=255)
                                                                                                    |                                                                                                  
GTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATCGATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCC  >  NC_000913/4551323‑4551520

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: