Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,993,889 T→A 16.8% 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,8890TA16.8% 521.7 / 90.7 214intergenic (+33/‑50)ygeI/pbluncharacterized protein/pseudogene, peptidoglycan‑binding enzyme family
Reads supporting (aligned to +/- strand):  ref base T (96/82);  new base A (27/9);  total (123/91)
Fisher's exact test for biased strand distribution p-value = 2.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

Only 100 of 224 total aligned reads displayed.

GCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATAT  >  NC_000913/2993790‑2993987
                                                                                                   |                                                                                                  
gcCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACtt                                                                                                    >  2:3517647/1‑100 (MQ=255)
 ccTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACttt                                                                                                   >  2:4371933/1‑100 (MQ=255)
  cTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCATATAAAGCATAACGAGATACtttt                                                                                                  <  1:1156088/100‑1 (MQ=255)
   ttttCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTta                                                                                                 <  2:3888541/100‑1 (MQ=255)
    tttCAAAATGATATGTATTCAGAAGATAAGCCTATCTCACCCTACTATAAAATTGAGTTGTAGACGTTTGACCAAATAAATCATAACTGAAAACTTTtat                                                                                                >  1:4075892/1‑100 (MQ=255)
      tCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataa                                                                                              >  2:4205563/1‑100 (MQ=255)
       ggcgatcatatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                             <  2:2806538/93‑1 (MQ=255)
       cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                             >  2:2367800/1‑100 (MQ=255)
       cAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTataat                                                                                             <  1:3544864/100‑1 (MQ=255)
         acaTGATATGTATTCAGAAGAGAATCCTATCTCATCTTACTATAAAATTGAGTAGTAGACTATTATCCAAACAAATCATAACGAGATACTTTTataat                                                                                             <  1:3569427/96‑1 (MQ=255)
           aTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa                                                                                         <  2:3234627/100‑1 (MQ=255)
           aTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa                                                                                         <  1:4355272/100‑1 (MQ=255)
           aTGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTa                                                                                         <  2:597284/100‑1 (MQ=255)
              atatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtct                                                                                      <  2:700708/100‑1 (MQ=255)
               tatGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTAtctc                                                                                     >  1:512527/1‑100 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGt                                                                                   <  2:2325373/100‑1 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                   <  2:1193672/100‑1 (MQ=255)
                 tGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGt                                                                                   <  2:689285/100‑1 (MQ=255)
                    aTTCAGAAGATAATCCTATCTCACCTTACTATAAAATTCAGTTGTAGACGATTGACCAAATAAATCACAACGAGATACTTTTATAATAGTATCTCGTTcc                                                                                <  1:2737740/100‑1 (MQ=255)
                        agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttc                                                                            >  1:1327171/1‑100 (MQ=255)
                        agaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAAAAGTATCTCGttatgatt                                                                            <  2:1783538/100‑7 (MQ=37)
                         gaagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttca                                                                           >  2:173636/1‑100 (MQ=255)
                          aagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcat                                                                          <  2:474886/100‑1 (MQ=255)
                          aagaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcat                                                                          <  1:2162758/100‑1 (MQ=255)
                           agaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcatt                                                                         <  2:4066228/100‑1 (MQ=255)
                           agaTAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTttcat                                                                          >  2:3610392/1‑99 (MQ=255)
                            gaaaaTCCTATATCACCATACTAAAAAATAGAGTGGAAGACGATTGACCAAATAAATAATAACGAGATCCTTTTATAATAGTATATCATTCATtacaatc                                                                        >  1:252372/4‑100 (MQ=255)
                              tAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCaa                                                                      >  2:4311542/1‑100 (MQ=255)
                              tAATCATATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACAAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCaa                                                                      >  1:411668/1‑100 (MQ=255)
                                aTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATc                                                                    >  2:3480341/1‑100 (MQ=255)
                                aTCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATc                                                                    >  2:1598815/1‑100 (MQ=255)
                                 tCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATcc                                                                   <  1:1842235/100‑1 (MQ=255)
                                     aTCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcac                                                               <  1:2817155/100‑1 (MQ=255)
                                     aTCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcac                                                               <  1:3832699/100‑1 (MQ=255)
                                      tctcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcaca                                                              >  2:3773767/1‑100 (MQ=255)
                                       ctcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcacac                                                             >  2:440745/1‑100 (MQ=255)
                                       ctcACCTTACTACAAAATTGTCTGATACCCCATCAACCAAATAAACCATAACGAGATATTTTTATAATACTATCTCGTTCCTTTCATTCAATCCTc                                                                 <  1:3271464/96‑1 (MQ=255)
                                        tcACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTcacaca                                                            >  2:2449749/1‑100 (MQ=255)
                                         cACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACAt                                                           >  1:1879684/1‑100 (MQ=255)
                                         cACCTTACTATAAAATTGAGTGGTAGACGATGGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCACATTCCTTTCAATCAATCCtatcactt                                                           >  1:700317/1‑93 (MQ=255)
                                          aCCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACAtn                                                          <  2:2522562/100‑2 (MQ=255)
                                          aCCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTATTATAATAGTATCTCGTTCCTTTAATTCAATCCTcaca                                                              >  1:3529081/1‑96 (MQ=255)
                                           ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGa                                                         >  2:1055402/1‑100 (MQ=255)
                                           ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGa                                                         >  1:1184238/1‑100 (MQ=255)
                                           ccTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGa                                                         >  2:3866833/1‑100 (MQ=255)
                                                cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaaaaa                                                    >  2:1510439/1‑100 (MQ=255)
                                                cTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGaaaaaa                                                    >  2:3035035/1‑100 (MQ=255)
                                                 tataAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAACAAat                                                   >  1:2834882/1‑100 (MQ=255)
                                                  aaaaaaTTGAGTGATAGACGATTGACCAAATGAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAata                                                  >  2:2967944/3‑100 (MQ=255)
                                                    aaaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATg                                                <  2:3665692/100‑1 (MQ=255)
                                                    aaaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATg                                                <  2:382980/100‑1 (MQ=255)
                                                     aaaTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgc                                               >  2:1947682/1‑100 (MQ=255)
                                                       aTTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATgcgc                                             <  2:1021838/100‑1 (MQ=255)
                                                         tGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTc                                           <  2:4036466/100‑1 (MQ=255)
                                                              ggTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaa                                      >  1:2596302/1‑100 (MQ=255)
                                                              ggTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCa                                       >  1:1069340/1‑99 (MQ=255)
                                                               gTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaaa                                     >  2:428780/1‑100 (MQ=255)
                                                                naGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCaaaa                                    >  2:3508875/2‑100 (MQ=255)
                                                                 aGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                   >  1:1789248/1‑100 (MQ=255)
                                                                  gACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAg                                   >  1:826267/1‑99 (MQ=255)
                                                                    cGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTa                                <  1:3480323/100‑1 (MQ=255)
                                                                      aTTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTaaa                              <  2:1260349/100‑1 (MQ=255)
                                                                       ttGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAAt                             <  1:3216814/100‑1 (MQ=255)
                                                                         gACCAAATAAATCATAACGAGATACTTTTATCATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                           >  1:213440/1‑100 (MQ=255)
                                                                         gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTATATTc                           >  1:143440/1‑100 (MQ=255)
                                                                         gACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTc                           >  2:3982745/1‑100 (MQ=255)
                                                                          aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCt                          >  2:885909/1‑100 (MQ=255)
                                                                          aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCt                          >  1:3360028/1‑100 (MQ=255)
                                                                          aCCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCt                          >  1:1190826/1‑100 (MQ=255)
                                                                           ccAAATAAATCATAACGAGACACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCtt                         >  1:2687766/1‑100 (MQ=255)
                                                                            cAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTa                        >  2:1035230/1‑100 (MQ=255)
                                                                              aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                      <  2:1256070/100‑1 (MQ=255)
                                                                              aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                      >  2:4342896/1‑100 (MQ=255)
                                                                              aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                      <  2:2152718/100‑1 (MQ=255)
                                                                              aataaatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAAc                      <  2:2734626/100‑1 (MQ=255)
                                                                                  gaaaggaaCGAGATACTATTATAAAAGTATCTAGTTCCTTTCATTCAATCTTCACACATGAAAAAATATGCGGTCACCAAAAGTTAAATTCTTAACCATc                  >  1:360401/7‑100 (MQ=37)
                                                                                  aatCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATc                  >  1:2940648/1‑100 (MQ=255)
                                                                                   atCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCt                 >  2:2032534/1‑100 (MQ=255)
                                                                                      ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAATATATGCGCTCACCAAAAGTTAAATACTTAACCATCTTCa              >  1:197531/3‑100 (MQ=38)
                                                                                      ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa              >  2:3398226/3‑100 (MQ=255)
                                                                                      ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa              >  1:2305249/3‑100 (MQ=255)
                                                                                      ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTAACACAGGAAAAACTATGCGCTCACCAAACGTAAAATTCTTAACCATCTTCa              >  1:3011556/3‑100 (MQ=37)
                                                                                      ggaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATCCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa              >  1:63986/3‑100 (MQ=38)
                                                                                      aTAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa              >  1:4158079/1‑100 (MQ=255)
                                                                                      aTAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCa              >  2:4137548/1‑100 (MQ=255)
                                                                                       taccttcaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTGCAc             <  1:3725618/93‑1 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTc               >  1:1028116/2‑98 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc             >  2:1121896/2‑100 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc             >  1:1296754/2‑100 (MQ=255)
                                                                                       gaaCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATACAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCAc             >  1:3254342/2‑100 (MQ=255)
                                                                                         aCGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt           >  1:915141/1‑100 (MQ=255)
                                                                                         aCGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt           <  2:1088326/100‑1 (MQ=255)
                                                                                         aCGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGt           >  2:3312012/1‑100 (MQ=255)
                                                                                          cGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTACCCATCTTCACGtt          >  1:3961015/1‑100 (MQ=39)
                                                                                          cGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt          >  1:2224904/1‑100 (MQ=255)
                                                                                          cGAGATACTATTATAAAAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGtt          >  2:3448166/1‑100 (MQ=255)
                                                                                             gaTACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTg       <  2:3928992/100‑1 (MQ=255)
                                                                                                aCTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtat    <  2:2222791/100‑1 (MQ=255)
                                                                                                 cTTTTATAATAGTATCTCGATTCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTATTAACCATCTGCACGATTTGta     >  1:1551536/1‑98 (MQ=255)
                                                                                                 cTTTTATAATAGTAAATCGTTCCTTTATTTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGtata   <  1:4205606/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GCCTTTTCAAAATGATATGTATTCAGAAGATAATCCTATCTCACCTTACTATAAAATTGAGTGGTAGACGATTGACCAAATAAATCATAACGAGATACTTTTATAATAGTATCTCGTTCCTTTCATTCAATCCTCACACATGAAAAAATATGCGCTCACCAAAAGTTAAATTCTTAACCATCTTCACGTTTTGTATAT  >  NC_000913/2993790‑2993987

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: