Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,519,062 A→G I176T (ATA→ACA)  yjhF ← putative transport system permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,519,0620AG99.0% 646.5 / ‑3.2 199I176T (ATA→ACA) yjhFputative transport system permease
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (119/78);  total (121/78)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 206 total aligned reads displayed.

TTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAGCAAA  >  NC_000913/4518972‑4519151
                                                                                          |                                                                                         
 tAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCa                                                                                        >  1:3018601/1‑93 (MQ=255)
 tAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCa                                                                                        >  1:4652020/1‑93 (MQ=255)
  aaaGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAg                                                                                       >  2:2171511/1‑93 (MQ=255)
  aaaGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAg                                                                                       <  2:2458528/93‑1 (MQ=255)
   aaGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGt                                                                                      >  2:2518936/1‑93 (MQ=255)
     ggCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTAc                                                                                    >  1:1791079/1‑93 (MQ=255)
     ggCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTAc                                                                                    <  2:3852632/93‑1 (MQ=255)
        gTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCaa                                                                                 >  1:2663838/1‑93 (MQ=255)
        gTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCaa                                                                                 >  1:3494397/1‑93 (MQ=255)
         tCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAAt                                                                                >  2:3155414/1‑93 (MQ=255)
         tCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAAt                                                                                <  1:4143684/93‑1 (MQ=255)
         tCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAAt                                                                                >  2:2086732/1‑93 (MQ=255)
          gcGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATc                                                                               <  2:1750428/92‑1 (MQ=255)
          ccGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATc                                                                               <  1:4289151/93‑1 (MQ=255)
          ccGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATc                                                                               >  1:1121860/1‑93 (MQ=255)
           cGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATcc                                                                              >  2:2034325/1‑93 (MQ=255)
           cGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATcc                                                                              >  1:2286842/1‑93 (MQ=255)
             gAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTa                                                                            >  1:1857602/1‑93 (MQ=255)
              aataatGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTAc                                                                           <  1:2039488/93‑1 (MQ=255)
              aataatGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTAc                                                                           >  1:1050245/1‑93 (MQ=255)
                 aatGCGCTGGCAGAATTTATTCAATACCGGTCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATc                                                                        <  1:2648680/93‑1 (MQ=255)
                  atGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCa                                                                       <  2:1767435/93‑1 (MQ=255)
                  atGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCa                                                                       >  1:1339032/1‑93 (MQ=255)
                   tGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAg                                                                      <  2:3720752/93‑1 (MQ=255)
                   tGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAg                                                                      >  2:3971031/1‑93 (MQ=255)
                       cTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGcc                                                                  <  2:2959299/93‑1 (MQ=255)
                         ggCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCTGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCat                                                                >  1:4403238/1‑93 (MQ=255)
                         ggCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCat                                                                >  1:861402/1‑93 (MQ=255)
                         ggCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCat                                                                >  1:4475054/1‑93 (MQ=255)
                         ggCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCat                                                                >  2:3908714/1‑93 (MQ=255)
                         ggCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCat                                                                <  1:2151472/93‑1 (MQ=255)
                          gCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCata                                                               <  2:458940/93‑1 (MQ=255)
                           cAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCatat                                                              >  1:1299957/1‑93 (MQ=255)
                           cAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCatat                                                              <  1:2946509/93‑1 (MQ=255)
                             gAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTc                                                            >  2:4302445/1‑93 (MQ=255)
                             gAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTc                                                            >  1:2271355/1‑93 (MQ=255)
                             gAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTc                                                            <  1:2928174/93‑1 (MQ=255)
                               aTTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTccc                                                          <  2:3529209/93‑1 (MQ=255)
                               aTTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTccc                                                          <  2:142068/93‑1 (MQ=255)
                                tttATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCt                                                         >  1:1394597/1‑93 (MQ=255)
                                 ttATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTg                                                        <  1:1894767/93‑1 (MQ=255)
                                 ttATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTg                                                        >  1:4643394/1‑93 (MQ=255)
                                  tATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCACTCCTACATCAGCGCCATATTCCCtac                                                       >  2:959135/1‑91 (MQ=255)
                                   aTTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGc                                                      <  2:1584476/93‑1 (MQ=255)
                                     tCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTTGCGa                                                    >  2:1729049/1‑93 (MQ=255)
                                      cAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGAt                                                   <  1:3214301/93‑1 (MQ=255)
                                       aaTACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATa                                                  <  1:2709693/93‑1 (MQ=255)
                                       aaTACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATa                                                  <  2:1270867/93‑1 (MQ=255)
                                         tACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAAc                                                >  2:3210979/1‑93 (MQ=255)
                                         tACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAAc                                                >  1:1724827/1‑93 (MQ=255)
                                          aCCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACa                                               <  1:534997/93‑1 (MQ=255)
                                           ccGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAg                                              >  2:821627/1‑93 (MQ=255)
                                            cGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGt                                             <  1:4611491/93‑1 (MQ=255)
                                             gggCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTc                                            <  2:435137/93‑1 (MQ=255)
                                             gggCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTc                                            <  2:686262/93‑1 (MQ=255)
                                             gggCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTc                                            <  1:1801008/93‑1 (MQ=255)
                                              ggCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCg                                           <  2:4172854/93‑1 (MQ=255)
                                               gCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCgg                                          <  2:1593981/93‑1 (MQ=255)
                                                ccGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCggc                                         >  1:3870071/1‑92 (MQ=255)
                                                 cGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGc                                        <  2:657765/93‑1 (MQ=255)
                                                 cGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGc                                        >  1:246576/1‑93 (MQ=255)
                                                  gCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGcc                                       >  1:1631314/1‑93 (MQ=255)
                                                  gCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGcc                                       >  1:1779617/1‑93 (MQ=255)
                                                   cAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCg                                      <  1:3161033/93‑1 (MQ=255)
                                                   cAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCg                                      >  2:1518281/1‑93 (MQ=255)
                                                   cAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCAATATTCCCTGGCGATAAAAGCAGGGCCg                                      >  1:1604708/1‑93 (MQ=255)
                                                    aaaGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCgg                                     >  2:3923898/1‑93 (MQ=255)
                                                    aaaGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCgg                                     <  2:504893/93‑1 (MQ=255)
                                                    aaaGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCgg                                     >  2:560302/1‑93 (MQ=255)
                                                    aaaGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCgg                                     <  1:2535573/93‑1 (MQ=255)
                                                    aaaGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCgg                                     >  2:59308/1‑93 (MQ=255)
                                                     aaGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCggg                                    >  2:884516/1‑93 (MQ=255)
                                                      aGAATAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCgggg                                   >  2:1188220/1‑93 (MQ=255)
                                                        aaTAAAGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTg                                 <  1:2852797/93‑1 (MQ=255)
                                                           aaaGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTgcgg                              >  1:1236769/1‑93 (MQ=255)
                                                            aaGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTgcggc                             <  1:3338669/93‑1 (MQ=255)
                                                             aGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTgcggcg                            >  1:2985984/1‑93 (MQ=255)
                                                             aGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTgcggcg                            >  2:2639493/1‑93 (MQ=255)
                                                             aGGTTGGAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTgcggcg                            >  2:557615/1‑93 (MQ=255)
                                                                  ggAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCaaa                       >  2:746129/1‑93 (MQ=255)
                                                                  ggAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCaaa                       >  2:2573778/1‑93 (MQ=255)
                                                                   gAATGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCaaaa                      >  1:3534313/1‑93 (MQ=255)
                                                                    aaTGCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCaaaaa                     >  2:487136/1‑93 (MQ=255)
                                                                       gCCGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAcc                  <  1:3122604/93‑1 (MQ=255)
                                                                        ccGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAcct                 >  1:4046480/1‑92 (MQ=255)
                                                                        ccGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCg                 >  2:3064458/1‑93 (MQ=255)
                                                                         cGACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCgt                >  2:1773710/1‑93 (MQ=255)
                                                                          gACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAAGCgtg               >  2:2873640/1‑93 (MQ=255)
                                                                          gACAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCgtg               <  1:916579/93‑1 (MQ=255)
                                                                           aCAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCgtgt              <  2:2297973/93‑1 (MQ=255)
                                                                            cAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGtgtg             >  1:395636/1‑93 (MQ=255)
                                                                            cAATGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGtgtg             <  1:4098895/93‑1 (MQ=255)
                                                                              aTGATACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCg           <  2:4661531/93‑1 (MQ=255)
                                                                                 aTACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGAcc        <  2:3476538/93‑1 (MQ=255)
                                                                                  tACCGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGGCGGGGTGCGGCGGCAAAACACCGTGTGCGACCa       >  1:4732127/1‑93 (MQ=255)
                                                                                    ccGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAGc     >  2:4033205/1‑93 (MQ=255)
                                                                                     cGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAGCa    >  1:3730529/1‑93 (MQ=255)
                                                                                     cGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAGCa    >  2:1120592/1‑93 (MQ=255)
                                                                                     cGTATGTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAGCa    >  1:1014881/1‑93 (MQ=255)
                                                                                       tataTCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAGCaaa  <  1:2703927/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATATTCCCTGGCGATAACAGTCGGGCCGGGGTGCGGCGGCAAAAAACCGTGTGCGACCAGCAAA  >  NC_000913/4518972‑4519151

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