Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,993,889 T→A 15.1% 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.1% 459.2 / 67.1 178intergenic (+33/‑50)ygeI/pbluncharacterized protein/pseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base T (82/69);  new base A (12/15);  total (94/84)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01

Only 100 of 187 total aligned reads displayed.

CCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATATT  >  NC_000913/2993791‑2993988
                                                                                                  |                                                                                                   
ccTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGCCCAAATAAATCATAACGAGATACtt                                                                                                     >  2:1114884/1‑99 (MQ=255)
   tttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtat                                                                                                 <  2:2448371/100‑1 (MQ=255)
    ttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTtata                                                                                                <  1:3030625/100‑1 (MQ=255)
      cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                              >  1:185846/1‑100 (MQ=255)
        aaaTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                                                                                            <  1:3425407/100‑1 (MQ=255)
        aaaTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAg                                                                                            >  1:2199826/1‑100 (MQ=255)
            gATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTATACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTAtc                                                                                        <  2:261066/100‑1 (MQ=255)
             atatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtct                                                                                       <  1:222550/100‑1 (MQ=255)
             atatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtct                                                                                       <  2:647030/100‑1 (MQ=255)
              tatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtctc                                                                                      <  1:3590450/100‑1 (MQ=255)
                tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATACATCATAACGAGATACTTTTATAATAGTAtctc                                                                                      >  2:1010062/1‑98 (MQ=255)
                tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                    >  1:3095487/1‑100 (MQ=255)
                tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                    <  2:2753661/100‑1 (MQ=255)
                  tATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTc                                                                                  <  2:1336388/100‑1 (MQ=255)
                  tATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGtta                                                                                  <  1:263605/100‑2 (MQ=255)
                   aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTcc                                                                                 >  1:335553/1‑100 (MQ=255)
                   aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTcc                                                                                 <  1:1658364/100‑1 (MQ=255)
                   aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttat                                                                                 <  1:3730516/100‑3 (MQ=255)
                   aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttat                                                                                 <  1:1993760/100‑3 (MQ=255)
                   aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttat                                                                                 <  1:3527696/100‑3 (MQ=255)
                    ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatg                                                                                <  1:2555110/100‑4 (MQ=38)
                    ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatg                                                                                <  2:2165965/100‑4 (MQ=38)
                    ttCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatg                                                                                <  1:469074/100‑4 (MQ=38)
                     gcAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGTCGATTGACCGAAAAAATCATAACGAGATACTATTATAAAAGTGTCTCGtta                                                                                  <  2:2064488/96‑2 (MQ=37)
                      cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                              >  2:1093556/1‑100 (MQ=255)
                      cAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCttt                                                                              >  2:1271531/1‑100 (MQ=255)
                       agaagaTAATCCTATCTCGCCTTACGATAAAATTGAGTGGTAGACGGTTGACCAAATAAATCATACCGAGATATTTTTATAATAGTATCTGGTTCCtct                                                                              >  2:2178559/1‑97 (MQ=255)
                       agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                             <  2:1433615/100‑1 (MQ=255)
                        gaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttca                                                                            <  1:2118380/100‑1 (MQ=255)
                          agaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                          >  1:3211297/1‑100 (MQ=255)
                          agaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                          <  2:2918106/100‑1 (MQ=255)
                          agaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                          <  1:1492598/100‑1 (MQ=255)
                               aTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATc                                                                     <  2:3185712/100‑1 (MQ=255)
                               aTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATc                                                                     >  1:2345034/1‑100 (MQ=255)
                               aTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATc                                                                     <  2:670878/100‑1 (MQ=255)
                                tCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                    <  2:2661573/100‑1 (MQ=255)
                                tCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                    <  2:2463625/100‑1 (MQ=255)
                                  cTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTc                                                                  <  2:2260466/100‑1 (MQ=255)
                                  cTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTc                                                                  <  2:2132066/100‑1 (MQ=255)
                                   tATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTca                                                                 <  1:2450993/100‑1 (MQ=255)
                                     tctcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcaca                                                               >  2:3647216/1‑100 (MQ=255)
                                      ctcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcacac                                                              <  2:233909/100‑1 (MQ=255)
                                          ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGa                                                          >  1:1799397/1‑100 (MQ=255)
                                           cTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaa                                                         <  1:2538269/100‑1 (MQ=255)
                                            ttACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaa                                                        <  1:2087410/100‑1 (MQ=255)
                                                tataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATGCAATCCTCACACATGAAAAAat                                                    <  2:2687438/100‑1 (MQ=255)
                                                 ataAAATTGAGTGGTAGACGCTTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAata                                                   >  2:2329395/1‑100 (MQ=255)
                                                 ataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAata                                                   >  1:42258/1‑100 (MQ=255)
                                                     aaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgcg                                               <  2:2592555/100‑1 (MQ=255)
                                                      aTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgcgc                                              <  2:42258/100‑1 (MQ=255)
                                                         gAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCa                                           >  1:1475430/1‑100 (MQ=255)
                                                             ggTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaa                                       >  1:794908/1‑100 (MQ=255)
                                                                aGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                    <  1:3510966/100‑1 (MQ=255)
                                                                aGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATAAACTCACAaaaa                                     >  2:3458784/1‑99 (MQ=255)
                                                                 gACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAga                                   >  2:400579/1‑99 (MQ=255)
                                                                  aCGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGtt                                  >  1:124563/1‑100 (MQ=255)
                                                                  aCGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCCCACACATGAAAAAATATGCGCTCACCAAAAGt                                   <  1:2324028/99‑1 (MQ=255)
                                                                   cGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGATCCTTTAAGTAAATCCTAAAACATGAAAAAATATGCGCTAACAaaaa                                     >  2:3601373/1‑96 (MQ=255)
                                                                     aTTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTaaa                               >  1:3739099/1‑100 (MQ=255)
                                                                        naCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                            >  1:1166751/2‑100 (MQ=255)
                                                                        gcccAAATAAATCATAACGAGATACTTTTATAATACTATCTCGTTCCTTTCATTCAATCCTCACACACGAAAAAATATGCTCTCAACAAAAGATAAATTc                            >  2:3637956/3‑100 (MQ=255)
                                                                        gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                            >  2:3458947/1‑100 (MQ=255)
                                                                        gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                            >  1:3236374/1‑100 (MQ=255)
                                                                        gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                            >  1:2983733/1‑100 (MQ=255)
                                                                        gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                            >  1:2850990/1‑100 (MQ=255)
                                                                        gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                            <  1:1191622/100‑1 (MQ=255)
                                                                         aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCt                           >  1:581509/1‑100 (MQ=255)
                                                                         aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCt                           >  2:2241093/1‑100 (MQ=255)
                                                                         aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCt                           <  1:257957/100‑1 (MQ=255)
                                                                          ccAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCtt                          >  1:744121/1‑100 (MQ=255)
                                                                          ccAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCtt                          >  1:641211/1‑100 (MQ=255)
                                                                           cAACTAAATCATCACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTGAATCCTCAAACATGAAAAAATATGCAGTCACCAAAAATTAAATTCTTa                         >  2:1583222/1‑100 (MQ=255)
                                                                           cAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTa                         >  2:181992/1‑100 (MQ=255)
                                                                           cAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCATCACACATGAAAAAATATGCGTTCATCAAAAGTTAAATTGTTa                         >  2:881286/1‑100 (MQ=255)
                                                                            aaataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTaa                        <  2:3316570/100‑1 (MQ=255)
                                                                              ataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAcc                      <  1:1129146/100‑1 (MQ=255)
                                                                               taaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCa                     <  1:3647243/100‑1 (MQ=255)
                                                                               taaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCa                     <  1:727478/100‑1 (MQ=255)
                                                                                 aatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATc                   >  2:519282/1‑100 (MQ=255)
                                                                                 aatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATc                   >  1:1499312/1‑100 (MQ=255)
                                                                                    cATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTc                >  1:76759/1‑100 (MQ=255)
                                                                                    cATAACGAGATACTTCTAGTATAGTATCTCGTTCCATACATTCACTCCTCCCAGATGCAAAAGTTTGAGCTAATCAAAAATCAAATTCTTAAc                       >  2:3491554/1‑93 (MQ=255)
                                                                                     aTAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa               <  1:438311/100‑1 (MQ=255)
                                                                                      tAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              <  2:3720418/100‑1 (MQ=255)
                                                                                      gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  2:2280725/2‑100 (MQ=255)
                                                                                      gaaCGAGATACTATTATAAAAGCATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc              >  1:1457755/2‑100 (MQ=39)
                                                                                       aaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg             >  1:2723368/1‑100 (MQ=255)
                                                                                       aaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACg             >  2:1042344/1‑100 (MQ=255)
                                                                                         ngagaTACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAACAAATATGCGCTGACCCAAAGTTAAATGTTTAACCATGTTCACttt           >  2:3789055/2‑100 (MQ=37)
                                                                                         cGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt           >  1:922010/1‑100 (MQ=255)
                                                                                            gaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg        <  1:36081/100‑1 (MQ=255)
                                                                                            gaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg        <  2:2018577/100‑1 (MQ=255)
                                                                                            gaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg        <  2:1186693/100‑1 (MQ=255)
                                                                                             aTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg        <  2:618799/99‑1 (MQ=255)
                                                                                              tACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGta      <  2:2389340/100‑1 (MQ=255)
                                                                                              tACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGta      <  2:794899/100‑1 (MQ=255)
                                                                                              tACTTGTATAATAGTATCTCGTTCCTTTCATTGAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGt       <  2:3019372/99‑1 (MQ=255)
                                                                                                cTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtata    >  1:2191603/1‑100 (MQ=255)
                                                                                                  tttATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCGCGTTTTGta      <  2:3122378/96‑1 (MQ=255)
                                                                                                  tttATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTAtatt  <  1:1383480/100‑1 (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: