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.

GTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCC  >  NC_000913/4551323‑4551520
                                                                                                    |                                                                                                  
gttgttATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAAt                                                                                                    <  1:1371022/100‑1 (MQ=255)
gttgttAGTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAAt                                                                                                    <  1:1216930/100‑1 (MQ=255)
   gttATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTg                                                                                                 >  2:1913006/1‑100 (MQ=255)
   gttATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTg                                                                                                 >  2:1063893/1‑100 (MQ=255)
      aTTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTgcag                                                                                              <  2:2847725/100‑1 (MQ=255)
       ttttGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTgcagc                                                                                             <  1:3325606/100‑1 (MQ=255)
       ttttGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTgcagc                                                                                             >  1:283253/1‑100 (MQ=255)
       ttttGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTgcagc                                                                                             <  1:2059198/100‑1 (MQ=255)
          taggagaaatgttatgtAAATGGGTCAACCATGGTTGCGATGAATGTTACATCCGCTGATCAATAACCATC‑GATTACACTTTGCTGCAAATTgaagc                                                                                             <  2:3098161/89‑1 (MQ=255)
          tGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAAATTGCAGcaaa                                                                                          <  2:1810318/100‑2 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                         <  1:1149669/100‑3 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                         <  2:2858069/100‑3 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                         <  2:2427940/100‑3 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                         <  1:3494020/100‑3 (MQ=255)
           gggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAAATTGCAGcaaag                                                                                         <  1:2667380/100‑3 (MQ=255)
            ggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGcaacaa                                                                                        >  1:1767608/1‑100 (MQ=255)
            ggCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAAATTGCAGcaaagg                                                                                        <  1:579079/100‑4 (MQ=255)
             gCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAAATTGCAGcaaaggg                                                                                       <  1:1634675/100‑5 (MQ=255)
              cGATATGTTATGTAAATTGGTCAACAATTGTTGCGATGAATGTCACATCCTCTGTTCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGcaacaa                                                                                        >  2:2310575/1‑98 (MQ=255)
               gATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAAcc                                                                                      <  2:1430778/99‑1 (MQ=255)
                atatGTTCTGTAGATCGGTCAACCCTTGTTGCGATGCATGTCCCATCCTCTGAGCAAGAACCCTC‑GATTCCCCTTTGCTGCAATTTGCATCAACAAcc                                                                                      <  2:3547862/98‑1 (MQ=255)
                atatGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCAt                                                                                    <  2:1409606/100‑1 (MQ=255)
                 tatgttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATg                                                                                   >  1:940908/1‑100 (MQ=255)
                 tatgttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAAATTGCAGcaaagggtaat                                                                                   <  1:3610911/100‑9 (MQ=39)
                    gttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATgaga                                                                                >  2:1712478/1‑100 (MQ=255)
                    gttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATgaga                                                                                >  2:1134871/1‑100 (MQ=255)
                     ttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGagag                                                                               <  1:810817/100‑1 (MQ=255)
                     ttatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGagag                                                                               <  1:189646/100‑1 (MQ=255)
                      tatgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGt                                                                              >  2:361459/1‑100 (MQ=255)
                        tgtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGa                                                                            >  1:3288665/1‑100 (MQ=255)
                         gtAAATTGGTCAACCATTGTTGCGATGAATGTCACATCTTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGaa                                                                           >  2:1533694/1‑100 (MQ=255)
                          tAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGaaa                                                                          >  1:1907014/1‑100 (MQ=255)
                            aaTTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAAtt                                                                        >  1:1969858/1‑100 (MQ=255)
                                ggTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtg                                                                    <  2:2639458/100‑1 (MQ=255)
                                 gTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTTCCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtg                                                                    <  2:2020616/99‑1 (MQ=255)
                                 gTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtgt                                                                   >  2:1166318/1‑100 (MQ=255)
                                 gTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTtgt                                                                   >  1:608060/1‑100 (MQ=255)
                                   cAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGa                                                                 >  2:1598416/1‑100 (MQ=255)
                                   cAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGa                                                                 >  1:2692034/1‑100 (MQ=255)
                                   cAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGa                                                                 >  1:2357151/1‑100 (MQ=255)
                                        aTTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTgg                                                            >  1:3566358/1‑100 (MQ=255)
                                         ttgttgCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGt                                                           <  2:2745775/100‑1 (MQ=255)
                                          tgttgCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGtt                                                          >  1:2246149/1‑100 (MQ=255)
                                           gttgCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTa                                                         <  1:1134880/100‑1 (MQ=255)
                                           gttgCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTa                                                         <  2:831101/100‑1 (MQ=255)
                                             tgCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAAc                                                       <  2:2859924/100‑1 (MQ=255)
                                               cGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAAc                                                       >  2:1398125/1‑98 (MQ=255)
                                               cGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCa                                                     >  2:2379871/1‑100 (MQ=255)
                                                gATGAATGCCCCATACTCTGGACAATAA‑AATCTGAGTACCCTTTGCTACAATTTGCAGCAACAACCAAGAGAGTAAAATTAACGTGATGTGGAAAACCa                                                     >  2:3630767/1‑99 (MQ=255)
                                                    aaTGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTGACCAATTTc                                                >  2:749920/1‑100 (MQ=255)
                                                      tGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAg                                              <  1:540405/100‑1 (MQ=255)
                                                      tGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAg                                              <  1:1718213/100‑1 (MQ=255)
                                                      tGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAg                                              >  1:1516373/1‑100 (MQ=255)
                                                      tGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAg                                              <  2:1016145/100‑1 (MQ=255)
                                                       gTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGa                                             >  1:766366/1‑100 (MQ=255)
                                                         cacaTCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAg                                              <  2:1350976/97‑1 (MQ=255)
                                                         cacaTCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAAt                                           >  2:3504192/1‑100 (MQ=255)
                                                             tCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCgg                                       >  2:2145676/1‑100 (MQ=255)
                                                              ccTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCggg                                      <  1:2219921/100‑1 (MQ=255)
                                                               ctctGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGt                                     <  2:1019710/100‑1 (MQ=255)
                                                               ctctGATCAATAACCAGC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATCCTTGTGATGTGGTTAACCAATTTCAGAATTCggg                                      <  2:1282190/99‑1 (MQ=255)
                                                                nctGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGtt                                    >  1:3623016/2‑100 (MQ=255)
                                                                 ctGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTg                                   <  2:3013047/100‑1 (MQ=255)
                                                                 ctGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTg                                   <  1:1198741/100‑1 (MQ=255)
                                                                   gATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGAc                                 <  2:638572/100‑1 (MQ=255)
                                                                   gATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGAc                                 >  2:635608/1‑100 (MQ=255)
                                                                       aaTAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGt                             >  1:2042500/1‑100 (MQ=255)
                                                                        aTAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTc                            <  1:3176001/100‑1 (MQ=255)
                                                                         tAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATCCTTGTGATGTGGTTAACAAATTTCAGAATACGGGTTGa                                  >  2:3328352/1‑93 (MQ=255)
                                                                          aaCCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCtt                          <  1:2520154/100‑1 (MQ=255)
                                                                          aaCCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCt                           <  1:1576330/99‑1 (MQ=255)
                                                                           aCCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTa                         >  2:1008153/1‑100 (MQ=255)
                                                                            ccATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTAc                        >  2:1928693/1‑100 (MQ=255)
                                                                            ccATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTAc                        >  1:2383508/1‑100 (MQ=255)
                                                                             cATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTAcc                       <  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)
                                                                                                    |                                                                                                  
GTTGTTATTTTGGGCGATATGTTATGTAAATTGGTCAACCATTGTTGCGATGAATGTCACATCCTCTGATCAATAACCATC‑GATTACCCTTTGCTGCAATTTGCAGCAACAACCATGAGAGTGAAATTCTTGTGATGTGGTTAACCAATTTCAGAATTCGGGTTGACATGTCTTACCAAAAGGTAGAACTTATACGCC  >  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: ‑