Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,993,889 T→A 19.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,8890TA19.2% 377.4 / 86.1 169intergenic (+33/‑50)ygeI/pbluncharacterized protein/pseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base T (72/64);  new base A (15/17);  total (88/81)
Fisher's exact test for biased strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.91e-01

Only 100 of 176 total aligned reads displayed.

GCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATT  >  NC_000913/2993790‑2993988
                                                                                                   |                                                                                                   
gcCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtt                                                                                                     <  2:2865005/100‑1 (MQ=255)
 ccTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACttt                                                                                                    >  1:3505223/1‑100 (MQ=255)
  cTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtttt                                                                                                   >  2:1628792/1‑100 (MQ=255)
   ttttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTta                                                                                                  >  2:3879361/1‑100 (MQ=255)
   ttttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCAGAACGAGATACTTTta                                                                                                  >  2:4200455/1‑100 (MQ=255)
    tttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtat                                                                                                 <  2:1289966/100‑1 (MQ=255)
     ttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtan                                                                                                 <  1:3089110/99‑2 (MQ=255)
      tCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataa                                                                                               >  2:1666908/1‑100 (MQ=255)
       cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                              >  1:3767263/1‑100 (MQ=255)
       cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                              >  1:329358/1‑100 (MQ=255)
       cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                              >  2:3496219/1‑100 (MQ=255)
       cAAAATGATATGTATCCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                              >  1:828315/1‑100 (MQ=255)
        aaaaTGATAATAATTTAGAAGATAATCGTATCTCACCTTACTATAAAATTGAGTGGTAGACTATTGACCACATAAATCATAACGAGATAATTTTaaaat                                                                                              <  1:796107/99‑1 (MQ=255)
           aTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGAGTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa                                                                                          <  1:1750091/100‑1 (MQ=255)
             gATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtc                                                                                        >  2:2294714/1‑100 (MQ=255)
              atatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATTAATCATAACGAGATACTTTTATAATAGTAtct                                                                                       <  2:4052950/100‑1 (MQ=255)
              atatGTATTCAGAAGATAATCCTATCTCACATTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtct                                                                                       <  1:911814/100‑1 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                    <  2:2657238/100‑1 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                    >  2:3105059/1‑100 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                    >  2:4112613/1‑100 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                    <  2:377994/100‑1 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                    <  1:2525965/100‑1 (MQ=255)
                  gTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGtt                                                                                   >  2:3708817/1‑100 (MQ=255)
                  gTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGtt                                                                                   <  2:4209037/100‑1 (MQ=255)
                   tATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTc                                                                                  >  1:2941033/1‑100 (MQ=255)
                   tATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGtta                                                                                  <  2:4243045/100‑2 (MQ=255)
                   tATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGtta                                                                                  <  2:2100885/100‑2 (MQ=255)
                   tATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGtta                                                                                  <  2:3154812/100‑2 (MQ=255)
                    aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttat                                                                                 <  2:3654058/100‑3 (MQ=255)
                      tCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCtt                                                                               <  2:2549676/100‑1 (MQ=255)
                       cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                              >  2:2548481/1‑100 (MQ=255)
                       cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatga                                                                               <  1:1685693/99‑5 (MQ=38)
                       cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGAGTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                              <  1:3263168/100‑1 (MQ=255)
                       cAGAAGATAATCCTATCGCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                              <  1:3388133/100‑1 (MQ=255)
                        agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                             >  2:1894493/1‑100 (MQ=255)
                        agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                             >  2:3223666/1‑100 (MQ=255)
                        agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                             <  1:1092651/100‑1 (MQ=255)
                        agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                             >  1:1677377/1‑100 (MQ=255)
                        agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatgatt                                                                             <  1:2467720/100‑7 (MQ=37)
                        agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatgatt                                                                             <  1:35628/100‑7 (MQ=37)
                           agaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                          >  2:3473355/1‑100 (MQ=255)
                                 tCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                    <  1:4157345/100‑1 (MQ=255)
                                 tCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                    >  2:1574100/1‑100 (MQ=255)
                                  ccTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCt                                                                   >  2:2348746/1‑100 (MQ=255)
                                    tATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTca                                                                 >  1:129661/1‑100 (MQ=255)
                                    tATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTca                                                                 >  1:3605790/1‑100 (MQ=255)
                                       ctcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcaca                                                               <  1:1379730/99‑1 (MQ=255)
                                        tcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcacaca                                                             >  1:1461097/1‑100 (MQ=255)
                                          aCCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATg                                                           <  2:2941045/100‑1 (MQ=255)
                                           ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGa                                                          <  1:2461092/100‑1 (MQ=255)
                                            cTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTAGTTCCTTTCATTCAATCCTCACACATGaa                                                         >  1:2262325/1‑100 (MQ=255)
                                              tACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaaa                                                       >  2:2040434/1‑100 (MQ=255)
                                                cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaaaaa                                                     <  1:1653805/100‑1 (MQ=255)
                                                  ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAata                                                   >  1:1946977/1‑100 (MQ=255)
                                                  ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAata                                                   <  1:897945/100‑1 (MQ=255)
                                                  ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAata                                                   <  1:4073424/100‑1 (MQ=255)
                                                  ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACAAGATACTTTTATAATAGTATCTAGTTCCTTTCCTTAAATCATCACACATGTAAAAat                                                    >  1:4267818/1‑99 (MQ=255)
                                                         tGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCtn                                            <  2:813323/100‑2 (MQ=255)
                                                          gAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAACAAATATGCGCTc                                            >  2:3564506/1‑99 (MQ=255)
                                                            gTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCAcc                                         >  2:101010/1‑100 (MQ=255)
                                                             ttgTAGACGATTTACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCAct                                         <  1:615919/97‑2 (MQ=255)
                                                             tGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCa                                        >  2:3442563/1‑100 (MQ=255)
                                                              ggTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaa                                       <  2:1753631/100‑1 (MQ=255)
                                                               gTAGACGATTGACCAAATAAATGATAACGAGATACTTTTGTAATAGTATGTCGTTCCTTTTATTCAATCCTCACACATGAAGAAATATgctatca                                           >  1:3538770/1‑90 (MQ=255)
                                                               gTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaaa                                      >  2:2219373/1‑100 (MQ=255)
                                                                tAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTAGTTCATTTCATTTAATTCTAACACATGAAAAAATATTGGCTCACaaaa                                      >  1:2253162/1‑99 (MQ=255)
                                                                 aGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                    <  2:1821397/100‑1 (MQ=255)
                                                                 aGACGATTGACCAAATAAATCATAACGAGATACTCTTATAATAGTATATCGTTCCTTTCATTCAATCCTCACACACGAAAAAATATGCGCACACCAAAAg                                    >  2:2990123/1‑100 (MQ=255)
                                                                   aCGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGtt                                  <  2:2574749/100‑1 (MQ=255)
                                                                      aTTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTaaa                               >  2:4211999/1‑100 (MQ=255)
                                                                       ttGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAAt                              <  2:2475572/100‑1 (MQ=255)
                                                                           ccAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCtt                          >  1:3821763/1‑100 (MQ=255)
                                                                           ccAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCtt                          >  1:1532921/1‑100 (MQ=255)
                                                                            cAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTa                         <  2:572851/100‑1 (MQ=255)
                                                                              aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                       <  2:2853027/100‑1 (MQ=255)
                                                                                  aatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATc                   >  1:2481026/1‑100 (MQ=255)
                                                                                     cATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTc                <  2:1294382/100‑1 (MQ=255)
                                                                                      ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa               >  2:1998606/3‑100 (MQ=255)
                                                                                      ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa               >  1:220378/3‑100 (MQ=255)
                                                                                      ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa               >  1:1329564/3‑100 (MQ=255)
                                                                                      aTAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCAAAAAAAAAAAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa               <  1:639932/100‑1 (MQ=255)
                                                                                       tAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              <  1:1406714/100‑1 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  2:4143947/2‑100 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  2:1519733/2‑100 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  2:1934807/2‑100 (MQ=255)
                                                                                        aaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg             >  2:1856651/1‑100 (MQ=255)
                                                                                         aCGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt            >  1:2661732/1‑100 (MQ=255)
                                                                                         aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg             >  2:2537911/1‑99 (MQ=255)
                                                                                         aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt            >  2:3310819/1‑100 (MQ=255)
                                                                                         aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt            >  2:1253467/1‑100 (MQ=255)
                                                                                            agaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtttt         >  1:2808755/1‑100 (MQ=255)
                                                                                             gaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCTCCAAAAGTTAAATTCTTAACCAACTTCACGTTTTg        <  2:1110537/100‑1 (MQ=255)
                                                                                             gaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg        <  1:200245/100‑1 (MQ=255)
                                                                                               tACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGta      >  1:4003182/1‑100 (MQ=255)
                                                                                               tACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGta      <  2:2819316/100‑1 (MQ=255)
                                                                                                aCTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtat     <  1:2296804/100‑1 (MQ=255)
                                                                                                aCTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtat     <  2:1890273/100‑1 (MQ=255)
                                                                                                aCTTTTATAATAGCAGCTCGTTCCTTTACTTCAATCCTCACACATGAAAAAATATGCGCTCACCAAATGTTAAATTCTTAACCATCTTCACGt            <  1:2946957/93‑1 (MQ=255)
                                                                                                 cTTTTATAATAGTATCTCGCTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtata    >  1:2101849/1‑100 (MQ=255)
                                                                                                  ttttATAATAGTATCTCGTTCCTTTCATTCAATCCTGACACATGAAAAAATATGGGCTCAATCAAAATTAAATTCTTAACCATCTTCATTTTTTGTatat   >  1:1044093/1‑100 (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: