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,0670TC99.8% 1305.2 / ‑5.3 410D121D (GAT→GACyidC60 KD inner‑membrane protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (234/175);  total (235/175)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.27e-01

Only 100 of 429 total aligned reads displayed.

GAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGAT  >  NC_000913/3882977‑3883156
                                                                                          |                                                                                         
 aaCTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAAc                                                                                        >  1:774078/1‑93 (MQ=255)
    tctACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCg                                                                                     >  2:9711280/1‑93 (MQ=255)
    tctACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCg                                                                                     >  1:6199162/1‑93 (MQ=255)
    tctACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCg                                                                                     >  2:8432738/1‑93 (MQ=255)
     ctACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCgg                                                                                    <  2:8019416/93‑1 (MQ=255)
        cccAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTa                                                                                 >  2:3086246/1‑93 (MQ=255)
         ccAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTaa                                                                                <  2:3589166/93‑1 (MQ=255)
          cAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAAc                                                                               <  1:9393571/93‑1 (MQ=255)
          cAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAAc                                                                               <  1:757395/93‑1 (MQ=255)
          cAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAAc                                                                               >  1:5867340/1‑93 (MQ=255)
          cAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGAGGTCTGACCGGTCGTGATGGCCCGCACAACCCGGCTAAc                                                                               >  1:3359190/1‑93 (MQ=255)
            gCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACgg                                                                             >  1:9837784/1‑93 (MQ=255)
              cGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGcc                                                                           <  1:2959218/93‑1 (MQ=255)
               gTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGAGAACCCGGCTAACGGccc                                                                          >  1:8130387/1‑93 (MQ=255)
                ttCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGcccc                                                                         >  2:7284197/1‑92 (MQ=255)
                ttCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCcg                                                                         <  1:6921964/93‑1 (MQ=255)
                ttCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCcg                                                                         >  2:5665820/1‑93 (MQ=255)
                   cAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGt                                                                      >  2:1886641/1‑93 (MQ=255)
                   cAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGt                                                                      <  2:819534/93‑1 (MQ=255)
                    cgCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTc                                                                     <  2:8343010/92‑1 (MQ=255)
                    aGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTc                                                                     >  1:2887882/1‑93 (MQ=255)
                       tGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGc                                                                  <  2:9061260/93‑1 (MQ=255)
                        gTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCt                                                                 >  1:8003852/1‑93 (MQ=255)
                          tGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGt                                                               >  2:1182862/1‑93 (MQ=255)
                          tGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAAACCGGCTAACGGCCCGCGTCCGCTGt                                                               >  2:3750481/1‑93 (MQ=255)
                           ggAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGta                                                              >  2:1012709/1‑93 (MQ=255)
                           ggAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGta                                                              <  1:4244461/93‑1 (MQ=255)
                           ggAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGta                                                              >  2:2058682/1‑93 (MQ=255)
                           ggAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGta                                                              <  1:2504390/93‑1 (MQ=255)
                              aaCTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATaa                                                           <  1:3232055/93‑1 (MQ=255)
                              aaCTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATaa                                                           >  1:886616/1‑93 (MQ=255)
                               aCTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAAc                                                          <  1:307866/93‑1 (MQ=255)
                                cTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACg                                                         <  1:1923734/93‑1 (MQ=255)
                                cTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACg                                                         >  2:1769493/1‑93 (MQ=255)
                                 ttCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGt                                                        <  1:6654363/93‑1 (MQ=255)
                                 ttCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGt                                                        <  1:1385626/93‑1 (MQ=255)
                                  tCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGtt                                                       <  2:9533088/93‑1 (MQ=255)
                                  tCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGtt                                                       >  2:6343262/1‑93 (MQ=255)
                                  tCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGtt                                                       <  1:5681693/93‑1 (MQ=255)
                                  tCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGtt                                                       <  2:7728766/93‑1 (MQ=255)
                                    aCCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCCTCCGCTGTATAACGTTGa                                                     >  2:4171056/1‑93 (MQ=255)
                                     ccGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaa                                                    >  2:470059/1‑93 (MQ=255)
                                     ccGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaa                                                    >  2:2719099/1‑93 (MQ=255)
                                      cGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaaa                                                   <  2:5938668/93‑1 (MQ=255)
                                      cGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaaa                                                   >  1:3626345/1‑93 (MQ=255)
                                        cAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGaaaaa                                                 >  2:2162375/1‑93 (MQ=255)
                                         aGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAg                                                >  1:1127256/1‑93 (MQ=255)
                                         aGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAg                                                <  1:4724264/93‑1 (MQ=255)
                                          gTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGa                                               >  1:10014663/1‑93 (MQ=255)
                                          gTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGa                                               <  2:7588110/93‑1 (MQ=255)
                                          gTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGa                                               >  2:3402445/1‑93 (MQ=255)
                                          gTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGa                                               <  2:784294/93‑1 (MQ=255)
                                           tttatttaTCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGAc                                              >  1:8658642/1‑93 (MQ=255)
                                           tttatttaTCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGAc                                              >  2:5037743/1‑93 (MQ=255)
                                           tttatttaTCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGAc                                              <  1:381121/93‑1 (MQ=255)
                                           tttatttaTCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGAc                                              >  2:2926057/1‑93 (MQ=255)
                                             tatttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACgt                                            >  2:4331058/1‑92 (MQ=255)
                                             tatttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGc                                            >  2:3843194/1‑93 (MQ=255)
                                             tatttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGc                                            >  2:7104226/1‑93 (MQ=255)
                                             tatttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGc                                            <  1:5164656/93‑1 (MQ=255)
                                              atttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCt                                           >  1:7930687/1‑93 (MQ=255)
                                              atttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCt                                           >  1:4527944/1‑93 (MQ=255)
                                               tttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCtt                                          >  2:113231/1‑93 (MQ=255)
                                                ttatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTa                                         <  1:4076638/93‑1 (MQ=255)
                                                 tatCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTAt                                        <  1:4353747/93‑1 (MQ=255)
                                                  atCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTAtg                                       <  2:8666206/93‑1 (MQ=255)
                                                    cAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTAtgtg                                     >  2:2325187/1‑93 (MQ=255)
                                                     aGGCACAGAGCTGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATCACGTTGAAAAAGACGCTTATGTgc                                    >  2:6770239/1‑93 (MQ=255)
                                                     aGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgc                                    <  1:7860590/93‑1 (MQ=255)
                                                      ggCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgct                                   <  2:6831911/93‑1 (MQ=255)
                                                       gCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctg                                  <  1:2778295/93‑1 (MQ=255)
                                                         acaGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTTTGTgctggg                                >  2:7395728/1‑92 (MQ=255)
                                                          caGAGCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTTTATAACGTTGAAAAAGACGCTTATGTgctggct                               >  1:9685009/1‑93 (MQ=255)
                                                           agagCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctggctg                              >  2:982437/1‑93 (MQ=255)
                                                           agagCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctggctg                              >  1:1617455/1‑93 (MQ=255)
                                                           agagCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTgctggctg                              >  2:5858905/1‑93 (MQ=255)
                                                            gagCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGa                             >  2:6154932/1‑93 (MQ=255)
                                                            gagCGGTCTGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGa                             >  1:9021898/1‑93 (MQ=255)
                                                                    tGACCGGTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGGAGGTCaaa                     >  2:8876247/1‑93 (MQ=255)
                                                                         ggTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGaa                <  2:5473011/93‑1 (MQ=255)
                                                                         ggTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGaa                >  1:1485845/1‑93 (MQ=255)
                                                                         ggTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGaa                <  1:1196159/93‑1 (MQ=255)
                                                                         ggTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGaa                >  1:8783058/1‑93 (MQ=255)
                                                                          gTCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAAAGAAc               >  2:3061594/1‑93 (MQ=255)
                                                                           tCGTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACt              >  1:1302163/1‑93 (MQ=255)
                                                                             gTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTTCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGGCAAAAAGAACTGc            >  2:1420008/1‑93 (MQ=255)
                                                                             gTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAATGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGc            <  2:6005907/93‑1 (MQ=255)
                                                                             gTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGc            <  1:8475864/93‑1 (MQ=255)
                                                                             gTGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGc            >  2:5271247/1‑93 (MQ=255)
                                                                              tGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCa           >  1:2340025/1‑93 (MQ=255)
                                                                              tGATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCa           >  1:2227362/1‑93 (MQ=255)
                                                                               gATGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAg          <  1:6298730/93‑1 (MQ=255)
                                                                                aTGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAgg         <  2:286492/93‑1 (MQ=255)
                                                                                 tGGCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGt        <  1:7142801/93‑1 (MQ=255)
                                                                                  ggCCCGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTg       <  1:6119727/93‑1 (MQ=255)
                                                                                     ccGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGGGCCg    >  2:2566340/1‑93 (MQ=255)
                                                                                      cGGATAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGa   >  2:9208178/1‑93 (MQ=255)
                                                                                      cGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGa   >  1:3185932/1‑93 (MQ=255)
                                                                                      cGGACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGa   >  1:3918919/1‑93 (MQ=255)
                                                                                       ggACAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGAt  <  2:1661811/93‑1 (MQ=255)
                                                                                          |                                                                                         
GAACTCTACCCAGCCGTTCCAGCTGTTGGAAACTTCACCGCAGTTTATTTATCAGGCACAGAGCGGTCTGACCGGTCGTGATGGCCCGGATAACCCGGCTAACGGCCCGCGTCCGCTGTATAACGTTGAAAAAGACGCTTATGTGCTGGCTGAAGGTCAAAACGAACTGCAGGTGCCGAT  >  NC_000913/3882977‑3883156

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