Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,993,889 T→A 15.2% intergenic (+33/‑50) ygeI → / → pbl uncharacterized protein/pseudogene, peptidoglycan‑binding enzyme family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,993,8890TA15.2% 286.5 / 41.2 112intergenic (+33/‑50)ygeI/pbluncharacterized protein/pseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base T (43/52);  new base A (12/5);  total (55/57)
Fisher's exact test for biased strand distribution p-value = 6.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

Only 100 of 118 total aligned reads displayed.

GCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATT  >  NC_000913/2993790‑2993988
                                                                                                   |                                                                                                   
gcCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtt                                                                                                     >  2:1844941/1‑100 (MQ=255)
 ccTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACttt                                                                                                    <  1:1167490/100‑1 (MQ=255)
  cTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtttt                                                                                                   >  2:107950/1‑100 (MQ=255)
  cTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtttt                                                                                                   >  1:1417513/1‑100 (MQ=255)
    tttCAAAATGATATGTATTCAGAAGATAATCCTATCTCTCCTTACTATAAAATTGAGTTGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtat                                                                                                 <  1:1588976/100‑1 (MQ=255)
    tttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtat                                                                                                 >  1:342963/1‑100 (MQ=255)
     ttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTGACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtata                                                                                                <  2:168780/100‑1 (MQ=255)
      tCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataa                                                                                               >  2:1533052/1‑100 (MQ=255)
       cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                              <  2:1511905/100‑1 (MQ=255)
       cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                              >  1:2126168/1‑100 (MQ=255)
       cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                              <  2:450086/100‑1 (MQ=255)
       cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                              >  2:1567508/1‑100 (MQ=255)
        aaaaTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGCCGATTGACCAAATAAATCATAACGAGATACTTTTataata                                                                                             >  2:2472109/1‑100 (MQ=255)
           aTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa                                                                                          <  2:2310996/100‑1 (MQ=255)
           aTGATATGTAATCAGAATATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGt                                                                                           <  1:2455465/99‑1 (MQ=255)
               tgtGTATTCAGAAAATAATCCTATCTCACCTTACTATAAAATCGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtctc                                                                                      >  2:1968689/3‑100 (MQ=255)
               tatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtctc                                                                                      <  2:1516673/100‑1 (MQ=255)
                atGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCg                                                                                     <  1:737692/100‑1 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                    <  2:2585026/100‑1 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                    >  2:2548516/1‑100 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGGATCTCGt                                                                                    <  2:1525804/100‑1 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                    <  1:866633/100‑1 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTAGTATACGATTAACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                    <  1:2779548/100‑1 (MQ=38)
                 tGTATGCAGAAGATAATCCTATCTCATCTTGCTATAAAATTGAGTGGTAGGCGATTGACCAATTAAAGCATAACGAGTTACTTTTATAATAGTATCTCg                                                                                     <  1:1576937/99‑1 (MQ=255)
                  gTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGtt                                                                                   >  1:685147/1‑100 (MQ=255)
                    aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTcc                                                                                 >  1:2272918/1‑100 (MQ=255)
                    aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttat                                                                                 <  2:1895672/100‑3 (MQ=255)
                    aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttat                                                                                 <  1:1013288/100‑3 (MQ=255)
                      tCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCtt                                                                               >  1:2786076/1‑100 (MQ=255)
                       cAGAAGATACTCCTATCTCAACTTACTATAAAATTGAGTAGCAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                              <  1:1008136/100‑1 (MQ=255)
                       cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                              >  2:2004926/1‑100 (MQ=255)
                       cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCtt                                                                               <  1:1118031/99‑1 (MQ=255)
                       cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTAGTTCCttt                                                                              >  1:1032819/1‑100 (MQ=255)
                       cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatgat                                                                              <  1:1787594/100‑6 (MQ=38)
                       cAGAAGATAATCCTATCTCACCTTACTATAAAAATGAGTGGTAGACGATTGACCAAATAAATCATAACGACATACTTTTATAATAGTATCTCGTTCCttt                                                                              >  1:1753848/1‑100 (MQ=255)
                        agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                             <  2:767423/100‑1 (MQ=255)
                        agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                             <  2:835329/100‑1 (MQ=255)
                         gaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttca                                                                            <  1:2597551/100‑1 (MQ=255)
                           agaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                          >  1:1790309/1‑100 (MQ=255)
                           agaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                          >  2:2308257/1‑100 (MQ=255)
                              tAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCaa                                                                       <  2:2258232/100‑1 (MQ=255)
                                 tCCTATCTCACCTTACTATAAAATTGAGTGGTAGACTATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                    <  1:219038/100‑1 (MQ=255)
                                 tCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                    <  1:665893/100‑1 (MQ=255)
                                 tCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATc                                                                     <  2:1582280/99‑1 (MQ=255)
                                  ccTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCt                                                                   >  2:713121/1‑100 (MQ=255)
                                  ccTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGAACAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCt                                                                   >  1:1594106/1‑100 (MQ=255)
                                    tATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTca                                                                 <  1:38303/100‑1 (MQ=255)
                                     aTCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcac                                                                <  2:731007/100‑1 (MQ=255)
                                     aTCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcac                                                                <  1:2411051/100‑1 (MQ=255)
                                        tcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcacaca                                                             <  1:1119829/100‑1 (MQ=255)
                                              tACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaaa                                                       <  1:384939/100‑1 (MQ=255)
                                                  ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAata                                                   >  1:113814/1‑100 (MQ=255)
                                                  ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAata                                                   >  1:2891903/1‑100 (MQ=255)
                                                  ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAat                                                    >  1:1518925/1‑99 (MQ=255)
                                                    aaaaTTGAGTGGTAGACGATTGACCAAATATATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATg                                                 <  1:809690/100‑1 (MQ=255)
                                                     aaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgc                                                >  1:1263689/1‑100 (MQ=255)
                                                     aaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgc                                                <  2:1763186/100‑1 (MQ=255)
                                                        ttGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCt                                             <  2:2712356/100‑1 (MQ=255)
                                                              tgTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaa                                       >  1:1301214/2‑100 (MQ=255)
                                                              ggTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaa                                       >  2:178244/1‑100 (MQ=255)
                                                                 aGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                    <  1:1975302/100‑1 (MQ=255)
                                                                  gACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGt                                   >  2:2463246/1‑100 (MQ=255)
                                                                    cGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTa                                 <  2:1611471/100‑1 (MQ=255)
                                                                         gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                            >  2:1783878/1‑100 (MQ=255)
                                                                          aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCt                           <  1:1640903/100‑1 (MQ=255)
                                                                          aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCt                           >  1:1105989/1‑100 (MQ=255)
                                                                            cAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTa                         >  1:48059/1‑100 (MQ=255)
                                                                              gaatgaaaggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                       >  2:1254655/11‑100 (MQ=25)
                                                                              aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                       <  2:1594367/100‑1 (MQ=255)
                                                                                taaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCa                     <  1:793367/100‑1 (MQ=255)
                                                                                 tgaaaggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCAt                    >  1:2423724/8‑100 (MQ=37)
                                                                                 aaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACGTTCAAACATTTGACATCACCAAAAGTTAAATTCTTAACCa                     >  1:227474/1‑99 (MQ=255)
                                                                                 aaaaaaTAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAATTATGCGCTCACCAAAAGTTAAATTCTTAAcc                      <  1:353523/93‑1 (MQ=255)
                                                                                  gaaaggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATc                   >  2:1094386/7‑100 (MQ=37)
                                                                                  aatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATc                   >  2:262872/1‑100 (MQ=255)
                                                                                     cATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTc                >  1:486289/1‑100 (MQ=255)
                                                                                      aTAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa               <  2:224581/100‑1 (MQ=255)
                                                                                       tAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  1:2792470/1‑100 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAGGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  1:970110/2‑100 (MQ=39)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa               >  1:789887/2‑99 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  2:1997083/2‑100 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  1:1881630/2‑100 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  2:433011/2‑100 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTAAACCATCTTCAc              >  1:2270756/2‑100 (MQ=38)
                                                                                         aCGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt            >  2:2536084/1‑100 (MQ=255)
                                                                                         aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt            >  1:1166973/1‑100 (MQ=255)
                                                                                             gaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg        <  1:1625092/100‑1 (MQ=255)
                                                                                             gaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg        <  2:347998/100‑1 (MQ=255)
                                                                                             gaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg        <  2:407489/100‑1 (MQ=255)
                                                                                              aTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGt       <  2:1949192/100‑1 (MQ=255)
                                                                                               tACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAGAAGTTAAATTCTTAACCATCTTCACGTTTTGta      <  1:1511274/100‑1 (MQ=255)
                                                                                               tACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGta      <  2:2355027/100‑1 (MQ=255)
                                                                                               tACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGta      <  2:1790317/100‑1 (MQ=255)
                                                                                               tACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGta      <  2:874841/100‑1 (MQ=255)
                                                                                               tACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGt       >  1:1709125/1‑99 (MQ=255)
                                                                                                aCTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtat     <  2:1569949/100‑1 (MQ=255)
                                                                                                 cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtata    >  1:2177765/1‑100 (MQ=255)
                                                                                                 cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtata    >  2:1373152/1‑100 (MQ=255)
                                                                                                   tttATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatt  >  2:641133/1‑100 (MQ=255)
                                                                                                   tttATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatt  <  2:1084490/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATT  >  NC_000913/2993790‑2993988

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: