Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,482,325 G→A R346H (CGC→CAC)  yheS → ABC‑F family protein predicted regulatory ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,482,3250GA99.2% 423.5 / ‑4.5 120R346H (CGC→CAC) yheSABC‑F family protein predicted regulatory ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (64/55);  total (64/56)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 120 total aligned reads displayed.

AGATGGAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGG  >  NC_000913/3482229‑3482421
                                                                                                |                                                                                                
aGATGGAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACaa                                                                                               >  1:1654658/1‑100 (MQ=255)
aGATGGAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACaa                                                                                               <  1:2132706/100‑1 (MQ=255)
aGATGGAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACaa                                                                                               <  1:1897128/100‑1 (MQ=255)
 gATGGAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGATCTGTTAGGCCACAAt                                                                                              <  2:1431438/100‑1 (MQ=255)
   tGGAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATgg                                                                                            >  1:783179/1‑100 (MQ=255)
   tGGAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATgg                                                                                            >  1:187530/1‑100 (MQ=255)
    ggAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGgc                                                                                           >  1:2630359/1‑100 (MQ=255)
    ggAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATCGGTCTGTTAGGCCACAATGgc                                                                                           >  1:3008694/1‑100 (MQ=255)
     gAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGgcg                                                                                          >  2:879705/1‑100 (MQ=255)
     gAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGgcg                                                                                          >  2:2945441/1‑100 (MQ=255)
     gAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGgcg                                                                                          >  2:694962/1‑100 (MQ=255)
     gAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGgcg                                                                                          >  2:695116/1‑100 (MQ=255)
       aaaaGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGcgcg                                                                                        <  2:1956245/100‑1 (MQ=255)
         aaGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCggg                                                                                      <  2:1964525/100‑1 (MQ=255)
         aaGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCggg                                                                                      <  2:448131/100‑1 (MQ=255)
           gTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTa                                                                                    <  1:3070670/100‑1 (MQ=255)
              aGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAAt                                                                                 >  1:1221706/1‑100 (MQ=255)
               gcgcGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATc                                                                                <  2:562283/100‑1 (MQ=255)
                  cgGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGAc                                                                             <  2:1666701/100‑1 (MQ=255)
                   gggCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACa                                                                            <  2:769977/100‑1 (MQ=255)
                   gggCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACa                                                                            <  2:1033323/100‑1 (MQ=255)
                       tATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTaa                                                                        <  1:1000468/100‑1 (MQ=255)
                        aTGGCGATCGCATTATTCTCGACTCGCTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAAt                                                                       >  2:2118787/1‑100 (MQ=255)
                        aTGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAAt                                                                       >  1:2504774/1‑100 (MQ=255)
                         tGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATc                                                                      <  1:2014050/100‑1 (MQ=255)
                          ggCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCACGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCa                                                                     >  1:2370583/1‑100 (MQ=255)
                           gCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCag                                                                    <  2:1594931/100‑2 (MQ=255)
                            cGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCaaa                                                                   >  2:866103/1‑100 (MQ=255)
                            cGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCaaa                                                                   >  1:3032679/1‑100 (MQ=255)
                            cGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTATGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCaaa                                                                   >  1:1662440/1‑100 (MQ=255)
                             gATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAAc                                                                  <  2:659734/100‑1 (MQ=255)
                             gATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAAc                                                                  <  2:659498/100‑1 (MQ=255)
                              aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACt                                                                 <  2:1630431/100‑1 (MQ=255)
                              aTCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACt                                                                 >  2:1660490/1‑100 (MQ=255)
                               tCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACtn                                                                <  2:306695/100‑2 (MQ=255)
                                cGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGCCGCGGGTGAATCGACATTCATCAAACTGt                                                               >  2:2769716/1‑100 (MQ=255)
                                 gCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGtt                                                              <  1:944510/100‑1 (MQ=255)
                                  cATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTa                                                             >  2:2574478/1‑100 (MQ=255)
                                  cATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTa                                                             >  1:1262314/1‑100 (MQ=255)
                                     tattCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGcc                                                          <  2:2541670/100‑1 (MQ=255)
                                       ttCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTCAATCGACATTAATCAAACTGTTAGCCgg                                                        <  1:145107/100‑1 (MQ=255)
                                       ttCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCgg                                                        <  1:2638770/100‑1 (MQ=255)
                                        tctcGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGt                                                       <  1:2706422/100‑1 (MQ=255)
                                        tctcGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGt                                                       >  2:623081/1‑100 (MQ=255)
                                         ctcgactcgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTg                                                      >  1:2239078/1‑100 (MQ=255)
                                               tcgaTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTg                                                <  1:256186/100‑1 (MQ=255)
                                                  aTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTgcgc                                             <  2:1327086/100‑1 (MQ=255)
                                                  aTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTgcgc                                             <  2:1327067/100‑1 (MQ=255)
                                                  aTTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTgcgc                                             >  1:417128/1‑100 (MQ=255)
                                                   ttAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGcc                                            <  2:1008553/100‑1 (MQ=255)
                                                   ttAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGcc                                            <  2:17727/100‑1 (MQ=255)
                                                      aaCTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCACACTGTTAGCCGGCGAACTTGCGCCAGt                                         >  2:1938686/1‑100 (MQ=255)
                                                      aaCTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGt                                         >  1:1049412/1‑100 (MQ=255)
                                                      aaCTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGt                                         >  1:1724282/1‑100 (MQ=255)
                                                      aaCTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGt                                         >  2:1939244/1‑100 (MQ=255)
                                                        cTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCa                                       >  1:781257/1‑100 (MQ=255)
                                                        cTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCa                                       >  1:781155/1‑100 (MQ=255)
                                                         tGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAg                                      >  2:486442/1‑100 (MQ=255)
                                                          ggaCCTGGTGCCCGGCTCGCGTATTGGTCTGTTGGGCCACGATGGCGCGGGTGAATCGACATTAGTCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGc                                     <  2:2568818/98‑1 (MQ=255)
                                                          gAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGc                                     <  2:3008733/100‑1 (MQ=255)
                                                          gAACCTGGTGCCCGGCTCGCGTATCGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGc                                     <  2:3008694/100‑1 (MQ=255)
                                                            aCCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCgg                                   <  2:799010/100‑1 (MQ=255)
                                                              cTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTg                                 >  2:2024624/1‑100 (MQ=255)
                                                               tGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGa                                >  2:510833/1‑100 (MQ=255)
                                                                 gTGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGaaa                              >  1:1727074/1‑100 (MQ=255)
                                                                  tGCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAAt                             >  1:2905398/1‑100 (MQ=255)
                                                                   gCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAAtt                            <  1:680138/100‑1 (MQ=255)
                                                                   gCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAAtt                            >  1:2064042/1‑100 (MQ=255)
                                                                   gCCCGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAAtt                            >  2:2630325/1‑100 (MQ=255)
                                                                    cccGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTg                           >  2:239754/1‑100 (MQ=255)
                                                                    cccGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGCTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTg                           >  2:2278196/1‑100 (MQ=255)
                                                                    cccGGCTCGCGTATTCGTCTCCTAGGCCACAATGGCGCGGGTAAATCGACCTTAACCAAACTGTTAGCCGGTGAACTTGCGCCGGTCATCGGTGAAATTg                           >  2:2278102/1‑100 (MQ=255)
                                                                     ccGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTgg                          <  1:3050723/100‑1 (MQ=255)
                                                                     ccGGCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTgg                          >  2:2765981/1‑100 (MQ=255)
                                                                        gCTCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCt                       >  2:170843/1‑100 (MQ=255)
                                                                          tCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTgg                     >  2:621657/1‑100 (MQ=255)
                                                                          tCGCGTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTgg                     <  1:2877019/100‑1 (MQ=255)
                                                                           cgcgTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGc                    <  1:2574478/100‑1 (MQ=255)
                                                                           cgcgTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGc                    <  1:2875205/100‑1 (MQ=255)
                                                                            gcgTATTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCg                   <  2:2239078/100‑1 (MQ=255)
                                                                                aTTGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAg               <  1:962216/100‑1 (MQ=255)
                                                                                 ttGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAgg              <  1:2026881/100‑1 (MQ=255)
                                                                                  tGGTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAggg             <  1:2687051/100‑1 (MQ=255)
                                                                                    gTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGAt           <  2:819787/100‑1 (MQ=255)
                                                                                    gTCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGAt           >  2:1823086/1‑100 (MQ=255)
                                                                                     tCTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATc          <  2:2421156/100‑1 (MQ=255)
                                                                                      cTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCa         <  1:267613/100‑1 (MQ=255)
                                                                                      cTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCa         >  2:1081933/1‑100 (MQ=255)
                                                                                      cTGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCa         >  2:2910682/1‑100 (MQ=255)
                                                                                       tGTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCaa        >  2:1906934/1‑100 (MQ=255)
                                                                                        gTTAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCaaa       >  1:2612525/1‑100 (MQ=255)
                                                                                         ttAGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAAc      >  1:1948494/1‑100 (MQ=255)
                                                                                          tAGGCCACAATTGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACt     <  1:2995545/100‑1 (MQ=255)
                                                                                           aGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTc    >  2:492033/1‑100 (MQ=255)
                                                                                           aGGCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTc    >  1:640904/1‑100 (MQ=255)
                                                                                            ggCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCg   <  1:2383878/100‑1 (MQ=255)
                                                                                            ggCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCg   >  1:658277/1‑100 (MQ=255)
                                                                                             gCCACAATGGCGCGGGTAACTCGACACTAATCGACCTGTTCGCGGGTGAACTTGCGCCAGTCAGCGGGGAAATTGGTCTGGCGAAAGGGATCAAACTCgg  <  2:2607666/100‑1 (MQ=255)
                                                                                             gCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCgg  >  1:1162199/1‑100 (MQ=255)
                                                                                             gCCACAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGAGCAAACTCgg  >  1:2667607/1‑100 (MQ=255)
                                                                                                |                                                                                                
AGATGGAAAAAGTCAGCGCGGGCTATGGCGATCGCATTATTCTCGACTCGATTAAACTGAACCTGGTGCCCGGCTCGCGTATTGGTCTGTTAGGCCGCAATGGCGCGGGTAAATCGACATTAATCAAACTGTTAGCCGGTGAACTTGCGCCAGTCAGCGGTGAAATTGGTCTGGCGAAAGGGATCAAACTCGG  >  NC_000913/3482229‑3482421

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