Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,676,230 T→C T197T (ACA→ACG yhjH ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,676,2300TC96.8% 826.4 / 8.2 279T197T (ACA→ACGyhjHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (7/2);  new base C (141/129);  total (148/131)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 292 total aligned reads displayed.

AACATCACGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATT  >  NC_000913/3676140‑3676322
                                                                                          |                                                                                            
aaCATCACGCCCCTCTTCCGGCGTTTCTACGCCCTCGACAATCGCCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAgcgt                                                                                            <  2:1260663/93‑1 (MQ=255)
      ccGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTcc                                                                                      <  1:4288654/92‑1 (MQ=255)
      aCGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTcc                                                                                      >  1:1190570/1‑93 (MQ=255)
      aCGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTcc                                                                                      >  2:3684532/1‑93 (MQ=255)
       cGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCt                                                                                     >  1:999055/1‑93 (MQ=255)
        gCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCtt                                                                                    <  1:323642/93‑1 (MQ=255)
          cACTCTTCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTcc                                                                                  <  2:3830416/93‑1 (MQ=255)
           aCTCTTCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCg                                                                                 <  2:3665022/93‑1 (MQ=255)
            ctctTCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCgg                                                                                <  1:4136194/93‑1 (MQ=255)
               ttCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGa                                                                             >  1:1451592/1‑93 (MQ=255)
               ttCCGGAGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCATGCGTCCTTCCGGcca                                                                             >  2:3617766/1‑91 (MQ=255)
                tCCGGCGTTGCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGAc                                                                            <  1:1513918/93‑1 (MQ=255)
                 ccGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACt                                                                           >  2:1536807/1‑93 (MQ=255)
                  cGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTg                                                                          <  1:1728623/93‑1 (MQ=255)
                    gCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGAc                                                                        >  1:3106465/1‑93 (MQ=255)
                    gCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGAc                                                                        >  2:522485/1‑93 (MQ=255)
                      gTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACgc                                                                      >  1:3964271/1‑93 (MQ=255)
                      gTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACgc                                                                      <  1:2110690/93‑1 (MQ=255)
                      gTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACgc                                                                      >  1:4376644/1‑93 (MQ=255)
                       tttCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACgca                                                                     >  2:3102724/1‑93 (MQ=255)
                       tttCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACgca                                                                     >  1:2186303/1‑93 (MQ=255)
                       tttCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACgca                                                                     >  1:939055/1‑93 (MQ=255)
                       tttCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACgca                                                                     <  2:4054016/93‑1 (MQ=255)
                        ttCTCCGCCCTCTACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACgcag                                                                    <  1:1063444/93‑1 (MQ=255)
                         tCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACgcagc                                                                   <  2:3392596/93‑1 (MQ=255)
                           tACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCAt                                                                 <  1:1298218/93‑1 (MQ=255)
                            aCGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATc                                                                >  1:4677725/1‑93 (MQ=255)
                            aCGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATc                                                                <  1:1557258/93‑1 (MQ=255)
                              gCCCTCGACAACCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATcac                                                              <  2:2818042/93‑1 (MQ=255)
                                ccTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACaa                                                            <  2:157777/93‑1 (MQ=255)
                                 cTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACaaa                                                           <  1:2633440/93‑1 (MQ=255)
                                 cTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACaaa                                                           <  2:1812113/93‑1 (MQ=255)
                                   cGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACa                                                         <  2:601928/93‑1 (MQ=255)
                                      cAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGtt                                                      >  2:144260/1‑93 (MQ=255)
                                      cAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGtt                                                      <  2:2880780/93‑1 (MQ=255)
                                       aaTCCCCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTtc                                                     <  2:2433498/93‑1 (MQ=255)
                                       aaTCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTtc                                                     >  1:2844779/1‑93 (MQ=255)
                                        aTCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTtct                                                    >  2:4412622/1‑93 (MQ=255)
                                        aTCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTtct                                                    <  2:3408583/93‑1 (MQ=255)
                                         tCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTtctc                                                   <  1:3687741/93‑1 (MQ=255)
                                           aCCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTcgc                                                 >  2:4239148/1‑93 (MQ=255)
                                            ccccGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTcgcg                                                >  2:719522/1‑93 (MQ=255)
                                             cccGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTcgcgc                                               <  2:690170/93‑1 (MQ=255)
                                              ccGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTcgcgcg                                              >  2:302664/1‑93 (MQ=255)
                                              ccGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTcgcgcg                                              >  2:1836264/1‑93 (MQ=255)
                                                gcgACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGAt                                            <  2:2348908/93‑1 (MQ=255)
                                                  gACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGAttt                                          <  1:590270/93‑1 (MQ=255)
                                                  gACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGAttt                                          >  2:3444097/1‑93 (MQ=255)
                                                   gcAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGAtttt                                         <  1:3550866/92‑1 (MQ=255)
                                                    cAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGTGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTg                                        >  1:2516927/1‑93 (MQ=255)
                                                    cAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTg                                        <  1:1568212/93‑1 (MQ=255)
                                                     aaTAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGa                                       <  1:760525/93‑1 (MQ=255)
                                                      aTAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGAt                                      >  1:3439675/1‑93 (MQ=255)
                                                       tAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATg                                     >  1:4066644/1‑93 (MQ=255)
                                                       tAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATg                                     <  1:2443044/93‑1 (MQ=255)
                                                       tAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATg                                     >  2:4516915/1‑93 (MQ=255)
                                                        aGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGt                                    >  2:3183928/1‑93 (MQ=255)
                                                        aGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGt                                    <  2:63901/93‑1 (MQ=255)
                                                        aGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGt                                    <  1:1245682/93‑1 (MQ=255)
                                                          cGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAg                                  <  2:752506/93‑1 (MQ=255)
                                                          cGATTCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAg                                  <  2:1297624/93‑1 (MQ=255)
                                                             ttCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCa                               <  2:852114/93‑1 (MQ=255)
                                                             ttCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCa                               >  2:4154234/1‑93 (MQ=255)
                                                             ttCATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGACGTAGTCa                               >  1:1924425/1‑93 (MQ=255)
                                                               cATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATa                             <  1:1952535/93‑1 (MQ=255)
                                                               cATTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATa                             >  2:1694692/1‑93 (MQ=255)
                                                                cttAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATaa                            <  2:3708877/92‑1 (MQ=255)
                                                                aTTAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATaa                            <  1:4421347/93‑1 (MQ=255)
                                                                 ttAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAAc                           <  2:2006053/93‑1 (MQ=255)
                                                                 ttAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAAc                           <  1:2346811/93‑1 (MQ=255)
                                                                 ttAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAAc                           <  2:4211111/93‑1 (MQ=255)
                                                                 ttAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAAc                           >  2:2604830/1‑93 (MQ=255)
                                                                 ttAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAAc                           >  2:1328859/1‑93 (MQ=255)
                                                                 ttAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAAc                           >  2:3961113/1‑93 (MQ=255)
                                                                 ttAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAAc                           >  1:4250214/1‑93 (MQ=255)
                                                                 ttAGATGTAAAAGCTAAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAAc                           <  2:2147936/93‑1 (MQ=255)
                                                                  tAGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACg                          >  2:2552751/1‑93 (MQ=255)
                                                                   aGATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGc                         >  2:1753967/1‑93 (MQ=255)
                                                                    gATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGca                        >  2:3285988/1‑93 (MQ=255)
                                                                    gATGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGca                        >  2:1171594/1‑93 (MQ=255)
                                                                     aTGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGcac                       <  1:3309864/93‑1 (MQ=255)
                                                                      tgaaaaaGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGcact                      >  1:2010901/4‑93 (MQ=255)
                                                                      tGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGcact                      >  2:1964591/1‑93 (MQ=255)
                                                                      tGTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGcact                      >  1:2936135/1‑93 (MQ=255)
                                                                       gTAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGcactt                     <  1:871243/93‑1 (MQ=255)
                                                                        tAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGcacttc                    <  1:2224967/93‑1 (MQ=255)
                                                                        tAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGcacttc                    <  2:3989760/93‑1 (MQ=255)
                                                                        tAAAAGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGcacttc                    <  1:3172968/93‑1 (MQ=255)
                                                                         aaaaGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGcacttca                   >  1:4352816/1‑93 (MQ=255)
                                                                         aaaaGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGcacttca                   <  1:1854031/93‑1 (MQ=255)
                                                                          aaaGCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGcacttcac                  <  2:1791609/93‑1 (MQ=255)
                                                                             gCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTa               <  1:1032970/93‑1 (MQ=255)
                                                                             gCTGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTa               <  2:3500301/93‑1 (MQ=255)
                                                                               tGAGAAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAgc             >  2:472083/1‑93 (MQ=255)
                                                                                  gaAAAGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCa          >  2:2298335/1‑93 (MQ=255)
                                                                                    aaaGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCaga        <  2:2459626/93‑1 (MQ=255)
                                                                                     aaGAGCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCagag       <  1:2236652/93‑1 (MQ=255)
                                                                                      agagTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAgaga      >  1:4158245/1‑93 (MQ=255)
                                                                                       gagCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGaa     >  1:2665951/1‑93 (MQ=255)
                                                                                       gagCGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGaa     >  1:3800104/1‑93 (MQ=255)
                                                                                          |                                                                                            
AACATCACGCCACTCTTCCGGCGTTTCTACGCCCTCGACAATCACCCCGCGACAATAGCGATTCATTAGATGTAAAAGCTGAGAAAAGAGTGTGCGTCCTTCCGGCGACTGACGCAGCATCACAAACAGTTCTCGCGCGATTTTGATGTAGTCATAACGCACTTCACTTAGCGCAGAGAAATT  >  NC_000913/3676140‑3676322

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