Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,326,072 T→G 8.8% intergenic (‑331/+327) yjdN ← / ← yjdM metalloprotein superfamily protein/zinc‑ribbon family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,326,0720TG 8.8% 119.5 / 2.4 60intergenic (‑331/+327)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base T (28/26);  new base G (4/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCG  >  NC_000913/4325924‑4326220
                                                                                                                                                    |                                                                                                                                                    
gCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGGCGGCACCGTTGCCGGATGCGGCGTAAACGCCt                                                                                                                                                      >  2:46863/1‑149 (MQ=37)
     gCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGTAATTTGTAGGCCTGATAAGACGCGCAAGCGTCTCATCAGGCAGTCGGCACCGTTGCCGTATGCGGCGTAAACGCCTTATcc                                                                                                                                                 <  1:353943/149‑1 (MQ=255)
      cATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                                                                                                >  1:920187/1‑149 (MQ=37)
       aTCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGGGTAAACGCCTTATCCgg                                                                                                                                               >  2:505794/1‑149 (MQ=37)
        tCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                              <  2:920187/149‑1 (MQ=37)
         cAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                             <  2:642702/149‑1 (MQ=37)
           ggCAGTCGGCACCATAGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                                                                           <  2:225080/149‑1 (MQ=255)
            gCAGTCGGCACTATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTTTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTGGCCGGATAAGGCGTTTACGCCGCATCCGGcaacg                                                                                                                                          <  2:215491/149‑5 (MQ=11)
            gCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCCTCAGGCAGTCGGCACCGTTGACGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                                          >  1:17550/1‑149 (MQ=255)
                tCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAAGGAGGGTGCGGGATTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAt                                                                                                                                      <  2:612953/149‑1 (MQ=255)
                   gCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCGTACGGAGGGTGCGGGTATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGCCAGTCGGCACCCTTGCCGGATGCGGCGTAAACGCCTTATCCGTCCTAGGGATGGt                                                                                                                                   >  1:409205/1‑149 (MQ=11)
                             cGGATTCTGCGTAAACGCCTTATCCGCCCTACGGAGGGTGCTGGAAGTTGTAGGCCTGATAAGACGCGCAATCGTCGCATCAGGCAGTCGGCACCGGTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACTGATGGTGCGGGAAttt                                                                                                                         <  1:323559/149‑1 (MQ=11)
                              ggATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCTTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTg                                                                                                                        <  2:729024/149‑1 (MQ=17)
                                 tGCGGCGTAACCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAg                                                                                                                     <  1:384978/149‑1 (MQ=17)
                                 tGCGGCGTAAACGCCTTATCCGGCCTACGGAGTGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAg                                                                                                                     <  2:264117/149‑1 (MQ=17)
                                    ggcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGAGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGcc                                                                                                                  <  2:442069/149‑1 (MQ=17)
                                     gcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCt                                                                                                                 <  2:422521/149‑1 (MQ=32)
                                       gTAAACGGCTATTCCGGCCTACGTAGGGTGCGGGAATTTGTAGGCCTGATACGACGCGCAAGCGTGGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGa                                                                                                               <  1:300502/149‑1 (MQ=11)
                                       gTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGa                                                                                                               >  1:778213/1‑149 (MQ=32)
                                        tAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGAt                                                                                                              <  1:983484/149‑1 (MQ=32)
                                         aaaCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATa                                                                                                             <  2:634650/149‑1 (MQ=32)
                                               cTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGGGGGAATTTGTAGGCGTGATAAGAcgc                                                                                                       >  1:594859/1‑149 (MQ=17)
                                                ttATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTTTAGGCCTGATAAGAcgcg                                                                                                      >  2:284036/1‑149 (MQ=17)
                                                     cGGCCGACTCAGGGTGCTGTAATTAGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGGCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTAGTAGGCCTGATAAGACGCGCAAGc                                                                                                 <  1:694251/149‑1 (MQ=11)
                                                        ccTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATTGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                                                                                              <  2:962811/149‑1 (MQ=17)
                                                             ggATGGTGCGGGATTCTGTAGGCGTGATAAGTCGCGCAAGCGTCGCATCAGGCAGGCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                                                                         <  2:592322/149‑1 (MQ=11)
                                                                gggTGCGGGAATTTGTAGGCCGGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAgg                                                                                      <  2:478299/149‑1 (MQ=17)
                                                                  gTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCa                                                                                    <  2:854960/149‑1 (MQ=32)
                                                                      gggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCCGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATATGTAGGCCTGATAAGACGCGCAAGCGTCGCAGCAGGCAGGTg                                                                                >  1:637284/1‑149 (MQ=255)
                                                                       ggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGGGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGCTgg                                                                               >  1:688433/1‑149 (MQ=17)
                                                                            ttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGAAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTg                                                                          >  1:804945/1‑149 (MQ=38)
                                                                                  gCCTGATAAGACGCGCAAGCGTCGCATCAGGGAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCg                                                                    <  1:235842/149‑1 (MQ=255)
                                                                                  gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACAGTTGCCg                                                                    >  2:173375/1‑149 (MQ=37)
                                                                                    cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGa                                                                  >  1:599600/1‑149 (MQ=39)
                                                                                    cTGATAAGACGCGCAAGCGTAGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATACGGCCTACGTATGGTGCCGGAATTTGTAGGCCTGATAATACTCGAAAGCGTCGCATAAGGCAGATGGCACAGTTTCCGGa                                                                  >  2:729173/1‑149 (MQ=11)
                                                                                        tAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCGGGCAGTTGGCACGGTTGCCGGATgcg                                                              >  2:275267/1‑149 (MQ=255)
                                                                                           gACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCATTTGGCACAGTTGCCGGATgcggcg                                                           >  1:606804/1‑149 (MQ=17)
                                                                                            aCGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGt                                                          >  1:1009189/1‑149 (MQ=21)
                                                                                              gcgcAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGGGAGGGAAGTTGTAGGCCAGATAAGACGCGCAAGCGTCGCATCAGGCAGGAGGAACTGTTGCAGGATgcggcggta                                                        >  2:226849/1‑146 (MQ=11)
                                                                                                   aGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGcc                                                   <  2:778213/149‑1 (MQ=39)
                                                                                                      gTCGCATCAGGCAGTCGGCACGGTTGCGGGATGCGGCGTAAACGCCGTATCCGGCCTACGGATGGTGCGGGTATTTGTTGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTa                                                <  2:600365/149‑1 (MQ=14)
                                                                                                      gTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTa                                                <  1:7930/149‑1 (MQ=39)
                                                                                                          cATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                            >  2:25278/1‑149 (MQ=39)
                                                                                                          cATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGGTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                            >  2:295684/1‑149 (MQ=255)
                                                                                                                 cAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACg                                     <  2:804945/149‑1 (MQ=39)
                                                                                                                      ggCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAggg                                >  1:524846/1‑149 (MQ=39)
                                                                                                                      ggCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGAATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGAAAACGCCTTAACCGGCCTACGGAggg                                >  2:460668/1‑149 (MQ=255)
                                                                                                                      ggCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGGGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCGGCGGAggg                                >  1:489017/1‑149 (MQ=255)
                                                                                                                       gCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGc                                                                              <  1:971876/102‑1 (MQ=17)
                                                                                                                               gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAt                       >  1:185403/1‑149 (MQ=11)
                                                                                                                               gCCGGATAAGGCGTTTACGCCGCATCCCGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCCCAGTTGCCGCATGCGGCGTAAACGCCATATCCGGCCTACGGAGGGTGCTGGAAt                       >  2:255854/1‑149 (MQ=11)
                                                                                                                                 cGGATAAGGCGTTTACGCCGCATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAgtt                     >  1:110507/1‑149 (MQ=11)
                                                                                                                                   gataaggcgTTTACGCCGCATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGt                   >  2:659985/6‑149 (MQ=11)
                                                                                                                                         gcgTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATcc                                             <  1:150847/117‑1 (MQ=39)
                                                                                                                                          cgTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTg            >  1:140742/1‑149 (MQ=37)
                                                                                                                                              aaCGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATaa        >  2:275793/1‑149 (MQ=37)
                                                                                                                                                 gCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCTGGAATTTGTAGGCCTGATAAGAc     >  1:197573/1‑149 (MQ=37)
                                                                                                                                                 gCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAc     >  1:1016254/1‑149 (MQ=37)
                                                                                                                                                 gCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAc     >  1:1012161/1‑149 (MQ=37)
                                                                                                                                                  ccTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGGTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGAGCGGGAAGTTGTAGGCCTGCTAAGAcg    >  2:952575/1‑149 (MQ=255)
                                                                                                                                                   cTTATCCGGCCTACGGCTGGTGCGGGAATTTGTAGGCCTGATATGACGCGCAAGCGTCGCATCAGGCAGTTGGCTCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGTAATTTGTAGGCCTGATAAGAcgc   <  2:409205/149‑1 (MQ=255)
                                                                                                                                                    ttATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgcg  >  2:690170/1‑149 (MQ=37)
                                                                                                                                                    |                                                                                                                                                    
GCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTTGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCG  >  NC_000913/4325924‑4326220

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