Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,310,728 T→C 53.6% I85V (ATT→GTT)  alsA ← D‑allose ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,310,7280TC53.6% 33.8 / 183.3 141I85V (ATT→GTT) alsAD‑allose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (32/33);  new base C (38/37);  total (71/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 146 total aligned reads displayed.

TAATATTGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTATTCCGGA  >  NC_000913/4310632‑4310827
                                                                                                |                                                                                                   
taataTTGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCATTCCGGTTAATTCATCAATAACGCTGAGTTCCTGATAaacaat                                                                                                  >  1:1067213/1‑100 (MQ=255)
  atatTGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATcc                                                                                                <  1:684625/100‑1 (MQ=255)
    atTGACGCCACAGCTTTTTTTCGTCAGATTACGACCAATATATAAATTTACCAGTCCGGTTAATTCATCAATAACGCTGAATTCCTGATAAACAATcc                                                                                                >  1:3228803/1‑98 (MQ=255)
    atTGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCg                                                                                              >  2:2014507/1‑100 (MQ=255)
    atTGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCg                                                                                              <  2:3721921/100‑1 (MQ=255)
     ttGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGa                                                                                             >  2:1997566/1‑100 (MQ=255)
      tGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGAt                                                                                            <  1:3745280/100‑1 (MQ=255)
       gACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGAt                                                                                            >  1:3770959/1‑99 (MQ=255)
       gACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATa                                                                                           >  2:79109/1‑100 (MQ=255)
          gCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATAcc                                                                                         >  2:4237695/1‑99 (MQ=255)
            cacaGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAg                                                                                      >  2:4102982/1‑100 (MQ=255)
            cacaGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAg                                                                                      >  1:552799/1‑100 (MQ=255)
                gATTTTTTTCGTCAGATGACGCCCAATATATCCATTTTCCAGTACGGTTACTTCATCAATAACGCTGAGTTCCTGATAAACAATCCCTATACCTAGTtg                                                                                   <  1:3998443/99‑1 (MQ=255)
                gATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATCCCGAGtt                                                                                    >  1:3692364/1‑98 (MQ=255)
                gATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTtgt                                                                                  <  1:1575772/100‑1 (MQ=255)
                     ttttCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGc                                                                             >  1:784125/1‑100 (MQ=255)
                       ttCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCgn                                                                             <  1:3500428/98‑2 (MQ=255)
                          gTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATCCCGAGTTGTGCGGCTAAt                                                                         >  2:3046319/1‑99 (MQ=255)
                          gTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAAtt                                                                        >  2:67985/1‑100 (MQ=255)
                          gTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAAtt                                                                        >  1:3133475/1‑100 (MQ=255)
                           tCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAAttt                                                                       <  1:1766776/100‑1 (MQ=255)
                            cAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTa                                                                      >  2:3796855/1‑100 (MQ=255)
                            cAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTa                                                                      >  2:3982880/1‑100 (MQ=255)
                            cAGATGACGACCAATATATAAATTTTCCAGTACGGCTAATTCACCAATAACGCTTACTTCCTGATAAATAATCCACATACCAACTTGTAGCGCTAAttac                                                                      >  1:4178131/1‑98 (MQ=255)
                             aGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTa                                                                      <  2:3352024/99‑1 (MQ=255)
                             aGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTAt                                                                     <  2:2536095/100‑1 (MQ=255)
                                tGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGAt                                                                  <  2:3376472/100‑1 (MQ=255)
                                tGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGAt                                                                  >  1:19789/1‑100 (MQ=255)
                                tGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGAt                                                                  <  2:3341792/100‑1 (MQ=255)
                                  acgacCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATcc                                                                <  1:1685837/100‑1 (MQ=255)
                                    gacCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAg                                                              <  1:3987219/100‑1 (MQ=255)
                                     acCAATATATAAATTTTCCAGTACGGTTAATTCATCAAGAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGc                                                             <  1:1923551/100‑1 (MQ=255)
                                      ccAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCt                                                            >  1:2201722/1‑100 (MQ=255)
                                        aaTATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTg                                                          >  1:3373986/1‑100 (MQ=255)
                                        aaTATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTg                                                          <  2:1589772/100‑1 (MQ=255)
                                        aaTATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTg                                                          <  1:79109/100‑1 (MQ=255)
                                         atatatAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGt                                                         <  1:2350550/100‑1 (MQ=255)
                                          tatataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGtt                                                        <  1:4218879/100‑1 (MQ=255)
                                            tataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGCGTTGTGCCGCTAATTTATGATCCAGCTTGCta                                                       >  1:1302797/1‑99 (MQ=255)
                                            tataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTtat                                                      >  2:1292366/1‑100 (MQ=255)
                                            tataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATCCCGAGTTGTGCCGCAAATTTATGATcc                                                                >  1:1115685/1‑90 (MQ=255)
                                             ataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTtata                                                     <  2:3119517/100‑1 (MQ=255)
                                             ataAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTtata                                                     >  2:3031819/1‑100 (MQ=255)
                                              taAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAg                                                    >  2:1484427/1‑100 (MQ=255)
                                              taAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAg                                                    >  1:1464938/1‑100 (MQ=255)
                                              taAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAg                                                    >  2:272112/1‑100 (MQ=255)
                                               aaaTTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGc                                                   >  2:942890/1‑100 (MQ=255)
                                               aaaTTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGc                                                   <  2:583442/100‑1 (MQ=255)
                                               aaaTTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGc                                                   <  1:14608/100‑1 (MQ=255)
                                                 aTTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTa                                                 >  1:2928986/1‑100 (MQ=255)
                                                 aTTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTa                                                 >  2:497934/1‑100 (MQ=255)
                                                 aTTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTa                                                 <  1:1331485/100‑1 (MQ=255)
                                                   tttCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTTATAAACAATCCCGATACCGAGGTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAAt                                               <  1:3005389/100‑1 (MQ=255)
                                                   tttCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAAt                                               >  1:2667156/1‑100 (MQ=255)
                                                    ttCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATg                                              >  2:1711709/1‑100 (MQ=255)
                                                      ccAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGtt                                            >  2:3820964/1‑100 (MQ=255)
                                                       cAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGtta                                           >  1:2203297/1‑100 (MQ=255)
                                                        aGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGttat                                          <  1:1930587/100‑1 (MQ=255)
                                                        aGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGttat                                          >  1:724810/1‑100 (MQ=255)
                                                         gTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGttatt                                         >  2:215362/1‑100 (MQ=255)
                                                          tACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGttatta                                        >  2:2012571/1‑100 (MQ=255)
                                                            cGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaat                                      <  1:3308575/100‑1 (MQ=255)
                                                            cGGTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaat                                      >  2:461586/1‑100 (MQ=255)
                                                             ggTTAATTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaatg                                     >  1:985280/1‑100 (MQ=255)
                                                                 aaTTCATCAATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaatggtaa                                 >  1:666744/1‑100 (MQ=255)
                                                                   ttCATCAATAACGCTGAGTTCCTTATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaatggtaatg                               <  1:134669/100‑1 (MQ=255)
                                                                    tcatcaATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaatggtaatgg                              <  2:4175862/100‑1 (MQ=255)
                                                                    tcatcaATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaatggtaatgg                              >  1:3206327/1‑100 (MQ=255)
                                                                    tcatcaATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATtaatggtaatgg                              >  2:3480851/1‑100 (MQ=255)
                                                                      atcaATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTg                            >  2:703906/1‑100 (MQ=255)
                                                                        caATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGcc                          <  2:2859200/100‑1 (MQ=255)
                                                                        caATAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGcc                          <  2:1749578/100‑1 (MQ=255)
                                                                         aaTAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCt                         >  2:1074918/1‑100 (MQ=255)
                                                                         aaTAACGCTGAGTTCCTGATAAAAAATCCCGATACCGAGATGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTGATTaatggtaatggt                             >  1:4146473/1‑96 (MQ=255)
                                                                          aTAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCtt                        >  2:2209345/1‑100 (MQ=255)
                                                                           tAACGCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCttt                       <  1:648025/100‑1 (MQ=255)
                                                                               gCTGAGTTCCTGATATACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTc                   <  2:4214169/100‑1 (MQ=255)
                                                                               gCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTc                   <  2:814326/100‑1 (MQ=255)
                                                                               gCTGAGTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTc                   <  1:865386/100‑1 (MQ=255)
                                                                                 tGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCgg                 >  2:2471925/1‑100 (MQ=255)
                                                                                  gAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGc                <  2:4178153/100‑1 (MQ=255)
                                                                                  gAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGc                <  2:2386552/100‑1 (MQ=255)
                                                                                  gAGGTCCTGATAAATAATCCGGATACCGAGTTGTGCCGGTAATGTATGATCCAGCTTGTTATAGCTAATGGTATTAATGGTAATGGTGCCTTTGGTCGGc                <  1:3015641/100‑1 (MQ=255)
                                                                                    gTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGGTCCAGCTTGTTATAGCTGATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTc              <  1:4070177/100‑1 (MQ=255)
                                                                                    gTTCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGc                >  1:2318567/1‑98 (MQ=255)
                                                                                      tCCTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCAt            >  2:3922552/1‑100 (MQ=255)
                                                                                       ccTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATg           <  1:2729269/100‑1 (MQ=255)
                                                                                       ccTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATg           <  1:2862295/100‑1 (MQ=255)
                                                                                       ccTGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATg           >  1:3297325/1‑100 (MQ=255)
                                                                                       ccTGATAAACAATCCCGATACAGAGGTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATg           >  1:3585396/1‑100 (MQ=255)
                                                                                        cTGATAAATAATCCCGATACCGAGTTTTGCCGCTAATTTATGTTCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATg           <  1:1845516/99‑1 (MQ=255)
                                                                                        cTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGt          >  1:482620/1‑100 (MQ=255)
                                                                                        cTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGt          >  2:1339130/1‑100 (MQ=255)
                                                                                         tGATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTa         >  1:3872909/1‑100 (MQ=255)
                                                                                          gATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTAt        >  2:2218918/1‑100 (MQ=255)
                                                                                          gATAAACAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGATAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATg           <  1:3869746/97‑1 (MQ=255)
                                                                                           ataaacaaTCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTAtt       >  2:3989993/1‑100 (MQ=255)
                                                                                            taaataaTCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTATTc      >  2:2819408/1‑100 (MQ=255)
                                                                                               ataatCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTATTCCgg   >  2:2263478/1‑100 (MQ=255)
                                                                                                taatCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTATTCCGGa  <  1:448037/100‑1 (MQ=255)
                                                                                                |                                                                                                   
TAATATTGACGCCACAGATTTTTTTCGTCAGATGACGACCAATATATAAATTTTCCAGTACGGTTAATTCATCAATAACGCTGAGTTCCTGATAAATAATCCCGATACCGAGTTGTGCCGCTAATTTATGATCCAGCTTGTTATAGCTAATGTTATTAATGGTAATGGTGCCTTTGGTCGGCTCATGTATTCCGGA  >  NC_000913/4310632‑4310827

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: