Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,998,039 A→G Y124H (TAC→CAC)  fliY ← putative periplasmic binding transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,998,0390AG99.2% 380.2 / NA 122Y124H (TAC→CAC) fliYputative periplasmic binding transport protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (64/57);  total (64/58)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

CCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGA  >  NC_000913/1997949‑1998126
                                                                                          |                                                                                       
cccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCg                                                                                       >  1:3991489/1‑93 (MQ=255)
 ccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCgg                                                                                      >  1:296106/1‑93 (MQ=255)
 ccccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCgg                                                                                      >  2:3356826/1‑93 (MQ=255)
  cccTCTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCggg                                                                                     <  1:1727990/93‑1 (MQ=255)
  cccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCggg                                                                                     >  2:2829047/1‑93 (MQ=255)
  cccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCggg                                                                                     >  1:870728/1‑93 (MQ=255)
   ccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCgggg                                                                                    <  1:4367975/93‑1 (MQ=255)
   ccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCgggg                                                                                    >  2:2487565/1‑93 (MQ=255)
     aCTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGtt                                                                                  >  1:1527881/1‑93 (MQ=255)
     aCTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGtt                                                                                  <  2:1098801/93‑1 (MQ=255)
      cTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTTGTTACCTTTTTTCACCAGCGGCCGAATACCTGAAATGGTGCGCGGGGCTg                                                                                 >  2:1225133/1‑93 (MQ=255)
       ttttttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTga                                                                                >  1:2501142/1‑93 (MQ=255)
         ttttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTgaga                                                                              >  1:1571566/1‑93 (MQ=255)
           ttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGaaa                                                                            <  2:970311/93‑1 (MQ=255)
           ttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGaaa                                                                            >  2:965884/1‑93 (MQ=255)
             gCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATc                                                                          <  2:1883156/93‑1 (MQ=255)
             gCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATc                                                                          <  1:1207083/93‑1 (MQ=255)
                tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTa                                                                       >  2:3112864/1‑93 (MQ=255)
                tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGGTGAGAAATCGTa                                                                       >  1:2261651/1‑93 (MQ=255)
                 ttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTAt                                                                      <  2:2658851/93‑1 (MQ=255)
                    aGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTAtttt                                                                   <  2:1850789/93‑1 (MQ=255)
                    aGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTAtttt                                                                   <  1:1437213/93‑1 (MQ=255)
                       tcatcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTc                                                                >  1:2562213/1‑93 (MQ=255)
                        catcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCt                                                               >  2:4521152/1‑93 (MQ=255)
                        catcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCt                                                               <  1:3849358/93‑1 (MQ=255)
                        catcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCt                                                               >  2:4269485/1‑93 (MQ=255)
                          tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTg                                                             >  2:3151699/1‑93 (MQ=255)
                          tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTg                                                             >  2:2740650/1‑93 (MQ=255)
                          tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTg                                                             >  2:2739028/1‑93 (MQ=255)
                          tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTg                                                             >  2:738029/1‑93 (MQ=255)
                          tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTg                                                             >  1:3026705/1‑93 (MQ=255)
                           cGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTgc                                                            <  1:3165768/93‑1 (MQ=255)
                             gCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTgcgc                                                          >  1:1698456/1‑93 (MQ=255)
                             gCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTgcgc                                                          >  1:4204097/1‑93 (MQ=255)
                               tGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtc                                                        >  1:565557/1‑93 (MQ=255)
                               tGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtc                                                        >  1:2420038/1‑93 (MQ=255)
                                gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtca                                                       <  2:4040562/93‑1 (MQ=255)
                                gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtca                                                       >  2:1125860/1‑93 (MQ=255)
                                gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtca                                                       <  2:2161830/93‑1 (MQ=255)
                                gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtca                                                       <  2:4040561/93‑1 (MQ=255)
                                 ttttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcat                                                      <  2:584254/93‑1 (MQ=255)
                                  tttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcatc                                                     <  1:1577080/93‑1 (MQ=255)
                                   ttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAATTCGTATTTTTTCTTGCGCtcatca                                                    <  2:4518312/93‑1 (MQ=255)
                                   ttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcatca                                                    <  2:3296868/93‑1 (MQ=255)
                                    tAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAg                                                   <  1:3049078/93‑1 (MQ=255)
                                      aTGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaa                                                 <  1:822193/93‑1 (MQ=255)
                                      aTGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaa                                                 >  1:383797/1‑93 (MQ=255)
                                       tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaaa                                                >  1:210065/1‑93 (MQ=255)
                                       tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaaa                                                >  1:2761554/1‑93 (MQ=255)
                                       tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaaa                                                >  2:2580084/1‑93 (MQ=255)
                                       tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGCGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaaa                                                <  1:2507254/93‑1 (MQ=255)
                                       tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATAGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaaa                                                <  2:3795981/93‑1 (MQ=255)
                                        ggTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAAt                                               >  1:2083974/1‑93 (MQ=255)
                                         gTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATg                                              >  2:2364074/1‑93 (MQ=255)
                                           gCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGt                                            <  1:2056424/93‑1 (MQ=255)
                                            ccTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTg                                           >  2:459237/1‑93 (MQ=255)
                                             cTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGa                                          <  2:3042617/93‑1 (MQ=255)
                                             cTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGa                                          >  2:3238739/1‑93 (MQ=255)
                                              ttCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAc                                         >  1:3392628/1‑93 (MQ=255)
                                               tCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAcc                                        <  2:769425/93‑1 (MQ=255)
                                               tCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAcc                                        <  1:1469068/93‑1 (MQ=255)
                                               tCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAcc                                        >  1:562417/1‑93 (MQ=255)
                                                cGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCt                                       >  1:521775/1‑93 (MQ=255)
                                                cGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCt                                       <  1:216285/93‑1 (MQ=255)
                                                 gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTg                                      >  2:1132605/1‑93 (MQ=255)
                                                 gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTg                                      >  1:1747263/1‑93 (MQ=255)
                                                  ttcccTTTTTTCGCCAGCGCCTGAATACCAGAAATGGTGTGCGGTGTTGAGAAATCGTATTTTTTCTGGCGCTCATCAGAAATGGTGACCTGa                                     <  2:2044379/90‑1 (MQ=255)
                                                  ttACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGGAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGa                                     >  1:746308/1‑93 (MQ=255)
                                                  ttACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGa                                     <  1:3381215/93‑1 (MQ=255)
                                                  ttACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGa                                     <  1:2406350/93‑1 (MQ=255)
                                                   tACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGAt                                    <  1:3040802/93‑1 (MQ=255)
                                                   tACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGAt                                    <  2:4674138/93‑1 (MQ=255)
                                                   tACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCCGAt                                    <  2:4063884/93‑1 (MQ=255)
                                                     ccTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTa                                  >  2:3602549/1‑93 (MQ=255)
                                                      cTTTTTTCACCAGCGCGTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTaa                                 <  2:1569766/93‑1 (MQ=255)
                                                       ttttttCCCCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt                                <  1:729721/93‑1 (MQ=255)
                                                       ttttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt                                <  2:2670966/93‑1 (MQ=255)
                                                       ttttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt                                >  1:2232700/1‑93 (MQ=255)
                                                       ttttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGGTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt                                >  2:1933929/1‑93 (MQ=255)
                                                        tttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATc                               <  1:543764/93‑1 (MQ=255)
                                                        tttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATc                               <  1:4180337/93‑1 (MQ=255)
                                                          tttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcac                             >  2:867667/1‑93 (MQ=255)
                                                          tttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcac                             <  1:279379/93‑1 (MQ=255)
                                                              aCCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACcaca                         >  1:1430559/1‑93 (MQ=255)
                                                              aCCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACcaca                         >  1:181417/1‑93 (MQ=255)
                                                                cAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATc                       <  2:4603167/93‑1 (MQ=255)
                                                                  gcgcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCaa                     <  2:3739501/93‑1 (MQ=255)
                                                                       tGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGt                >  2:947450/1‑93 (MQ=255)
                                                                       tGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGt                >  1:4472508/1‑93 (MQ=255)
                                                                          aTACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTa             >  2:3320805/1‑93 (MQ=255)
                                                                           tACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTag            >  1:1599306/1‑93 (MQ=255)
                                                                             ccAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTagag          <  1:2865799/93‑1 (MQ=255)
                                                                             ccAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTagag          <  2:3023767/93‑1 (MQ=255)
                                                                               aGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc        <  2:2821170/93‑1 (MQ=255)
                                                                               aGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc        <  2:3832486/93‑1 (MQ=255)
                                                                                gAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTcc       >  2:2587933/1‑93 (MQ=255)
                                                                                gAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTcc       <  1:3846927/93‑1 (MQ=255)
                                                                                 aaaTGTTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCa      <  1:2940804/93‑1 (MQ=255)
                                                                                     ggTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCGGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGa  >  1:2512474/1‑93 (MQ=255)
                                                                                     ggTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGa  <  1:357791/93‑1 (MQ=255)
                                                                                          |                                                                                       
CCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGA  >  NC_000913/1997949‑1998126

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