Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,325,948 A→G 6.6% intergenic (‑207/+451) 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,325,9480AG 6.6% 208.2 / 4.9 90intergenic (‑207/+451)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base A (36/47);  new base G (4/2);  total (40/50)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

GAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGG  >  NC_000913/4325800‑4326096
                                                                                                                                                    |                                                                                                                                                    
gagGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCa                                                                                                                                                      >  1:634470/1‑149 (MQ=39)
  ggcggcCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCAtt                                                                                                                                                    <  2:1001881/149‑1 (MQ=39)
   gcggcCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTg                                                                                                                                                   >  1:543493/1‑149 (MQ=38)
   gcggcCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCGACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTg                                                                                                                                                   >  1:66047/1‑149 (MQ=38)
    cggcCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGc                                                                                                                                                  <  1:243217/149‑1 (MQ=11)
     ggcCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGcc                                                                                                                                                 <  1:41651/149‑1 (MQ=17)
      gcCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCg                                                                                                                                                >  2:949672/1‑149 (MQ=17)
       ccTGATGCGACGCTCGCGCGTCTTATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGGAAACGCCTGATCCGGCCTACGGAGGGGGCGGGACTTGGAAGGCCTGATAAGCCGCGCAAGCGTCGCATGAGGCAGTCGGCACCGTTGCCgg                                                                                                                                               >  2:768988/1‑149 (MQ=11)
       ccTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCgg                                                                                                                                               <  2:54594/149‑1 (MQ=37)
        cTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGa                                                                                                                                              <  2:249185/149‑1 (MQ=37)
         tGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGAt                                                                                                                                             >  2:466891/1‑149 (MQ=37)
            tAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCAGAAAACACCCGCAAGCCTCCAAACAGTCACAACGCACCATAGCCGGAAgcg                                                                                                                                          >  1:819767/1‑149 (MQ=11)
               gACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcggcg                                                                                                                                       >  1:1064622/1‑149 (MQ=37)
                 cgcgcgAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTa                                                                                                                                     <  1:378435/149‑1 (MQ=37)
                  gcgcgAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTaa                                                                                                                                    >  1:304600/1‑149 (MQ=17)
                      gAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATAAGGCGTTTACGc                                                                                                                                <  2:634470/149‑1 (MQ=255)
                       aGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCGTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGcc                                                                                                                               <  1:704979/149‑1 (MQ=21)
                         cGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATAGCCGGATGCGGCGTAAACGCCtt                                                                                                                             <  2:739736/149‑1 (MQ=21)
                            gtcATCAGGCAGTCGGCTCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATc                                                                                                                          <  2:1066369/147‑1 (MQ=21)
                                tCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                      <  2:1110071/149‑1 (MQ=34)
                                    gCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                                                                  >  2:901002/1‑149 (MQ=34)
                                           gCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGt                                                                                                           >  2:588726/1‑149 (MQ=32)
                                             aCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGc                                                                                                         <  2:741274/149‑1 (MQ=32)
                                              cTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                        <  2:318239/149‑1 (MQ=32)
                                              cTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                        <  1:901002/149‑1 (MQ=32)
                                              cTGTGGCCGGAGGCGGCGTAACCGCCGTCTCCGGCCTACGGAGGGGGCGGGAATCTGTAGGCCTGCTCAGACGCGCAAGCGTCGCATCAGGCAGTCGGCCCCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                        <  1:768988/149‑1 (MQ=11)
                                               tGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCgg                                                                                                       <  2:543493/149‑1 (MQ=32)
                                               tGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCgg                                                                                                       <  2:1141686/149‑1 (MQ=32)
                                                   gCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAt                                                                                                   >  2:571862/1‑149 (MQ=32)
                                                    ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAtt                                                                                                  >  2:69333/1‑149 (MQ=32)
                                                       gATGCGGCGTAACCGCCTTATCCGGCCTACGGAGTGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTTCGGAGGGTGCGGGAATTTGt                                                                                               <  2:211151/149‑1 (MQ=12)
                                                       gATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATACGACGCGCAAGCGTCGCAGCAGGCAGTCGGCACCATTGCCGGATGCGGCGTCTCCGCCTTATCCGGGCGCCGGCGGGTGCGGGCATATGt                                                                                               >  1:789182/1‑149 (MQ=11)
                                                        aTGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTa                                                                                              <  1:258345/149‑1 (MQ=32)
                                                          tcggcgTAAACGCCTTATCCGTCATACGTAGGGTCCGGGAATTTGCATGCCGGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGAGGCGTAAACGCCTTATCCTGCCTACGGAGGGTGCGGGAATTTGTAgg                                                                                            <  1:1066127/148‑1 (MQ=11)
                                                               gTAAACGCCTTATCCGGCCTACGGAGGGTGCGGTAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGa                                                                                       <  2:549300/149‑1 (MQ=17)
                                                                      ccttcgccGGCCTCCGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGGCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcg                                                                                <  1:610555/143‑1 (MQ=255)
                                                                       cTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                               >  1:748735/1‑149 (MQ=32)
                                                                       cTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTAGGCACAATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGTCGTGTGCGGTAATTTGTAGGCCTTCTAAGAcgc                                                                               >  2:857064/1‑149 (MQ=11)
                                                                       cTTATCCGGCCTACGGAGGGAGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGAAGTAGGCACCATTGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACGTAGGGTGCGGGATTGTGTGGGGGTGATAAGAcgc                                                                               >  2:517882/1‑149 (MQ=11)
                                                                             cGGCCTACGTAGGGTCCGGGAATTTGTAGGCCTGATACGACGCGCTAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACTGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGc                                                                         <  2:72373/149‑1 (MQ=11)
                                                                                ccgaCGGAGGGCGCGGGCAGTTGTAGGCCTGATAAGACGCGCTAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                                                                      <  2:782796/146‑1 (MQ=11)
                                                                                ccTACGTAGGGTGCGGTAAATTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTc                                                                      <  2:1037213/149‑1 (MQ=12)
                                                                                ccTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCATGATAATACGCGCAAGCGtg                                                                      >  2:160597/1‑148 (MQ=17)
                                                                                     ggAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                                                 <  1:1030679/149‑1 (MQ=32)
                                                                                        gggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAgg                                                              <  1:519292/149‑1 (MQ=32)
                                                                                        gggTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAgg                                                              <  1:466891/149‑1 (MQ=32)
                                                                                          ctGCGGGAATTGGCAGCCCTGAGACGACGCCCAAGCGCCGCATCAGGCAGGCGGCACCATAGCCGGAAGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATGCGACGCTTGCGCGTCTTATCAGGcc                                                            <  1:1100729/148‑2 (MQ=11)
                                                                                              gggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGGGGGAATTTGTAGGCCTGATAAGACGCGAAAGCGTCGCATCAGGCAGTCg                                                        >  2:132465/1‑149 (MQ=25)
                                                                                               ggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGTTGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCgg                                                       >  1:1123105/1‑149 (MQ=14)
                                                                                                  aTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCAc                                                    >  1:970963/1‑149 (MQ=37)
                                                                                                   tttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCAcc                                                   <  2:66047/149‑1 (MQ=37)
                                                                                                      gTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTACCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGGAGGCCTGATTAGACGCGCTAGCGGCGGATCAGGCAGTCGGCGCCGtt                                                >  1:1096818/1‑149 (MQ=11)
                                                                                                       gaGGCCTGATCAGCCGCGCAAGCGTCGCATCAGGCAGTCGGCACCAGTGCCGGATGCGGCGTAAACGCCTTAGCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTg                                               <  1:5352/148‑1 (MQ=255)
                                                                                                        aGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGc                                              >  2:546722/1‑149 (MQ=37)
                                                                                                         ggCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCTGTCGGCACCGTTGcc                                             >  1:115514/1‑149 (MQ=255)
                                                                                                          gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCg                                            <  2:375227/149‑1 (MQ=37)
                                                                                                          gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCg                                            <  1:571862/149‑1 (MQ=37)
                                                                                                          gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCg                                            >  2:999605/1‑149 (MQ=37)
                                                                                                               aTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgc                                       >  2:423935/1‑149 (MQ=37)
                                                                                                               aTAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgc                                       >  2:418796/1‑149 (MQ=37)
                                                                                                                tacgACGCGCAAGCGTCGCATCAGGCAGTCGTCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcg                                      <  2:879811/146‑1 (MQ=255)
                                                                                                                 aaGACGCGCAACCTTCGCATCATGCAGTCGGCACCATTGGCGGATGCGGCGTAAACGCATTATCCGGCCTACGGAGGGTGCTGTAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcgg                                     <  2:185006/149‑1 (MQ=11)
                                                                                                                  ataCGCGGAAGCGTCGCATGAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGGCTTTTCCGTCATACTGAGGGTGTGTGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcggc                                    <  2:56776/147‑1 (MQ=11)
                                                                                                                  aGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcggc                                    <  2:422632/149‑1 (MQ=37)
                                                                                                                   gACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                                                    >  2:296736/1‑68 (MQ=35)
                                                                                                                   gACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcggcg                                   >  2:1143711/1‑149 (MQ=37)
                                                                                                                   gACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGACGGATGCGGCGTAAACGCCTTATCCGGCCGACGGAGGGTGCGGGAATATGTAGGCCTGATAAGACGCGCAAGCGTAGAATCAGGCAGTAGGCAACGTTGCCGGATgcggcg                                   >  1:951189/1‑149 (MQ=11)
                                                                                                                        gcAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATTAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAAc                              >  1:353599/1‑149 (MQ=37)
                                                                                                                        gcAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACgcgc                                                                             >  2:768124/1‑102 (MQ=32)
                                                                                                                            gCGTCGCATCAGGCAGTCGGCACTATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCt                          <  2:338683/149‑1 (MQ=17)
                                                                                                                               tCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCaa                                                                           >  1:770643/1‑97 (MQ=33)
                                                                                                                               tCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCaa                                                                           <  2:770643/97‑1 (MQ=33)
                                                                                                                                  cATCAGGCAGTCGGCACCTTTGCCGGTTGCGGCGTACACGCCTTATCCGGCCGACGGAGGGGGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCg                    <  1:221280/149‑1 (MQ=11)
                                                                                                                                  cATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTg                                                                                                >  2:1009476/1‑73 (MQ=34)
                                                                                                                                   aTCAGGCATTCGGCAACAGTGCCGGAAGCGGCGTAAACGCCGTATCCTGCCTATGCAGGGTGCGGGAATTTGTCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGGCGGCTCCGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                   <  2:789182/149‑1 (MQ=11)
                                                                                                                                   aTCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTAGTAGGCCTGATAAGACGCGCAAGCGTCGCAGCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTGATCCgg                   >  1:569223/1‑149 (MQ=255)
                                                                                                                                    tCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGcc                                                                                          <  2:854681/77‑1 (MQ=34)
                                                                                                                                    tCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGcc                                                                                          >  1:854681/1‑77 (MQ=34)
                                                                                                                                    tCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                  <  1:949672/149‑1 (MQ=37)
                                                                                                                                         cAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACg             <  1:546722/149‑1 (MQ=37)
                                                                                                                                          aGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACgg            >  2:412676/1‑149 (MQ=37)
                                                                                                                                           gTCGGCACCATTGCCGGATGCGGCGTAACCGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa           <  1:517882/149‑1 (MQ=37)
                                                                                                                                           gTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAg                                                                                             >  2:542671/1‑67 (MQ=35)
                                                                                                                                              ggCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGa                                                                                       >  1:68947/1‑70 (MQ=35)
                                                                                                                                              ggCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGa                                                                                       <  2:68947/70‑1 (MQ=35)
                                                                                                                                               gCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGc                                                                                           >  1:521751/1‑65 (MQ=35)
                                                                                                                                               gCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGc                                                                                           <  2:521751/65‑1 (MQ=35)
                                                                                                                                               gCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGt       <  2:748735/149‑1 (MQ=37)
                                                                                                                                                  ccATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACTGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCg    <  2:653845/149‑1 (MQ=21)
                                                                                                                                                   cATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGc                                                                                           <  2:1095608/61‑1 (MQ=33)
                                                                                                                                                   cATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGc                                                                                           >  1:1095608/1‑61 (MQ=33)
                                                                                                                                                    aTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCggg  <  1:69333/149‑1 (MQ=32)
                                                                                                                                                    |                                                                                                                                                    
GAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGG  >  NC_000913/4325800‑4326096

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: