Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 751,748 C→T G91R (GGA→AGA)  ybgD ← putative fimbrial‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913751,7480CT99.6% 880.5 / ‑4.9 225G91R (GGA→AGA) ybgDputative fimbrial‑like protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (99/125);  total (100/125)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

Only 100 of 226 total aligned reads displayed.

CGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGG  >  NC_000913/751658‑751837
                                                                                          |                                                                                         
 gCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCt                                                                                        >  2:4147825/1‑93 (MQ=255)
 gCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCt                                                                                        <  2:6468905/93‑1 (MQ=255)
 cccTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCt                                                                                        <  1:1120705/92‑1 (MQ=255)
    tGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTAcc                                                                                     >  2:1733684/1‑93 (MQ=255)
    tGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTAcc                                                                                     >  2:693192/1‑93 (MQ=255)
       aCTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAATGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGtt                                                                                  <  1:3179286/93‑1 (MQ=255)
       aCTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCCTTCTGCTACCGtt                                                                                  <  1:6200067/93‑1 (MQ=255)
        cTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTg                                                                                 >  2:4686250/1‑93 (MQ=255)
        cTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTg                                                                                 >  2:4887123/1‑93 (MQ=255)
         tGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGt                                                                                <  2:4941654/93‑1 (MQ=255)
         tGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGt                                                                                <  2:5588608/93‑1 (MQ=255)
          ggTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTc                                                                               <  2:5594358/93‑1 (MQ=255)
          ggTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTc                                                                               <  2:2987314/93‑1 (MQ=255)
          ggTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTc                                                                               <  2:4415722/93‑1 (MQ=255)
              ttGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGaa                                                                           <  1:1151022/93‑1 (MQ=255)
              ttGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGaa                                                                           <  1:3957060/93‑1 (MQ=255)
                gCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAgc                                                                         >  2:3474960/1‑93 (MQ=255)
                gCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAgc                                                                         <  1:1384595/93‑1 (MQ=255)
                    aacaaAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAgcaggc                                                                     <  1:3357736/93‑1 (MQ=255)
                    aacaaAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAgcaggc                                                                     >  2:5446169/1‑93 (MQ=255)
                    aacaaAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAgcaggc                                                                     >  1:1793934/1‑93 (MQ=255)
                     acaaAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAgcaggca                                                                    <  2:3589138/93‑1 (MQ=255)
                      caaAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAgcaggcag                                                                   >  1:1829509/1‑93 (MQ=255)
                        aaGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGAt                                                                 <  1:2628182/93‑1 (MQ=255)
                          ggCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATca                                                               <  2:5476220/93‑1 (MQ=255)
                            cGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATcaca                                                             <  1:461225/93‑1 (MQ=255)
                             gTAGCGCCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAg                                                            >  1:4186743/1‑93 (MQ=255)
                             gTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAg                                                            <  1:900064/93‑1 (MQ=255)
                             gTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAg                                                            >  2:1246518/1‑93 (MQ=255)
                              tAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGt                                                           <  1:4084512/93‑1 (MQ=255)
                              tAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGt                                                           <  2:6376221/93‑1 (MQ=255)
                               aGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGtt                                                          <  2:1088935/93‑1 (MQ=255)
                               aGCACCAGTTGTCTTAGCCGTGCTATCGAAGGGTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGtt                                                          >  1:5924029/1‑93 (MQ=255)
                                gCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTg                                                         <  1:3402859/93‑1 (MQ=255)
                                 cACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGa                                                        <  2:4625846/93‑1 (MQ=255)
                                 cACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGa                                                        >  2:6070402/1‑93 (MQ=255)
                                 cACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGa                                                        <  1:3301330/93‑1 (MQ=255)
                                  aCCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGAt                                                       >  1:5851382/1‑93 (MQ=255)
                                  aCCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGAt                                                       <  2:4104952/93‑1 (MQ=255)
                                  aCCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGAt                                                       >  2:1131392/1‑93 (MQ=255)
                                  aCCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGAt                                                       >  1:2903063/1‑93 (MQ=255)
                                   ccAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATc                                                      >  1:4437915/1‑93 (MQ=255)
                                    cAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCa                                                     >  1:5043854/1‑93 (MQ=255)
                                      gTTGTCTTAGCCGTGCTATCGTAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAg                                                   <  2:5657601/93‑1 (MQ=255)
                                      gTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAg                                                   >  1:5072953/1‑93 (MQ=255)
                                       ttGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGc                                                  >  1:2164563/1‑93 (MQ=255)
                                        tGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCg                                                 <  1:5787092/93‑1 (MQ=255)
                                         gTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGa                                                >  2:5412321/1‑93 (MQ=255)
                                          tCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGaa                                               >  1:2276259/1‑93 (MQ=255)
                                          tCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATCCTGCTACCGTTGTCAGAAGCAGTCAGATCACAGTTGATCAAGCGaa                                               <  1:5182230/93‑1 (MQ=255)
                                           cTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAAt                                              <  1:5979355/93‑1 (MQ=255)
                                           cTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAAt                                              <  1:5339938/93‑1 (MQ=255)
                                           cTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAAt                                              <  1:3029976/93‑1 (MQ=255)
                                            ttAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATg                                             >  2:4722181/1‑93 (MQ=255)
                                             tAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGt                                            >  1:2513930/1‑93 (MQ=255)
                                               gCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCg                                          >  1:1722438/1‑93 (MQ=255)
                                               gCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCg                                          >  1:6117424/1‑93 (MQ=255)
                                               gCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCg                                          >  1:3779344/1‑93 (MQ=255)
                                                ccGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGa                                         <  1:667929/93‑1 (MQ=255)
                                                 cGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACGGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGAGAATGTCGAc                                        <  2:3824553/93‑1 (MQ=255)
                                                 cGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGAc                                        <  2:4647232/93‑1 (MQ=255)
                                                  gTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACg                                       <  1:6215284/93‑1 (MQ=255)
                                                  gTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACg                                       >  2:5494504/1‑93 (MQ=255)
                                                   tGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACgg                                      <  2:2585304/93‑1 (MQ=255)
                                                     cTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGcc                                    <  2:198172/93‑1 (MQ=255)
                                                      tATCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCa                                   <  2:3137494/93‑1 (MQ=255)
                                                       aTCGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCAc                                  <  2:1926293/93‑1 (MQ=255)
                                                         cGAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACtt                                <  1:3590907/93‑1 (MQ=255)
                                                          gAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACttt                               <  2:4467432/93‑1 (MQ=255)
                                                          gAAGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACttt                               >  2:1299275/1‑93 (MQ=255)
                                                           aaGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTa                              >  2:402915/1‑93 (MQ=255)
                                                           aaGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTa                              <  1:1015984/93‑1 (MQ=255)
                                                            aGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTAt                             <  2:429857/93‑1 (MQ=255)
                                                            aGGTTACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTAt                             >  2:841201/1‑93 (MQ=255)
                                                                tACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCt                         <  1:1515128/93‑1 (MQ=255)
                                                                tACGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCt                         <  2:2858056/93‑1 (MQ=255)
                                                                 aCGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGCCGACGGCCACTTTATTGCtt                        >  1:1662039/1‑92 (MQ=255)
                                                                  cGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTat                       <  2:5943958/93‑1 (MQ=255)
                                                                  cGCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTat                       <  1:2493031/93‑1 (MQ=255)
                                                                   gCCAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatg                      >  2:5079694/1‑93 (MQ=255)
                                                                    ccAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatga                     <  2:3409959/93‑1 (MQ=255)
                                                                    ccAACTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatga                     <  1:483126/93‑1 (MQ=255)
                                                                      aaCTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatgatg                   >  1:5106842/1‑93 (MQ=255)
                                                                      aaCTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTatgatg                   >  1:2048615/1‑93 (MQ=255)
                                                                       aCTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGc                  <  2:753325/93‑1 (MQ=255)
                                                                       aCTTTGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGc                  <  1:2028913/93‑1 (MQ=255)
                                                                         tttGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTc                >  2:2546589/1‑93 (MQ=255)
                                                                          ttGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTcc               >  2:3559539/1‑93 (MQ=255)
                                                                          ttGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTcc               <  2:6453777/93‑1 (MQ=255)
                                                                           tGGATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTccc              <  1:5557387/93‑1 (MQ=255)
                                                                             gATACCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCgg            >  1:5694070/1‑93 (MQ=255)
                                                                                aCCGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGttt         >  2:1093546/1‑93 (MQ=255)
                                                                                 ccGGCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTta        <  2:2788836/93‑1 (MQ=255)
                                                                                    gCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTatat     <  1:1080268/93‑1 (MQ=255)
                                                                                    gCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTatat     >  1:1541595/1‑93 (MQ=255)
                                                                                    gCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTatat     <  2:4637446/93‑1 (MQ=255)
                                                                                    gCATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTatat     >  1:5526622/1‑93 (MQ=255)
                                                                                     cgttCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATg    <  1:4756943/91‑1 (MQ=255)
                                                                                     cATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATg    <  1:190732/93‑1 (MQ=255)
                                                                                     cATTCTGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATg    <  1:4985302/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGCCTGCACTGGTGTTGCTCAACAAAGGCGTAGCACCAGTTGTCTTAGCCGTGCTATCGAAGGTTACGCCAACTTTGGATACCGGCATTCCGCTACCGTTGTCAGAAGCAGGCAGATCACAGTTGATCAAGCGAATGTCGACGGCCACTTTATTGCTATGATGCTCCCGGTTTATATGGG  >  NC_000913/751658‑751837

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