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.6% 725.5 / ‑5.2 232Y124H (TAC→CAC) fliYputative periplasmic binding transport protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (114/117);  total (115/117)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 237 total aligned reads displayed.

CCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACG  >  NC_000913/1997949‑1998128
                                                                                          |                                                                                         
  cccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCggg                                                                                       <  2:5440667/93‑1 (MQ=255)
  cccACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCggg                                                                                       >  1:3057530/1‑93 (MQ=255)
  cccACTTTTTTGCCTTTCAGATCATCGGCTCTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCggg                                                                                       <  1:949476/93‑1 (MQ=255)
    cACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGt                                                                                     <  2:3600493/93‑1 (MQ=255)
       ttttttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTga                                                                                  >  1:2073162/1‑93 (MQ=255)
         ttttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTgaga                                                                                <  1:3266194/93‑1 (MQ=255)
           ttGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGaaa                                                                              <  1:4810345/93‑1 (MQ=255)
            tGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAAt                                                                             <  2:6621717/93‑1 (MQ=255)
            tGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAAt                                                                             <  1:4008600/93‑1 (MQ=255)
             gCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATc                                                                            <  1:4899535/93‑1 (MQ=255)
             gCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATc                                                                            >  1:4082806/1‑93 (MQ=255)
              ccTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCg                                                                           <  1:834440/93‑1 (MQ=255)
               cTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGt                                                                          >  2:4367174/1‑93 (MQ=255)
                tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTa                                                                         >  2:1032376/1‑93 (MQ=255)
                tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTa                                                                         <  2:4422220/93‑1 (MQ=255)
                tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTa                                                                         >  2:1732634/1‑93 (MQ=255)
                tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTa                                                                         >  1:229806/1‑93 (MQ=255)
                tttCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTGTTTTCACCAGCGCCTGAACACCAGAAATGGTGTGCGGGGTTGAGAAATCGTa                                                                         <  1:6732815/93‑1 (MQ=255)
                  tCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTAtt                                                                       <  2:4978460/93‑1 (MQ=255)
                    aGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTAtttt                                                                     >  1:3843815/1‑93 (MQ=255)
                    aGATCATCGGCTGTTTTAATGGTGACTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTAtttt                                                                     <  2:2507337/93‑1 (MQ=255)
                       tcatcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTc                                                                  >  1:1145953/1‑93 (MQ=255)
                        catcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCt                                                                 <  1:6203921/93‑1 (MQ=255)
                         atcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCtt                                                                <  1:5095921/93‑1 (MQ=255)
                          tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTg                                                               <  2:2136462/93‑1 (MQ=255)
                          tcGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTg                                                               >  2:5570464/1‑93 (MQ=255)
                           cGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTgc                                                              <  1:1810949/93‑1 (MQ=255)
                           cGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTgc                                                              <  1:6170859/93‑1 (MQ=255)
                           cGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTgc                                                              <  1:5804391/93‑1 (MQ=255)
                              cTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCt                                                           <  2:1934172/93‑1 (MQ=255)
                              cTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCt                                                           >  2:845814/1‑93 (MQ=255)
                               tGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtc                                                          >  1:416569/1‑93 (MQ=255)
                               tGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtc                                                          >  1:2185783/1‑93 (MQ=255)
                                gTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtca                                                         >  1:556056/1‑93 (MQ=255)
                                 ttttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcat                                                        >  2:1640964/1‑93 (MQ=255)
                                 ttttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcat                                                        <  2:7355273/93‑1 (MQ=255)
                                  tttAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCtcatc                                                       >  1:1301947/1‑93 (MQ=255)
                                       tGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGaaa                                                  >  2:608758/1‑93 (MQ=255)
                                         gTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATg                                                <  1:5584922/93‑1 (MQ=255)
                                         gTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATg                                                <  2:5792113/93‑1 (MQ=255)
                                         gTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATg                                                <  1:3756782/93‑1 (MQ=255)
                                          tGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATgg                                               <  1:6530021/93‑1 (MQ=255)
                                          tGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATgg                                               <  1:6929333/93‑1 (MQ=255)
                                            ccTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTg                                             >  1:4491115/1‑93 (MQ=255)
                                             cTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGa                                            >  1:174925/1‑93 (MQ=255)
                                              ttCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGAc                                           <  2:4870785/93‑1 (MQ=255)
                                                 gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTg                                        >  1:601074/1‑93 (MQ=255)
                                                 gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTg                                        <  1:5405802/93‑1 (MQ=255)
                                                 gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTg                                        >  1:7216061/1‑93 (MQ=255)
                                                 gTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTg                                        >  2:3098123/1‑93 (MQ=255)
                                                  ttACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGa                                       >  1:1383622/1‑93 (MQ=255)
                                                  ttACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGa                                       >  2:7481789/1‑93 (MQ=255)
                                                  ttACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGa                                       <  2:3136013/93‑1 (MQ=255)
                                                   tACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGAt                                      >  2:7236295/1‑93 (MQ=255)
                                                   tACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGAt                                      >  1:7764763/1‑93 (MQ=255)
                                                    aCCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGAtt                                     <  2:7836315/93‑1 (MQ=255)
                                                     ccTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTa                                    >  2:6057887/1‑93 (MQ=255)
                                                     ccTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTa                                    <  2:3234947/93‑1 (MQ=255)
                                                      cTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTaa                                   >  2:7112162/1‑93 (MQ=255)
                                                       ttttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt                                  >  1:1613507/1‑93 (MQ=255)
                                                       ttttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAAt                                  >  1:3498670/1‑93 (MQ=255)
                                                         ttttCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATca                                >  1:5526818/1‑93 (MQ=255)
                                                         tgttCACCAGCGCCTGAGTACCATAAATGGTGTGCGGGGTTGAGAAATCTTGTTTTTTCTTGTGGTCATCAGAAATGGTGACCTGATTAATca                                <  1:2967462/91‑1 (MQ=255)
                                                            tCACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcacca                             <  2:1668829/93‑1 (MQ=255)
                                                             cACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcaccac                            <  1:2474458/93‑1 (MQ=255)
                                                             cACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcaccac                            <  1:6964244/93‑1 (MQ=255)
                                                             cACCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATcaccac                            >  1:7559582/1‑93 (MQ=255)
                                                              aCCAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACcaca                           >  2:3804126/1‑93 (MQ=255)
                                                                cAGCGCCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATc                         >  2:2234355/1‑93 (MQ=255)
                                                                  gcgcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCaa                       <  1:861569/93‑1 (MQ=255)
                                                                  acgcCTGAATACCAGAAATGGTGAGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCaa                       <  2:5667871/92‑1 (MQ=255)
                                                                   cgcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAAt                      <  1:4383515/93‑1 (MQ=255)
                                                                    gcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATa                     >  1:846810/1‑93 (MQ=255)
                                                                    gcCTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATa                     <  1:5287050/93‑1 (MQ=255)
                                                                     ccTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATAc                    <  1:6246920/93‑1 (MQ=255)
                                                                     ccTGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATAc                    <  1:1945197/93‑1 (MQ=255)
                                                                       tGAATACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGt                  >  1:3955806/1‑93 (MQ=255)
                                                                        gAATACCAGACATGGCGTGCGGGGTTGAGAAATCGCATTTTTTCTTGCGCTCATCTGAAATGGTGACCTGATTAATAACCACATCAATACGtt                 >  2:6270529/1‑93 (MQ=255)
                                                                         aaTACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGttt                <  2:7653227/93‑1 (MQ=255)
                                                                          aTACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTa               <  1:2154369/93‑1 (MQ=255)
                                                                          aTACCAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTa               <  2:5947270/93‑1 (MQ=255)
                                                                              cAGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGt           >  1:7257998/1‑93 (MQ=255)
                                                                               aGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc          >  2:3649127/1‑93 (MQ=255)
                                                                               aGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc          <  2:4680027/93‑1 (MQ=255)
                                                                               aGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc          <  1:2400449/93‑1 (MQ=255)
                                                                               aGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc          <  2:5300492/93‑1 (MQ=255)
                                                                               aGAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTc          >  1:2749801/1‑93 (MQ=255)
                                                                                gAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTcc         <  2:584743/93‑1 (MQ=255)
                                                                                gAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTcc         <  1:3118037/93‑1 (MQ=255)
                                                                                gAAATGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTcc         >  2:4779436/1‑93 (MQ=255)
                                                                                 aaaTGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCa        <  2:4749152/93‑1 (MQ=255)
                                                                                 aaaTGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCa        <  2:1167701/93‑1 (MQ=255)
                                                                                   aTGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGc      >  2:3426955/1‑93 (MQ=255)
                                                                                    tGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCg     <  1:6653655/93‑1 (MQ=255)
                                                                                    tGGTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCg     >  1:738266/1‑93 (MQ=255)
                                                                                     ggTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATAGGTTTAGAGTCCAGCGa    <  2:4733373/93‑1 (MQ=255)
                                                                                     ggTGTGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGa    <  1:3207477/93‑1 (MQ=255)
                                                                                      gtgtGCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGAc   <  2:1770877/93‑1 (MQ=255)
                                                                                       tgtgCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACg  >  2:2877150/1‑93 (MQ=255)
                                                                                       tgtgCGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACg  >  2:6357134/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCCCCACTTTTTTGCCTTTCAGATCATCGGCTGTTTTAATGGTGCCTTCGTTACCTTTTTTCACCAGCGCCTGAATACCAGAAATGGTGTACGGGGTTGAGAAATCGTATTTTTTCTTGCGCTCATCAGAAATGGTGACCTGATTAATCACCACATCAATACGTTTAGAGTCCAGCGACG  >  NC_000913/1997949‑1998128

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