Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,666,143 A→G I258I (ATT→ATC hcaR ← HcaR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,666,1430AG98.8% 541.0 / ‑3.7 169I258I (ATT→ATChcaRHcaR
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (52/115);  total (53/116)
Fisher's exact test for biased strand distribution p-value = 5.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 172 total aligned reads displayed.

CTGCACAATGGCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCAT  >  NC_000913/2666053‑2666232
                                                                                          |                                                                                         
cTGCACAATGGCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGAt                                                                                         <  1:969839/93‑1 (MQ=255)
 tGCACAATGGCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCTATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATg                                                                                        <  2:4192396/93‑1 (MQ=255)
  gCACAATGGCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATgg                                                                                       <  1:112682/93‑1 (MQ=255)
  gCACAATGGCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATgg                                                                                       <  2:1272709/93‑1 (MQ=255)
  gCACAATGGCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATgg                                                                                       <  1:3930439/93‑1 (MQ=255)
   cacaATGGCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATggg                                                                                      <  2:3179327/93‑1 (MQ=255)
    acaATGGCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGa                                                                                     <  2:3854545/93‑1 (MQ=255)
    acaATGGCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGa                                                                                     <  2:1978418/93‑1 (MQ=255)
     caATGGCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGAc                                                                                    <  2:2669672/93‑1 (MQ=255)
     caATGGCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGAc                                                                                    >  1:2840296/1‑93 (MQ=255)
     caATGGCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGAc                                                                                    >  1:2840295/1‑93 (MQ=255)
         ggCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGaaa                                                                                <  2:3669727/93‑1 (MQ=255)
          gCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGaaaa                                                                               <  1:1764243/93‑1 (MQ=255)
          gCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGaaaa                                                                               <  1:3544686/93‑1 (MQ=255)
           cGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGaaaaa                                                                              <  2:4682698/93‑1 (MQ=255)
            gATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAaac                                                                             >  2:757622/1‑93 (MQ=255)
            gATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAaac                                                                             >  2:4525350/1‑93 (MQ=255)
             aTGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACTCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAaaca                                                                            >  1:3937568/1‑93 (MQ=255)
               gAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAaacaac                                                                          <  2:4621000/93‑1 (MQ=255)
                aaaaCGCTTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAAcc                                                                         <  2:3978977/93‑1 (MQ=255)
                aaaTCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAAcc                                                                         <  1:3646654/93‑1 (MQ=255)
                aaaTCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAAcc                                                                         >  1:1540216/1‑93 (MQ=255)
                aaaTCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAAcc                                                                         <  1:2715486/93‑1 (MQ=255)
                 aaTCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCt                                                                        <  1:4177929/93‑1 (MQ=255)
                  aTCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTg                                                                       >  2:1753559/1‑93 (MQ=255)
                  aTCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTg                                                                       <  1:613256/93‑1 (MQ=255)
                  aTCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTg                                                                       >  1:4618487/1‑93 (MQ=255)
                  aTCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTg                                                                       >  2:467799/1‑93 (MQ=255)
                  aTCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTg                                                                       >  2:571747/1‑93 (MQ=255)
                   tCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGt                                                                      <  1:860791/93‑1 (MQ=255)
                   tCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAAAAACCTGt                                                                      >  1:2423174/1‑93 (MQ=255)
                    cgcgTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTc                                                                     <  2:485990/93‑1 (MQ=255)
                     gcgTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTcc                                                                    >  2:3021069/1‑93 (MQ=255)
                      cgTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCg                                                                   <  2:4124269/93‑1 (MQ=255)
                        tAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGt                                                                 <  1:1752801/93‑1 (MQ=255)
                        tAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGt                                                                 <  2:1378496/93‑1 (MQ=255)
                        tAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGt                                                                 <  1:3460142/93‑1 (MQ=255)
                         aaTGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTa                                                                <  1:2138443/93‑1 (MQ=255)
                          aTGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTAt                                                               <  1:1938000/93‑1 (MQ=255)
                           tGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTAtt                                                              <  1:1934945/93‑1 (MQ=255)
                            gCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTa                                                             <  1:3583891/93‑1 (MQ=255)
                             ccGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTaa                                                            <  2:704319/93‑1 (MQ=255)
                               ggCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTaaaa                                                          <  1:2067990/93‑1 (MQ=255)
                                 cTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAACatt                                                        >  1:3719138/1‑93 (MQ=255)
                                 cTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAatt                                                        >  1:3719137/1‑93 (MQ=255)
                                  ttCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAatta                                                       <  2:2286575/93‑1 (MQ=255)
                                  ttCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAatta                                                       >  2:2718099/1‑93 (MQ=255)
                                   tCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAattat                                                      <  2:901133/93‑1 (MQ=255)
                                   tCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAattat                                                      <  1:3403132/93‑1 (MQ=255)
                                     ctctcCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTc                                                    <  2:4556423/92‑1 (MQ=255)
                                      tctcCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCa                                                   <  1:1206723/93‑1 (MQ=255)
                                      tctcCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCa                                                   >  1:1350610/1‑93 (MQ=255)
                                       ctcCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCat                                                  >  1:1471591/1‑93 (MQ=255)
                                        tcCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCata                                                 >  2:2377053/1‑93 (MQ=255)
                                        tcCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCata                                                 <  2:2840642/93‑1 (MQ=255)
                                          cccTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCAtata                                               >  1:312799/1‑93 (MQ=255)
                                           ccTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATaa                                              <  1:2958143/93‑1 (MQ=255)
                                            cTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAAc                                             <  2:162925/93‑1 (MQ=255)
                                            cTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAAc                                             <  2:166062/93‑1 (MQ=255)
                                             tttCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAAcc                                            >  2:4013349/1‑93 (MQ=255)
                                             tttCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAAcc                                            >  2:1127300/1‑93 (MQ=255)
                                             tttCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAAcc                                            >  2:1127299/1‑93 (MQ=255)
                                               tcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCgg                                          >  2:145468/1‑93 (MQ=255)
                                               tcttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCgg                                          >  2:145469/1‑93 (MQ=255)
                                                 ttcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGt                                        <  1:324610/93‑1 (MQ=255)
                                                  tcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTa                                       <  2:4398886/93‑1 (MQ=255)
                                                  tcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTa                                       <  1:2378181/93‑1 (MQ=255)
                                                  tcCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTa                                       <  1:2286380/93‑1 (MQ=255)
                                                    cACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATc                                     >  2:644535/1‑93 (MQ=255)
                                                    cACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATc                                     <  1:3314771/93‑1 (MQ=255)
                                                       gCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCaaa                                  <  2:2186428/93‑1 (MQ=255)
                                                       gCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCaaa                                  <  1:4063557/93‑1 (MQ=255)
                                                        ccATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAg                                 >  2:3889677/1‑93 (MQ=255)
                                                         cATCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGt                                >  2:2961285/1‑93 (MQ=255)
                                                          aTCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTg                               <  2:513971/93‑1 (MQ=255)
                                                           tCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGa                              <  2:2406919/93‑1 (MQ=255)
                                                           tCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGa                              <  2:3493533/93‑1 (MQ=255)
                                                           tCAGTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGCATTAAAATTATTCATATAACCGGGTATCAAAGTGa                              >  1:3901388/1‑93 (MQ=255)
                                                              gTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGc                           <  2:4695985/93‑1 (MQ=255)
                                                              gTAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGc                           <  1:2423231/93‑1 (MQ=255)
                                                               tAAAGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGcc                          >  1:921017/1‑93 (MQ=255)
                                                                aaaGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGccc                         <  1:4455977/93‑1 (MQ=255)
                                                                 aaGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCa                        <  2:2961910/93‑1 (MQ=255)
                                                                 aaGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCa                        <  1:3600184/93‑1 (MQ=255)
                                                                  aGCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAg                       >  1:2339992/1‑93 (MQ=255)
                                                                   gCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGc                      >  1:1893503/1‑93 (MQ=255)
                                                                   gCAATGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGc                      <  2:2653645/93‑1 (MQ=255)
                                                                     aaTGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGccc                    <  2:601434/93‑1 (MQ=255)
                                                                       tGGAAGGCACATTTCCGGCGATGGGACGATAAACACCCTGTCCGGTATTAACATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCa                  <  1:403111/93‑1 (MQ=255)
                                                                       tGGAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCa                  <  1:403109/93‑1 (MQ=255)
                                                                       aggAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCa                  <  1:3387506/92‑1 (MQ=255)
                                                                        ggAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCAt                 >  1:1831042/1‑93 (MQ=255)
                                                                        ggAAGGCACATTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCAt                 <  1:3778295/93‑1 (MQ=255)
                                                                               acaTTTCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCccacc          <  1:3235624/93‑1 (MQ=255)
                                                                                  tttCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGa       <  1:554168/93‑1 (MQ=255)
                                                                                   ttCCGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGAt      <  1:1820416/93‑1 (MQ=255)
                                                                                      cGGCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCa   <  2:3553460/93‑1 (MQ=255)
                                                                                       ggCGATGGGACTAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCAt  <  1:2485374/93‑1 (MQ=255)
                                                                                       ggCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCAt  <  2:4512639/93‑1 (MQ=255)
                                                                                       ggCGATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCAt  >  1:4198817/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTGCACAATGGCGATGAAATCGCGTAATGCCGGCTTCATCTCCCCTTTCTTCCACGCCATCAGTAAAGCAATGGAAGGCACATTTCCGGCAATGGGACGAAAAACAACCTGTCCGGTATTAAAATTATTCATATAACCGGGTATCAAAGTGACGCCCAGCCCCATGCCCACCAGATTCAT  >  NC_000913/2666053‑2666232

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: