Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 752,007 T→C G4G (GGA→GGG ybgD ← putative fimbrial‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913752,0070TC99.4% 506.4 / ‑4.7 157G4G (GGA→GGGybgDputative fimbrial‑like protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (57/99);  total (57/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

Only 100 of 163 total aligned reads displayed.

CTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGTCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCTAGTTTTTAGTT  >  NC_000913/751917‑752096
                                                                                          |                                                                                         
ctgaaaGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGccc                                                                                         <  2:2365110/90‑1 (MQ=255)
cTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGccc                                                                                         >  2:3494820/1‑93 (MQ=255)
 ttAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCt                                                                                        <  2:1442301/93‑1 (MQ=255)
 ttAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCt                                                                                        >  1:3796043/1‑93 (MQ=255)
   aaaGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCttt                                                                                      >  1:2337250/1‑93 (MQ=255)
      gTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTaaa                                                                                   <  2:2768464/93‑1 (MQ=255)
       tGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAAc                                                                                  <  1:3237833/93‑1 (MQ=255)
       tGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAAc                                                                                  <  1:2745351/93‑1 (MQ=255)
        gAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACa                                                                                 <  1:4654108/93‑1 (MQ=255)
        gAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACa                                                                                 <  1:285450/93‑1 (MQ=255)
        gAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACa                                                                                 <  1:285589/93‑1 (MQ=255)
          aTTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATa                                                                               >  2:2932919/1‑93 (MQ=255)
          aTTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATa                                                                               <  2:705524/93‑1 (MQ=255)
          aTTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATa                                                                               <  1:650060/93‑1 (MQ=255)
           tttCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATaa                                                                              <  1:3686061/93‑1 (MQ=255)
           tttCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATaa                                                                              <  1:1142328/93‑1 (MQ=255)
           tttCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATaa                                                                              >  1:404856/1‑93 (MQ=255)
            ttCGCCAGAACCTTCAGCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGTCCAATGTTTTTTGCCCTTTAAACATAAc                                                                             <  2:424908/93‑1 (MQ=255)
              cGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACag                                                                           <  1:2820244/93‑1 (MQ=255)
               gCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACaga                                                                          <  2:2246651/93‑1 (MQ=255)
               gCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACaga                                                                          <  2:2246506/93‑1 (MQ=255)
               gCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACaga                                                                          <  2:1329255/93‑1 (MQ=255)
                ccAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACagag                                                                         <  1:1774922/93‑1 (MQ=255)
                 cAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACTGCCAGTGCGGCCAATTTTTTTTGCCCTTTAAACATAACAGAGt                                                                        <  1:1155355/93‑1 (MQ=255)
                 cAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGt                                                                        <  2:182584/93‑1 (MQ=255)
                  aGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTc                                                                       >  1:2029632/1‑93 (MQ=255)
                   gAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTcc                                                                      <  2:3917087/93‑1 (MQ=255)
                   gAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTcc                                                                      >  1:3465647/1‑93 (MQ=255)
                    aaCCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCt                                                                     <  2:127556/93‑1 (MQ=255)
                      ccTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCttt                                                                   <  2:2918404/93‑1 (MQ=255)
                       cTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTa                                                                  >  1:1934582/1‑93 (MQ=255)
                        ttCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTaa                                                                 <  2:81458/93‑1 (MQ=255)
                         tcatcaGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAg                                                                <  1:902599/93‑1 (MQ=255)
                         tcatcaGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAg                                                                >  1:755334/1‑93 (MQ=255)
                         tcatcaGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAg                                                                >  1:2075003/1‑93 (MQ=255)
                          catcaGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAgg                                                               <  1:1923255/93‑1 (MQ=255)
                             cagcagCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGata                                                            <  1:852166/93‑1 (MQ=255)
                             cagcagCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGata                                                            <  2:2529282/93‑1 (MQ=255)
                             cagcagCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGata                                                            <  1:540/93‑1 (MQ=255)
                             cagcagCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGata                                                            <  1:4744022/93‑1 (MQ=255)
                              agcagcATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGatat                                                           >  1:3173059/1‑93 (MQ=255)
                              agcagcATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGatat                                                           >  2:4646953/1‑93 (MQ=255)
                               gcagcaTAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGAtata                                                          >  1:3528968/1‑93 (MQ=255)
                               gcagcaTAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGAtata                                                          <  1:3905804/93‑1 (MQ=255)
                               gcagcaTAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGAtata                                                          <  1:2493345/93‑1 (MQ=255)
                                  gcaTAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGaa                                                       >  2:1112098/1‑93 (MQ=255)
                                  gcaTAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGaa                                                       <  1:124996/93‑1 (MQ=255)
                                   caTAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAAt                                                      <  2:3785575/93‑1 (MQ=255)
                                   caTAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAAt                                                      >  2:4551540/1‑93 (MQ=255)
                                    aTAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATa                                                     >  2:4066317/1‑93 (MQ=255)
                                     tAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAg                                                    <  1:1043413/93‑1 (MQ=255)
                                     tAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAg                                                    <  1:3008103/93‑1 (MQ=255)
                                        aCAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAgggg                                                 <  2:4104940/93‑1 (MQ=255)
                                          aGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGta                                               <  1:1672318/93‑1 (MQ=255)
                                           ggTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGtat                                              <  1:1187059/93‑1 (MQ=255)
                                           ggTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGtat                                              >  1:1300839/1‑93 (MQ=255)
                                           ggTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGtat                                              >  2:2074027/1‑93 (MQ=255)
                                                aGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTa                                         >  2:4405842/1‑93 (MQ=255)
                                                 gTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTAc                                        >  1:1076517/1‑93 (MQ=255)
                                                    aaCAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGcc                                     <  2:772041/93‑1 (MQ=255)
                                                     aCAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCa                                    <  2:3405141/93‑1 (MQ=255)
                                                       agcagAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGa                                  <  1:4141438/93‑1 (MQ=255)
                                                         cagAGATACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAat                                <  2:724539/93‑1 (MQ=255)
                                                           gagaTACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAatat                              <  1:2124143/93‑1 (MQ=255)
                                                           gagaTACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAatat                              >  2:1550869/1‑93 (MQ=255)
                                                             gaTACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCg                            <  1:1646533/93‑1 (MQ=255)
                                                             gaTACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCg                            <  1:2119215/93‑1 (MQ=255)
                                                              aTACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGt                           >  2:3000827/1‑93 (MQ=255)
                                                              aTACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGt                           >  1:1063234/1‑93 (MQ=255)
                                                              aTACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGt                           <  1:2432871/93‑1 (MQ=255)
                                                               tACGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTa                          <  1:1640033/93‑1 (MQ=255)
                                                                 cGGCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTAtt                        <  2:1087783/93‑1 (MQ=255)
                                                                  ggCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTAttt                       >  2:1777833/1‑93 (MQ=255)
                                                                   gCCAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTg                      <  1:915787/93‑1 (MQ=255)
                                                                    ccAGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGa                     <  1:2410665/93‑1 (MQ=255)
                                                                      aGTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGatt                   >  2:1278536/1‑93 (MQ=255)
                                                                       gTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGatta                  <  1:446502/93‑1 (MQ=255)
                                                                       gTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGatta                  >  1:808793/1‑93 (MQ=255)
                                                                       gTGCGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGatta                  >  1:808792/1‑93 (MQ=255)
                                                                          cGGCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTg               <  2:1203582/93‑1 (MQ=255)
                                                                           ggCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGc              <  1:3197597/93‑1 (MQ=255)
                                                                           ggCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGc              <  1:2100681/93‑1 (MQ=255)
                                                                           ggCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGc              <  2:2913152/93‑1 (MQ=255)
                                                                           ggCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGc              >  2:2362886/1‑93 (MQ=255)
                                                                            gCCAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCt             <  1:631539/93‑1 (MQ=255)
                                                                             ccAATGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCta            >  2:2509469/1‑93 (MQ=255)
                                                                               aaTGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagt          <  2:2394872/93‑1 (MQ=255)
                                                                                 tGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagttt        >  1:1120466/1‑93 (MQ=255)
                                                                                 tGTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagttt        >  2:4653106/1‑93 (MQ=255)
                                                                                  gTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagtttt       <  1:4380733/93‑1 (MQ=255)
                                                                                  gTTTTTTGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagtttt       <  2:2111171/93‑1 (MQ=255)
                                                                                   ttttttGCCCTTTAAACATAACAGAGTCCTTTAAGTATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagttttt      <  1:2238545/93‑1 (MQ=255)
                                                                                   ttttttGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagttttt      >  2:3149113/1‑93 (MQ=255)
                                                                                   ttttttGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagttttt      >  2:3044844/1‑93 (MQ=255)
                                                                                    tttttGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagttttta     >  2:1638469/1‑93 (MQ=255)
                                                                                    tttttGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagttttta     <  1:2775899/93‑1 (MQ=255)
                                                                                     ttttGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagtttttag    <  1:4173954/93‑1 (MQ=255)
                                                                                      tttGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagtttttagt   <  2:218643/93‑1 (MQ=255)
                                                                                      gttGCCCTTTAAACATAACAGAGTCCTGTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagtttttagt   <  2:3775303/92‑1 (MQ=255)
                                                                                       ttGCCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCtagtttttagtt  >  1:2104673/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTTAAAGTGAATTTCGCCAGAACCTTCATCAGCAGCATAAACAGGTGCAGTGAACAGCAGAGATACGGCCAGTGCGGCCAATGTTTTTTGTCCTTTAAACATAACAGAGTCCTTTAAGGATATAGAATAGGGGTATAGCTACGCCAGAATATCGTATTTGATTATTGCTAGTTTTTAGTT  >  NC_000913/751917‑752096

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