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.

CCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATT  >  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)
                                          ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑a                                                          >  1:1499058/1‑100 (MQ=255)
                                          cattACTATAAAACTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑a                                                          >  1:2676468/3‑100 (MQ=255)
                                           cTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑aa                                                         <  1:273693/100‑1 (MQ=255)
                                            ttACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑aaa                                                        >  2:3386422/1‑100 (MQ=255)
                                              aCTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑aaaaa                                                      <  1:2352254/100‑1 (MQ=255)
                                               cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑aaaaaa                                                     <  1:3475004/100‑1 (MQ=255)
                                               cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑aaaaaa                                                     >  1:3315742/1‑100 (MQ=255)
                                                tataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAat                                                    >  2:3940566/1‑100 (MQ=255)
                                                 ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAata                                                   >  1:1783147/1‑100 (MQ=255)
                                                   aaaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATg                                                 >  2:2754962/1‑100 (MQ=255)
                                                   aaaaTTAAGTGGTAGACGATTAACAAAATAAATCATAACGAGATAGTTTTATAATAGTATCTCGTTCTTTTACTTCAATTCTCATAAATG‑AAAAAATATg                                                 >  2:1906480/1‑100 (MQ=255)
                                                     aaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATgcg                                               <  1:3677431/100‑1 (MQ=255)
                                                      aTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATgcgc                                              >  2:2386889/1‑100 (MQ=255)
                                                           gTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCAcc                                         >  1:3002206/1‑100 (MQ=255)
                                                            tGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCa                                        <  1:2445378/100‑1 (MQ=255)
                                                              gTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCaaa                                      <  1:3084206/100‑1 (MQ=255)
                                                               tAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCaaaa                                     >  1:937873/1‑100 (MQ=255)
                                                                aGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAg                                    <  2:3854836/100‑1 (MQ=255)
                                                                  aCGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCaaa                                      <  2:225010/96‑1 (MQ=255)
                                                                   cGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTACTTTCATTAAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTa                                 >  2:1922760/1‑100 (MQ=255)
                                                                    gATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATACTCACAAATG‑AAAAAATTTGTGTTCACAAAAAGTTaa                                >  2:964142/1‑100 (MQ=255)
                                                                      ttGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAAt                              <  2:1783138/100‑1 (MQ=255)
                                                                      ttGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAAt                              <  2:872194/100‑1 (MQ=255)
                                                                        gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTc                            >  1:2499611/1‑100 (MQ=255)
                                                                          ccAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCGCCAAAAGTTAAATTCtt                          >  1:2215760/1‑100 (MQ=255)
                                                                             aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                       >  1:3147534/1‑100 (MQ=255)
                                                                             aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                       <  2:388907/100‑1 (MQ=255)
                                                                              ataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                       <  2:402153/100‑1 (MQ=255)
                                                                              ataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAcc                      <  1:2563348/100‑1 (MQ=255)
                                                                                 gaaaggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATc                   >  1:2914882/7‑100 (MQ=37)
                                                                                     ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa               >  2:855058/3‑100 (MQ=255)
                                                                                     ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa               >  1:2174361/3‑100 (MQ=255)
                                                                                      gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAACATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  2:947530/2‑100 (MQ=38)
                                                                                      gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  1:3907845/2‑100 (MQ=255)
                                                                                       aaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg             >  2:1535198/1‑100 (MQ=255)
                                                                                        aCGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt            >  1:1958206/1‑100 (MQ=255)
                                                                                        aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt            >  1:3351281/1‑100 (MQ=255)
                                                                                        aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt            >  2:3653209/1‑100 (MQ=255)
                                                                                        aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAAACCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt            >  2:2155799/1‑100 (MQ=39)
                                                                                        aCGAGATACTATTATAAAAGTATCTCGTTCCTTTAATTCAATCCTCACACATG‑AAAAAATATGAGCTCACCAAAAGTTAAATTCTTAACAATCTTCACg             >  2:3222828/1‑99 (MQ=38)
                                                                                         cGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt           >  1:2120580/1‑100 (MQ=255)
                                                                                         cGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt           >  2:670339/1‑100 (MQ=255)
                                                                                         cGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt           >  1:38174/1‑100 (MQ=255)
                                                                                          gagaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGttt          >  1:1004568/1‑100 (MQ=255)
                                                                                          gagaTACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGttt          >  1:880104/1‑100 (MQ=255)
                                                                                           agaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtttt         >  1:3390239/1‑100 (MQ=255)
                                                                                           agaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtttt         >  1:1821169/1‑100 (MQ=255)
                                                                                            gaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg        <  1:525986/100‑1 (MQ=255)
                                                                                                cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtata    >  2:1035441/1‑100 (MQ=255)
                                                                                                cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtata    >  2:1093798/1‑100 (MQ=255)
                                                                                                cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtata    >  2:3099393/1‑100 (MQ=255)
                                                                                                  tttATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatt  >  1:711149/1‑100 (MQ=255)
                                                                                                  |                                                                                                    
CCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATG‑AAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATT  >  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: ‑