Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 531,453 T→C intergenic (‑8/‑222) ybbS ← / → ybbT putative transcriptional regulator LYSR‑type/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913531,4530TC99.2% 776.5 / ‑5.1 240intergenic (‑8/‑222)ybbS/ybbTputative transcriptional regulator LYSR‑type/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (131/107);  total (133/107)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 241 total aligned reads displayed.

TTACATAATCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAATTTGTGGACT  >  NC_000913/531363‑531542
                                                                                          |                                                                                         
ttACATAATCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCAt                                                                                         >  2:4188320/1‑93 (MQ=255)
 tACATAATCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATc                                                                                        >  1:342459/1‑93 (MQ=255)
 tACATAATCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATc                                                                                        >  1:352690/1‑93 (MQ=255)
  aCATAATCGTTCTGCCGCTTTTGAAAAGCTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATcc                                                                                       <  1:4613078/93‑1 (MQ=255)
  aCATAATCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGCACATCTTCAGGCATcc                                                                                       >  2:2137400/1‑93 (MQ=255)
  aCATAATCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATcc                                                                                       <  2:1473491/93‑1 (MQ=255)
    aTAATCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATcccc                                                                                     <  1:3105370/93‑1 (MQ=255)
    aTAATCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATcccc                                                                                     >  2:3102965/1‑93 (MQ=255)
    aTAATCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATcccc                                                                                     >  1:4177674/1‑93 (MQ=255)
     tAATCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCt                                                                                    >  1:1771737/1‑93 (MQ=255)
     tAATCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCt                                                                                    >  1:2209207/1‑93 (MQ=255)
     tAATCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCt                                                                                    >  1:3932186/1‑93 (MQ=255)
      aaTCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCtt                                                                                   <  2:3009944/93‑1 (MQ=255)
        tCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCtttt                                                                                 <  1:539415/93‑1 (MQ=255)
        tCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCtttt                                                                                 <  1:2711172/93‑1 (MQ=255)
          gTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTaa                                                                               >  1:1433301/1‑93 (MQ=255)
            tCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAAt                                                                             <  1:2874044/93‑1 (MQ=255)
             cTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATc                                                                            <  2:597584/93‑1 (MQ=255)
              tGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATcc                                                                           <  2:385864/93‑1 (MQ=255)
              tGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATcc                                                                           <  2:1509190/93‑1 (MQ=255)
              tGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATcc                                                                           >  2:424754/1‑93 (MQ=255)
               gCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCg                                                                          <  2:103799/93‑1 (MQ=255)
                  gCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCaa                                                                       <  2:4521898/93‑1 (MQ=255)
                  gCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCaa                                                                       >  1:2448568/1‑93 (MQ=255)
                  gCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCaa                                                                       >  2:4085959/1‑93 (MQ=255)
                  gCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCaa                                                                       >  1:1741177/1‑93 (MQ=255)
                    ttttGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGt                                                                     <  2:3283029/93‑1 (MQ=255)
                     tttGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGtt                                                                    <  1:4033316/93‑1 (MQ=255)
                      ttGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTg                                                                   >  1:3215909/1‑93 (MQ=255)
                         aaaaaCTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGt                                                                <  1:4518990/93‑1 (MQ=255)
                          aaaaCTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTg                                                               >  1:3494976/1‑93 (MQ=255)
                           aaaCTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGa                                                              >  2:2282844/1‑93 (MQ=255)
                           aaaCTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGa                                                              >  1:1334549/1‑93 (MQ=255)
                            aaCTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGAt                                                             <  2:1637469/93‑1 (MQ=255)
                             aCTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGAtt                                                            >  1:2267695/1‑93 (MQ=255)
                             aCTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGAtt                                                            <  1:1134404/93‑1 (MQ=255)
                              cTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTCAATCCGCAAGTTGCGTGCttt                                                           >  2:3593100/1‑93 (MQ=255)
                              cTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGAttt                                                           >  1:2778890/1‑93 (MQ=255)
                              cTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGAttt                                                           >  2:3631425/1‑93 (MQ=255)
                              cTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGAttt                                                           >  1:747689/1‑93 (MQ=255)
                              cTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGAttt                                                           <  2:4728888/93‑1 (MQ=255)
                                tCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTc                                                         >  2:1843268/1‑93 (MQ=255)
                                 ccTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCt                                                        <  2:1180778/93‑1 (MQ=255)
                                 ccTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCt                                                        <  1:1961627/93‑1 (MQ=255)
                                   tGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCttt                                                      >  2:875749/1‑92 (MQ=255)
                                   tGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTa                                                      <  2:2194995/93‑1 (MQ=255)
                                     tttCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTCAATCCGCAAGTTGCGTGATTTTCTTATc                                                    >  1:1986356/1‑93 (MQ=255)
                                      ttCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATcc                                                   >  2:1614817/1‑93 (MQ=255)
                                         aGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCctct                                                <  2:916678/93‑1 (MQ=255)
                                         aGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCctct                                                <  1:1898845/93‑1 (MQ=255)
                                         aGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCctct                                                <  1:357460/93‑1 (MQ=255)
                                          gCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTg                                               >  2:3703972/1‑93 (MQ=255)
                                          gCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTg                                               >  2:4148715/1‑93 (MQ=255)
                                           cAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGa                                              <  2:1239052/93‑1 (MQ=255)
                                            aaCCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGAt                                             >  2:2771197/1‑93 (MQ=255)
                                            aaCCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGAt                                             >  2:4448077/1‑93 (MQ=255)
                                              ccGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGAttt                                           >  2:2438348/1‑93 (MQ=255)
                                               cGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTa                                          <  2:2278789/93‑1 (MQ=255)
                                                gCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTAt                                         <  1:1038168/93‑1 (MQ=255)
                                                 cAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTGAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATc                                        <  2:3493996/93‑1 (MQ=255)
                                                  aaTGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCa                                       <  1:260494/93‑1 (MQ=255)
                                                    tGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTAATTTTCTTATCCTCTGATTTATCAGt                                     >  1:2112081/1‑93 (MQ=255)
                                                     gAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTa                                    <  1:1192676/93‑1 (MQ=255)
                                                     gAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTa                                    >  2:2122931/1‑93 (MQ=255)
                                                       aaGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTAtt                                  <  2:1108761/93‑1 (MQ=255)
                                                         gTCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTAtttt                                <  2:3439070/93‑1 (MQ=255)
                                                          tctgCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTAttttt                               <  1:89461/90‑1 (MQ=255)
                                                          tCCGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTAttttt                               >  2:3594105/1‑93 (MQ=255)
                                                           ccGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTa                              <  2:526398/93‑1 (MQ=255)
                                                           ccGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTa                              <  2:756257/93‑1 (MQ=255)
                                                            cGCAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTAc                             <  1:3467279/93‑1 (MQ=255)
                                                              cAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACAt                           <  2:3333958/93‑1 (MQ=255)
                                                              cAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACAt                           <  1:37571/93‑1 (MQ=255)
                                                              cAAGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACAt                           <  1:782327/93‑1 (MQ=255)
                                                               aaGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATg                          >  2:709709/1‑93 (MQ=255)
                                                                aGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGa                         >  2:3144400/1‑93 (MQ=255)
                                                                aGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGa                         <  1:2396624/93‑1 (MQ=255)
                                                                aGGTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGa                         >  2:2690986/1‑93 (MQ=255)
                                                                 ggTTTCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGAt                        >  1:3013754/1‑93 (MQ=255)
                                                                    ttCTGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAAc                     <  1:4336893/93‑1 (MQ=255)
                                                                       tGGATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCt                  >  1:3076437/1‑93 (MQ=255)
                                                                        ggATCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTg                 <  2:4015439/93‑1 (MQ=255)
                                                                          aTCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGtt               <  2:2264522/93‑1 (MQ=255)
                                                                          aTCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGtt               <  2:4315659/93‑1 (MQ=255)
                                                                           tCGAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTc              <  2:1584120/93‑1 (MQ=255)
                                                                             gAACATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCaa            <  1:2006317/93‑1 (MQ=255)
                                                                              aaCATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAAt           >  1:2631639/1‑93 (MQ=255)
                                                                              aaCATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAAt           >  1:2956946/1‑93 (MQ=255)
                                                                               aCATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAAtt          >  2:2500133/1‑93 (MQ=255)
                                                                               aCATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAAtt          >  1:1560901/1‑93 (MQ=255)
                                                                               aCATCTTCAGGCATCCCCTTTTAAATCCAGAAGGCGCGTGGTTTTCTTCTGCTCTTATTTCTCAGTATTTGTACATTATCACCCTTATCAAtt          >  1:4306920/1‑93 (MQ=255)
                                                                                cATCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAAttt         >  2:3293019/1‑93 (MQ=255)
                                                                                 aTCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAATTtg        >  2:2923981/1‑93 (MQ=255)
                                                                                  tCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAATTtgt       <  1:3845148/93‑1 (MQ=255)
                                                                                  tCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAATTtgt       >  2:399704/1‑93 (MQ=255)
                                                                                  tCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAATTtgt       <  2:3875637/93‑1 (MQ=255)
                                                                                  tCTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAATTtgt       <  1:659782/93‑1 (MQ=255)
                                                                                   cTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAATTtgtg      <  1:4465850/93‑1 (MQ=255)
                                                                                   cTTCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAATTtgtg      >  1:1007539/1‑93 (MQ=255)
                                                                                     tCAGGCATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAATTTGTGGa    <  2:4172158/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTACATAATCGTTCTGCCGCTTTTGAAAAACTTCCTGTTTCAGCAACCGCAATGAAAGTCCGCAAGGTTTCTGGATCGAACATCTTCAGGTATCCCCTTTTAAATCCGCAAGTTGCGTGATTTTCTTATCCTCTGATTTATCAGTATTTTTACATGATAACCCTGTTCAATTTGTGGACT  >  NC_000913/531363‑531542

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: