Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,993,889 T→A 19.6% 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,8890TA19.6% 431.3 / 106.4 197intergenic (+33/‑50)ygeI/pbluncharacterized protein/pseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base T (82/75);  new base A (24/14);  total (108/89)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-02

Only 100 of 201 total aligned reads displayed.

CCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATT  >  NC_000913/2993791‑2993988
                                                                                                  |                                                                                                    
ccTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACttt                                                                                                     <  1:999342/100‑1 (MQ=255)
    ttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtata                                                                                                 <  2:3227707/100‑1 (MQ=255)
    ttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtata                                                                                                 <  2:438815/100‑1 (MQ=255)
    ttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtata                                                                                                 >  1:1647875/1‑100 (MQ=255)
      cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                               <  1:77533/100‑1 (MQ=255)
      cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                               <  2:282258/100‑1 (MQ=255)
       aaaaTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataata                                                                                              >  2:1249642/1‑100 (MQ=255)
        aaaTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                                                                                             <  2:1208437/100‑1 (MQ=255)
          aTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa                                                                                           >  1:1690457/1‑100 (MQ=255)
           tGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAt                                                                                          >  2:3474979/1‑100 (MQ=255)
           tGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAt                                                                                          >  1:909063/1‑100 (MQ=255)
            gATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc                                                                                         >  1:3950348/1‑100 (MQ=255)
               atGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCg                                                                                      <  1:1453684/100‑1 (MQ=255)
                tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                     >  1:1292403/1‑100 (MQ=255)
                tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                     <  2:16987/100‑1 (MQ=255)
                tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                     <  1:2829705/100‑1 (MQ=255)
                tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                     <  1:3614332/100‑1 (MQ=255)
                tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                     <  1:2183957/100‑1 (MQ=255)
                tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAAAAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                     >  1:2766201/1‑100 (MQ=255)
                  tATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTc                                                                                   >  2:2654795/1‑100 (MQ=255)
                  tATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGtta                                                                                   <  1:2264530/100‑2 (MQ=255)
                  tATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGtta                                                                                   <  1:1039363/100‑2 (MQ=255)
                  tATTCAGAAGATAAACCTATCTTACCTTACAATAAAATTGAGAGGTAGACGATTGACCAAATAAAACAAAACGAGATACTTTTATAATAGTATCACGtt                                                                                    >  2:3285209/1‑99 (MQ=255)
                   aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTcc                                                                                  >  1:491228/1‑100 (MQ=255)
                   aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttat                                                                                  <  1:3606246/100‑3 (MQ=255)
                   aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttat                                                                                  <  2:1984601/100‑3 (MQ=255)
                    ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAATTCATAACGAGATACTATTATAAAAGTATCTCGttatg                                                                                 <  2:2108496/100‑4 (MQ=38)
                       agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                              <  1:215218/100‑1 (MQ=255)
                       agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                              <  2:2942349/100‑1 (MQ=255)
                       agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                              >  2:3330453/1‑100 (MQ=255)
                         aagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcat                                                                            <  2:1473912/100‑1 (MQ=255)
                           gaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattc                                                                          <  1:3464869/100‑1 (MQ=255)
                            aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattca                                                                         <  2:606397/100‑1 (MQ=255)
                            aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattca                                                                         <  2:1582994/100‑1 (MQ=255)
                            aTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcattca                                                                         >  1:1078905/1‑100 (MQ=255)
                             tAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCaa                                                                        >  2:3677403/1‑100 (MQ=255)
                              aaTCCTATCTTACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAAt                                                                       <  1:116627/100‑1 (MQ=255)
                               aTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATc                                                                      <  2:549051/100‑1 (MQ=255)
                                tCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                     <  2:2523262/100‑1 (MQ=255)
                                tCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                     <  1:1548332/100‑1 (MQ=255)
                                 ccTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCt                                                                    <  1:3443613/100‑1 (MQ=255)
                                 ccTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCt                                                                    >  1:2379134/1‑100 (MQ=255)
                                 ccTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCt                                                                    <  2:2871592/100‑1 (MQ=255)
                                 ccTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCt                                                                    <  2:2889476/100‑1 (MQ=255)
                                    aTCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcac                                                                 <  1:978614/100‑1 (MQ=255)
                                    aTCTAAACTTACTATAAAATAAAGTGGAAGACGATTGACCAAATAAATCATAACCAGATTCTCTTCTAATAGTATCTCGATCATTTCAATCAATcc                                                                     >  2:3496080/1‑96 (MQ=255)
                                        cACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCTTCACACAt                                                             >  2:3461957/1‑100 (MQ=255)
                                         aCCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATg                                                            <  2:231222/100‑1 (MQ=255)
                                          ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGa                                                          >  1:1499058/1‑100 (MQ=255)
                                          cattACTATAAAACTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGa                                                          >  1:2676468/3‑100 (MQ=255)
                                           cTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaa                                                         <  1:273693/100‑1 (MQ=255)
                                            ttACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaa                                                        >  2:3386422/1‑100 (MQ=255)
                                              aCTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaaaa                                                      <  1:2352254/100‑1 (MQ=255)
                                               cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaaaaa                                                     <  1:3475004/100‑1 (MQ=255)
                                               cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaaaaa                                                     >  1:3315742/1‑100 (MQ=255)
                                                tataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAat                                                    >  2:3940566/1‑100 (MQ=255)
                                                 ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAata                                                   >  1:1783147/1‑100 (MQ=255)
                                                   aaaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATg                                                 >  2:2754962/1‑100 (MQ=255)
                                                   aaaaTTAAGTGGTAGACGATTAACAAAATAAATCATAACGAGATAGTTTTATAATAGTATCTCGTTCTTTTACTTCAATTCTCATAAATGAAAAAATATg                                                 >  2:1906480/1‑100 (MQ=255)
                                                     aaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgcg                                               <  1:3677431/100‑1 (MQ=255)
                                                      aTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgcgc                                              >  2:2386889/1‑100 (MQ=255)
                                                           gTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCAcc                                         >  1:3002206/1‑100 (MQ=255)
                                                            tGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCa                                        <  1:2445378/100‑1 (MQ=255)
                                                              gTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaaa                                      <  1:3084206/100‑1 (MQ=255)
                                                               tAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaaaa                                     >  1:937873/1‑100 (MQ=255)
                                                                aGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                    <  2:3854836/100‑1 (MQ=255)
                                                                  aCGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaaa                                      <  2:225010/96‑1 (MQ=255)
                                                                   cGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTACTTTCATTAAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTa                                 >  2:1922760/1‑100 (MQ=255)
                                                                    gATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATACTCACAAATGAAAAAATTTGTGTTCACAAAAAGTTaa                                >  2:964142/1‑100 (MQ=255)
                                                                      ttGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAAt                              <  2:1783138/100‑1 (MQ=255)
                                                                      ttGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAAt                              <  2:872194/100‑1 (MQ=255)
                                                                        gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                            >  1:2499611/1‑100 (MQ=255)
                                                                          ccAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCGCCAAAAGTTAAATTCtt                          >  1:2215760/1‑100 (MQ=255)
                                                                             aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                       >  1:3147534/1‑100 (MQ=255)
                                                                             aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                       <  2:388907/100‑1 (MQ=255)
                                                                              ataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                       <  2:402153/100‑1 (MQ=255)
                                                                              ataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAcc                      <  1:2563348/100‑1 (MQ=255)
                                                                                 gaaaggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATc                   >  1:2914882/7‑100 (MQ=37)
                                                                                     ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa               >  2:855058/3‑100 (MQ=255)
                                                                                     ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa               >  1:2174361/3‑100 (MQ=255)
                                                                                      gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAACATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  2:947530/2‑100 (MQ=38)
                                                                                      gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  1:3907845/2‑100 (MQ=255)
                                                                                       aaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg             >  2:1535198/1‑100 (MQ=255)
                                                                                        aCGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt            >  1:1958206/1‑100 (MQ=255)
                                                                                        aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt            >  1:3351281/1‑100 (MQ=255)
                                                                                        aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt            >  2:3653209/1‑100 (MQ=255)
                                                                                        aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAAACCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt            >  2:2155799/1‑100 (MQ=39)
                                                                                        aCGAGATACTATTATAAAAGTATCTCGTTCCTTTAATTCAATCCTCACACATGAAAAAATATGAGCTCACCAAAAGTTAAATTCTTAACAATCTTCACg             >  2:3222828/1‑99 (MQ=38)
                                                                                         cGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt           >  1:2120580/1‑100 (MQ=255)
                                                                                         cGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt           >  2:670339/1‑100 (MQ=255)
                                                                                         cGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt           >  1:38174/1‑100 (MQ=255)
                                                                                          gagaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGttt          >  1:1004568/1‑100 (MQ=255)
                                                                                          gagaTACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGttt          >  1:880104/1‑100 (MQ=255)
                                                                                           agaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtttt         >  1:3390239/1‑100 (MQ=255)
                                                                                           agaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtttt         >  1:1821169/1‑100 (MQ=255)
                                                                                            gaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg        <  1:525986/100‑1 (MQ=255)
                                                                                                cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtata    >  2:1035441/1‑100 (MQ=255)
                                                                                                cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtata    >  2:1093798/1‑100 (MQ=255)
                                                                                                cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtata    >  2:3099393/1‑100 (MQ=255)
                                                                                                  tttATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatt  >  1:711149/1‑100 (MQ=255)
                                                                                                  |                                                                                                    
CCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATT  >  NC_000913/2993791‑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: