Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,883,067 T→C D121D (GAT→GAC yidC → 60 KD inner‑membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,883,0670TC98.8% 511.5 / ‑2.6 166D121D (GAT→GACyidC60 KD inner‑membrane protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (86/78);  total (88/78)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.58e-01

Only 100 of 177 total aligned reads displayed.

AACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGATGA  >  NC_000913/3882978‑3883158
                                                                                         |                                                                                           
 aCTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAAcc                                                                                         >  2:10024/1‑93 (MQ=255)
   tctACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCg                                                                                       >  1:3343671/1‑93 (MQ=255)
       cccAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTa                                                                                   >  1:2072863/1‑93 (MQ=255)
        ccAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTGGTGATGGCCCGGACAACCCGGCTaa                                                                                  >  1:1628927/1‑93 (MQ=255)
         cAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAAc                                                                                 <  1:645487/93‑1 (MQ=255)
           gCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACgg                                                                               >  1:3267747/1‑93 (MQ=255)
            ccGGTCCAGCTGTTAGAAACGTCACCGCAGTTTACTTATCAGACCCAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGc                                                                              <  2:3969109/93‑1 (MQ=255)
              gTTCCAGCTGTTGGAAACTTCACCGCAGTTTCTTTATCAGGCACCGAGCGGTCTCACCAGTCGTGATGGCCCGGACCACCCGGCTAACCGccc                                                                            >  1:1255397/1‑93 (MQ=255)
              gTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGccc                                                                            >  1:4619304/1‑93 (MQ=255)
              gTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGccc                                                                            >  1:155431/1‑93 (MQ=255)
              gTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGccc                                                                            >  2:2485559/1‑93 (MQ=255)
               ttCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCcg                                                                           >  2:2986580/1‑93 (MQ=255)
               ttCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCcg                                                                           >  2:4114682/1‑93 (MQ=255)
               ttCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCcg                                                                           >  1:2641651/1‑93 (MQ=255)
                 ccAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCcgcg                                                                         >  2:587569/1‑93 (MQ=255)
                  cAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGTGCGTGATGGCCCGGACAACCCGGCTACCGGCCcgcgc                                                                        >  1:2595131/1‑92 (MQ=255)
                  cAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGt                                                                        >  1:82436/1‑93 (MQ=255)
                  cAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGCCCCGGACAACCCGGCTAACGTCCCGCGt                                                                        >  1:1778018/1‑93 (MQ=255)
                   aGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTc                                                                       >  2:2843003/1‑93 (MQ=255)
                    gCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTcc                                                                      <  1:2025805/93‑1 (MQ=255)
                      tGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGc                                                                    >  2:1778625/1‑93 (MQ=255)
                       gTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCt                                                                   >  1:3520199/1‑93 (MQ=255)
                        ttGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTg                                                                  >  2:1795479/1‑93 (MQ=255)
                        ttGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTg                                                                  >  2:1093949/1‑93 (MQ=255)
                          ggAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGta                                                                <  2:4354504/93‑1 (MQ=255)
                             aaCTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATaa                                                             >  2:1576846/1‑93 (MQ=255)
                             aaCTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATaa                                                             >  2:337132/1‑93 (MQ=255)
                              aCTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAAc                                                            >  2:4256067/1‑93 (MQ=255)
                               cTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACg                                                           >  1:1442346/1‑93 (MQ=255)
                               cTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACg                                                           <  1:2330568/93‑1 (MQ=255)
                                  cACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTg                                                        >  2:1145599/1‑93 (MQ=255)
                                   cccGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGa                                                       <  2:1191449/92‑1 (MQ=255)
                                    ccGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATCACGTTGaa                                                      >  2:3505084/1‑93 (MQ=255)
                                     cGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGCTGGGCCGGACAACCCGGCTAACGACCCACGCCCCCTGTTCAACGTTCaac                                                     >  1:4169428/1‑92 (MQ=255)
                                      gcggTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaaaa                                                    <  2:54766/90‑1 (MQ=255)
                                      gCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaaaa                                                    >  2:1371511/1‑93 (MQ=255)
                                      gCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaaaa                                                    >  1:3034849/1‑93 (MQ=255)
                                       cAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaaaaa                                                   <  1:1767232/93‑1 (MQ=255)
                                        aGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAg                                                  <  2:849232/93‑1 (MQ=255)
                                        aGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGACCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAg                                                  <  1:2347376/93‑1 (MQ=255)
                                         gTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGa                                                 <  1:1387996/93‑1 (MQ=255)
                                         gTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGa                                                 <  1:2833963/93‑1 (MQ=255)
                                         gTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGa                                                 <  1:4028034/93‑1 (MQ=255)
                                          tttatttaTCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGAc                                                <  2:1739363/93‑1 (MQ=255)
                                          tttatttaTCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGAc                                                >  1:274388/1‑93 (MQ=255)
                                            tatttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGc                                              <  2:2333412/93‑1 (MQ=255)
                                            tatttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGc                                              <  2:4553670/93‑1 (MQ=255)
                                             atttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCt                                             >  2:1480931/1‑93 (MQ=255)
                                             atttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAAACCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCt                                             <  2:1878297/93‑1 (MQ=255)
                                              tttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCtt                                            <  1:1664645/93‑1 (MQ=255)
                                               ttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTa                                           <  1:2477943/93‑1 (MQ=255)
                                                tatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTAt                                          <  1:4679425/93‑1 (MQ=255)
                                                 atCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTAtg                                         <  2:3557778/93‑1 (MQ=255)
                                                   cAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTAtgtg                                       <  1:4448040/93‑1 (MQ=255)
                                                    aGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgc                                      <  2:2628094/93‑1 (MQ=255)
                                                    aGGCACAGAGCGGGCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCCGGAGAACGGTGAAAAAGACGCTCACGTgc                                      >  1:3601963/1‑93 (MQ=255)
                                                     ggCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgct                                     >  2:1202956/1‑93 (MQ=255)
                                                     ggCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgct                                     <  2:4561507/93‑1 (MQ=255)
                                                     ggCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgct                                     >  2:684675/1‑93 (MQ=255)
                                                      gCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctg                                    >  1:3067994/1‑93 (MQ=255)
                                                      gCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctg                                    <  2:127566/93‑1 (MQ=255)
                                                      gCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctg                                    <  2:127565/93‑1 (MQ=255)
                                                      gCACAGAGCGGTCTGACCGGTCATGATGGCACGGAAAACCCTGCCCACGGCTCGCGTCCGCTGTATAATGTTGAAAATGGCGCATATGCgctg                                    >  2:4429399/1‑93 (MQ=255)
                                                        acaGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctggc                                  >  2:229995/1‑93 (MQ=255)
                                                        acaGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctggc                                  <  1:708017/93‑1 (MQ=255)
                                                         caGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctggct                                 >  1:1970854/1‑93 (MQ=255)
                                                           gagCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGa                               >  1:4192913/1‑93 (MQ=255)
                                                             gCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAg                             <  2:2315841/93‑1 (MQ=255)
                                                             gCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCGAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAg                             <  1:545782/93‑1 (MQ=255)
                                                               ggTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGt                           <  2:4722135/93‑1 (MQ=255)
                                                               ggTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGt                           <  1:632104/93‑1 (MQ=255)
                                                                gTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTc                          <  2:3855082/93‑1 (MQ=255)
                                                                  cTGACCGGTCGTGATGGCCCGGATAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAACAAGACGCTTATGTGCTGGCTAAAGGTaaa                        >  1:4036644/1‑93 (MQ=255)
                                                                   tGACCGGTCGTGAAGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCaaa                       >  2:887069/1‑93 (MQ=255)
                                                                    gACCGGTCGTGATGACCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCaaaa                      >  1:2100018/1‑93 (MQ=255)
                                                                      ccGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACg                    >  1:1895002/1‑93 (MQ=255)
                                                                       cGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGATGAAAAAGACACTTATGTGCTGGCTGAATGTTGAAAcca                   >  2:160183/1‑91 (MQ=255)
                                                                        ggTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGaa                  >  2:3272074/1‑93 (MQ=255)
                                                                        ggTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGaa                  >  2:3288933/1‑93 (MQ=255)
                                                                        ggTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGaa                  <  1:2095842/93‑1 (MQ=255)
                                                                        ggTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGaa                  >  2:4101986/1‑93 (MQ=255)
                                                                         gTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAAc                 >  1:1157212/1‑93 (MQ=255)
                                                                         gTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAAc                 <  1:1196812/93‑1 (MQ=255)
                                                                           cGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTg               >  1:3296936/1‑93 (MQ=255)
                                                                            gTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCCCTGTATAACGTTGAAAAAGACGCTTATGTGCTCGCTGAAGGTCAAAACGAACTGc              >  1:3766422/1‑93 (MQ=255)
                                                                             tGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCa             <  2:1968334/93‑1 (MQ=255)
                                                                             tGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCa             <  2:4011689/93‑1 (MQ=255)
                                                                             tGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCa             >  2:3225273/1‑93 (MQ=255)
                                                                              gATGGCCCGTACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAg            <  2:2533362/93‑1 (MQ=255)
                                                                              gATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAg            <  2:3341805/93‑1 (MQ=255)
                                                                              gATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAg            >  1:4112046/1‑93 (MQ=255)
                                                                              gATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAg            <  2:4674159/93‑1 (MQ=255)
                                                                              gATGGCCCGGACAACCCGGCTAACGGCACGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAg            <  2:2533361/93‑1 (MQ=255)
                                                                               aTGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAgg           >  2:3102864/1‑93 (MQ=255)
                                                                               aTGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAgg           >  1:4722370/1‑93 (MQ=255)
                                                                                 ggCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTg         >  2:1632844/1‑93 (MQ=255)
                                                                                    ccGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCg      <  1:75130/93‑1 (MQ=255)
                                                                                     cGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGa     <  2:2390527/93‑1 (MQ=255)
                                                                                     cGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGa     >  2:319647/1‑93 (MQ=255)
                                                                                        aTAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGATGa  >  1:3136148/1‑93 (MQ=255)
                                                                                         |                                                                                           
AACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGATGA  >  NC_000913/3882978‑3883158

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