Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,175,525 (C)6→7 coding (402/630 nt) yqiE ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,175,5191.C91.1% 366.3 / 28.6 124coding (408/630 nt)yqiEhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (6/5);  new base C (59/54);  total (65/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

CGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATC  >  NC_000913/3175427‑3175608
                                                                                             |                                                                                            
cgAATATCTTCGTTTTCATCAGCCAGCCCGTAAATCCCGCTTGCGGTCGTGGCGTCCACTTCCCCCACCATAATTGACGAACGCTCACTGGTg                                                                                               <  2:3532503/93‑1 (MQ=255)
cgAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTg                                                                                               >  1:933825/1‑93 (MQ=255)
 gAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑c                                                                                            >  2:2596171/1‑93 (MQ=255)
 gAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑c                                                                                            >  1:2056078/1‑93 (MQ=255)
 gAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑c                                                                                            >  2:2401157/1‑93 (MQ=255)
 gAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑c                                                                                            >  1:893004/1‑93 (MQ=255)
  aaTATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑cc                                                                                          >  2:1518286/1‑93 (MQ=255)
   atatCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑ccc                                                                                        >  1:3944588/1‑93 (MQ=255)
     tttttCGTTTGCATCAGCCAGCCCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑ccccc                                                                                      <  2:1336150/90‑1 (MQ=255)
      tCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑cccccc                                                                                     >  1:1220783/1‑93 (MQ=255)
      tCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑cccccc                                                                                     >  1:1581635/1‑93 (MQ=255)
            ttttTACCAGCCAGCCCGTGAATCCCGTTTGCGGCCGTGGCTTCCGCTTCCCCCACCATAATTGACGACCGCTCACTGGTG‑‑CCACCCCGGGCt                                                                                <  2:4187755/93‑1 (MQ=255)
            ttttCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCCCTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCt                                                                                <  2:2135005/93‑1 (MQ=255)
             tttCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCtt                                                                               >  2:673446/1‑93 (MQ=255)
             ttccatcaGCCAGACCGTGAATCCCGCTTGCGGTCGTGGCGTCCACTCCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCtt                                                                               <  2:289360/90‑1 (MQ=255)
              ttCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTg                                                                              >  2:1198424/1‑93 (MQ=255)
                catcaGCCAGACGGTGAATCCCGCTTGCGGTCGTGCCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGcc                                                                            <  1:3789776/93‑1 (MQ=255)
                 atcaGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTGGCCCACCATAATTAACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCa                                                                           >  1:4392952/1‑93 (MQ=255)
                 atcaGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCCGGGCTTGcc                                                                            >  2:1334090/1‑93 (MQ=255)
                 atcaGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCa                                                                           >  1:686821/1‑93 (MQ=255)
                   caGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAgg                                                                         >  1:1721632/1‑93 (MQ=255)
                    aGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGa                                                                        <  2:84430/93‑1 (MQ=255)
                     gCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGaa                                                                       <  2:2217425/93‑1 (MQ=255)
                      ccAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGaaa                                                                      >  1:976150/1‑93 (MQ=255)
                       cAGCCCGTGAATACCCCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAAc                                                                     <  1:3472898/93‑1 (MQ=255)
                       cAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAAc                                                                     >  1:2484657/1‑93 (MQ=255)
                        aGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACt                                                                    >  1:1621417/1‑93 (MQ=255)
                        aGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGCG‑‑CTCCCCCGGGCGTGCCAGGTAACt                                                                    >  1:1869164/1‑93 (MQ=255)
                        aGACCGTGAATACCGCTTGCGGTCGGGGCGTGCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACt                                                                    <  2:2679898/93‑1 (MQ=255)
                         gACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACtt                                                                   >  2:1036060/1‑93 (MQ=255)
                         gACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACtt                                                                   >  2:2780559/1‑93 (MQ=255)
                           cggTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTaa                                                                 <  2:3071486/91‑1 (MQ=255)
                             atGAATACCGAGTGCGGGCGGGGCTTCCAATGCGCCCGCCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAaca                                                               <  2:251010/92‑1 (MQ=255)
                              tGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAacac                                                              >  1:3370713/1‑93 (MQ=255)
                                aaTACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGACCCCTCACTCGTG‑‑CACCCCCGGGCTTGCCAGGACACTTAACACCg                                                            >  1:2301924/1‑93 (MQ=255)
                                aaTACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCTGTTGC‑‑CCCCCCCGGCGTGCCCGGAAACTTAACACCCGt                                                          >  1:519071/1‑93 (MQ=255)
                                aaTACCGCTGGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCg                                                            <  2:1890062/93‑1 (MQ=255)
                                  tcccGCTTGCGGTCGTGGCGTCCCCTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGt                                                          <  2:330254/91‑1 (MQ=255)
                                  tACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGt                                                          >  2:3302242/1‑93 (MQ=255)
                                   aCCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACAACGGtt                                                         >  1:4676507/1‑93 (MQ=255)
                                     cGCTTGCGGCCCGGGCCTCCCCTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTg                                                       <  1:349627/93‑1 (MQ=255)
                                       cTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGt                                                     >  1:106984/1‑93 (MQ=255)
                                       cTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGt                                                     >  1:106983/1‑93 (MQ=255)
                                        ttGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTc                                                    <  1:4418186/93‑1 (MQ=255)
                                          gCGGCCGGGGCGCCCACTCCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCg                                                  <  2:3797190/93‑1 (MQ=255)
                                            ggTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGtt                                                <  1:749958/93‑1 (MQ=255)
                                              gcGTGGCGTCCCCTCCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTg                                              <  2:3882261/92‑1 (MQ=255)
                                                gTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGAc                                            >  2:3740420/1‑93 (MQ=255)
                                                 tGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACt                                           <  1:4476797/93‑1 (MQ=255)
                                                  ggCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTa                                          <  1:4627500/93‑1 (MQ=255)
                                                   gCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTAt                                         >  2:3372259/1‑93 (MQ=255)
                                                     gTCCCCTCCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCa                                       <  1:1172711/93‑1 (MQ=255)
                                                     gTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCa                                       >  2:2571516/1‑93 (MQ=255)
                                                     gTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCa                                       >  1:3117296/1‑93 (MQ=255)
                                                      tCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAg                                      >  1:2137970/1‑93 (MQ=255)
                                                       ccACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAgc                                     >  1:3959967/1‑92 (MQ=255)
                                                       ccACTTCCCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGt                                     <  2:4470824/93‑1 (MQ=255)
                                                        cgcTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTc                                    <  1:2517598/91‑1 (MQ=255)
                                                        cACTCCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTc                                    <  2:4272363/93‑1 (MQ=255)
                                                         cctccGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTcc                                   <  2:4612980/89‑1 (MQ=255)
                                                         cctccGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTcc                                   <  1:2120761/89‑1 (MQ=255)
                                                         aCTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTcc                                   >  1:646784/1‑93 (MQ=255)
                                                          cTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTccc                                  <  1:2434494/93‑1 (MQ=255)
                                                          cTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTccc                                  >  2:1643317/1‑93 (MQ=255)
                                                            tCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGc                                >  1:4285099/1‑93 (MQ=255)
                                                            ccccccACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGc                                <  2:3051509/90‑1 (MQ=255)
                                                              gCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCct                              >  1:4466984/1‑93 (MQ=255)
                                                              gCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCct                              >  2:3355394/1‑93 (MQ=255)
                                                                ccaccaTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctct                            <  1:2190443/93‑1 (MQ=255)
                                                                ccaccaTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctct                            >  1:1894766/1‑93 (MQ=255)
                                                                ccaccaTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctct                            <  2:1067369/93‑1 (MQ=255)
                                                                ccaccaTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCctct                            >  1:3594364/1‑93 (MQ=255)
                                                                 caccaTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCtt                           >  2:872770/1‑93 (MQ=255)
                                                                 caccaTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCtt                           >  1:4084492/1‑93 (MQ=255)
                                                                  cccaAAATTGACGAACCCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCTTTTGACTATCAGTCCCGCCTCTTc                          <  2:4068255/92‑1 (MQ=255)
                                                                  accaTAATTGACGAACGCTCACTGGTGCCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCtt                           >  1:1735672/1‑93 (MQ=255)
                                                                   ccaTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCa                         >  2:4585727/1‑93 (MQ=255)
                                                                    caTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCaa                        >  2:40186/1‑93 (MQ=255)
                                                                    caTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCaa                        <  1:3885496/93‑1 (MQ=255)
                                                                     aTAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAAt                       <  1:126563/93‑1 (MQ=255)
                                                                      tAATTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCTTTTGACTATCAGTCCCGCCTCTTCAATc                      <  1:4052415/93‑1 (MQ=255)
                                                                      tAATTGACGAACGCTCACTGGTG‑‑CCCCCCGGGCTTGCCACGAAACTTAACAGCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCg                     >  1:1994641/1‑93 (MQ=255)
                                                                       aaTTGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCg                     >  2:1251952/1‑93 (MQ=255)
                                                                          tGACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCtt                  <  1:3121341/93‑1 (MQ=255)
                                                                           gACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCATTCCCGCCTCTTCAATCGCTTc                 <  2:4671666/93‑1 (MQ=255)
                                                                           gACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTc                 <  1:334057/93‑1 (MQ=255)
                                                                           gACGAACGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTc                 <  2:648392/93‑1 (MQ=255)
                                                                           gACGAACCCGCATTGGTG‑‑CCCCCCGGGCTTGCCAGGAAACTTAACCCCGGTTTGGTCCGTTTGCCTATCAGTCCCGCCTCTTCAATCGCTTcg                <  2:4495971/93‑1 (MQ=255)
                                                                             cGAACGCCCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTcgc               <  2:451826/93‑1 (MQ=255)
                                                                               aaCGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGa             >  1:2173373/1‑93 (MQ=255)
                                                                               aaCGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGa             >  1:2048026/1‑93 (MQ=255)
                                                                                aCGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGAc            >  2:137252/1‑93 (MQ=255)
                                                                                 cGCTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACg           >  1:2979995/1‑93 (MQ=255)
                                                                                  gcgcACTGGTGCCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACg           <  2:1435237/90‑1 (MQ=255)
                                                                                   cTCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACggg         <  1:4559097/93‑1 (MQ=255)
                                                                                    tCACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGc        >  1:3031573/1‑93 (MQ=255)
                                                                                    tCACTGGGGCCCCCCCGGGCTTGCCAGGAAACTTACCACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGc        <  1:769756/93‑1 (MQ=255)
                                                                                     cACTGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCa       >  1:2273503/1‑93 (MQ=255)
                                                                                     cACTGGTGCCCCCCCGGGCTGGCCAGCAACCTTAACACCGGTTTGCGCCGTTGGACTACCAGGCCCGCCTCTTCAATCGCTTCGCGACGGGCa       <  2:2386781/93‑1 (MQ=255)
                                                                                        tGGTGCCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACa    <  2:3228223/93‑1 (MQ=255)
                                                                                             |                                                                                            
CGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTG‑‑CCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATC  >  NC_000913/3175427‑3175608

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