Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,325,948 A→G 6.3% 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.3% 804.8 / 31.4 277intergenic (‑207/+451)yjdN/yjdMmetalloprotein superfamily protein/zinc‑ribbon family protein
Reads supporting (aligned to +/- strand):  ref base A (94/156);  new base G (3/14);  total (97/180)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-01

Only 100 of 279 total aligned reads displayed.

AGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCG  >  NC_000913/4325799‑4326094
                                                                                                                                                     |                                                                                                                                                  
  agGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCAGGATGCGGCTTTAACGCCTTATCCGGCCTACGGAGGGTTCGGTAATTTTTAGTCCTGATAAGACGCGCAATCTTCGCATCAGTCATTCGGCACCAtt                                                                                                                                                  <  1:2647899/150‑1 (MQ=14)
    gcggcCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGc                                                                                                                                                <  1:3633647/150‑1 (MQ=18)
    gcggcCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGc                                                                                                                                                <  2:2840838/150‑1 (MQ=38)
     cggcCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGcc                                                                                                                                               <  1:1197886/150‑1 (MQ=38)
     cggcCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGcc                                                                                                                                               <  1:919298/150‑1 (MQ=38)
     cggcCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGcc                                                                                                                                               <  2:1632322/150‑1 (MQ=38)
     cggcCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGcc                                                                                                                                               <  2:2261217/150‑1 (MQ=38)
          tGATTAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATg                                                                                                                                          <  1:4111767/150‑1 (MQ=37)
           gATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgc                                                                                                                                         <  2:448550/150‑1 (MQ=17)
            aTAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcg                                                                                                                                        >  2:2898383/1‑150 (MQ=17)
            aTAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATgcg                                                                                                                                        <  1:351787/150‑1 (MQ=37)
                gACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcggc                                                                                                                                      <  1:2615607/148‑1 (MQ=17)
                 aCGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCgg                                                                                                                                             >  1:3987433/1‑140 (MQ=17)
                 aCGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCgg                                                                                                                                             <  2:3987433/140‑1 (MQ=17)
                   gcgcgAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTaaa                                                                                                                                 <  1:3081007/150‑1 (MQ=37)
                      cgAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGc                                                                                                                              <  2:1794860/150‑1 (MQ=35)
                       gAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGcc                                                                                                                             <  2:491205/150‑1 (MQ=35)
                             cGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                       <  1:2130728/150‑1 (MQ=34)
                               cATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                     <  2:4070103/150‑1 (MQ=34)
                                  cAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                                                                  <  1:714452/150‑1 (MQ=34)
                                     gCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCATAACCTGCCAACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGAGGCGCATCTGGCAGTCGGCACCAATGCCGGATGCGGCGTACAGGCCTTATCCggggtggg                                                                                                               >  2:1430848/1‑144 (MQ=11)
                                       aGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGa                                                                                                             <  2:670658/150‑1 (MQ=34)
                                        ctCGGAACTGTTGCCGGATTCGTCGTAAACGACATATCAGGCCTAACGATGGTGCTGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAg                                                                                                            <  2:3524046/149‑1 (MQ=11)
                                         tCGGCACTGTTGCCGGATGCGGCGTAAACGGGTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAgg                                                                                                           <  1:3733321/150‑1 (MQ=18)
                                         tCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAgg                                                                                                           >  2:860276/1‑150 (MQ=33)
                                         tCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAgg                                                                                                           <  2:2101375/150‑1 (MQ=33)
                                         tCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAgg                                                                                                           <  1:1741647/150‑1 (MQ=33)
                                              aCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                      <  1:112789/150‑1 (MQ=32)
                                              aCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                      <  1:3772134/150‑1 (MQ=32)
                                              aCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                      <  1:2937872/150‑1 (MQ=32)
                                              aCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                      <  2:2018905/150‑1 (MQ=32)
                                              aCTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                      <  2:540431/150‑1 (MQ=32)
                                                   tGCCGGATGCGGCGTAAACGCCTTATCCGGCCTCCGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAt                                                                                                 <  2:2232345/150‑1 (MQ=17)
                                                     ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAttt                                                                                               <  1:1242354/150‑1 (MQ=32)
                                                        gATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTa                                                                                            <  1:2495607/150‑1 (MQ=32)
                                                        gATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTa                                                                                            <  1:1571911/150‑1 (MQ=32)
                                                        gATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTa                                                                                            <  2:295303/150‑1 (MQ=32)
                                                            cggcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAttt                                                                                               <  2:2193594/143‑1 (MQ=32)
                                                            cggcgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAttt                                                                                               >  1:2193594/1‑143 (MQ=32)
                                                               cgTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGTCGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGa                                                                                     >  1:3603440/1‑150 (MQ=17)
                                                                gTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGAt                                                                                    <  1:1898210/150‑1 (MQ=32)
                                                                     cGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGGGGGAATTTGTAGGCCTGATAAGAc                                                                               <  2:3544819/150‑1 (MQ=17)
                                                                       ccTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                             <  1:1327934/150‑1 (MQ=32)
                                                                            tCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTg                                                                                                                                                 <  2:3371614/77‑1 (MQ=34)
                                                                            tCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTg                                                                                                                                                 >  1:3371614/1‑77 (MQ=34)
                                                                              cGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg                                                                      >  2:1128102/1‑150 (MQ=32)
                                                                              cGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCg                                                                      >  2:2752376/1‑150 (MQ=32)
                                                                                   tACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCa                                                                 >  2:469639/1‑150 (MQ=32)
                                                                                     cGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                      <  2:188732/111‑1 (MQ=32)
                                                                                     cGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                                               <  2:3284385/150‑1 (MQ=32)
                                                                                      ggAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa                                                              <  2:1080431/150‑1 (MQ=32)
                                                                                       gAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTaaa                                                                                                                                 >  2:4087552/1‑82 (MQ=33)
                                                                                       gAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                             <  2:660587/150‑1 (MQ=32)
                                                                                           gttcGGGATTTTGTCGGCCTGATAAGACGCGCTAGCGTCGCATCAGGCAGTCGGCACCTTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAg                                                         <  2:1549060/147‑1 (MQ=11)
                                                                                                ggAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAg                                                                                           <  2:3442255/111‑1 (MQ=32)
                                                                                                 gAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATc                                                               >  2:2384141/1‑138 (MQ=32)
                                                                                                 gAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGCAAGAGACGTATCCTGCCTCCGGAGGGTGCGGGCAGTAGAAGGCCTGGTAAGACGCGCACGCGGCGCTTCAGGCAgaaggcc                                                   >  2:515955/1‑144 (MQ=11)
                                                                                                  aaTTTGTCGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCTGTCGGCAc                                                  <  2:3113462/150‑1 (MQ=255)
                                                                                                   aTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGACTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCAcc                                                 >  1:3049104/1‑150 (MQ=37)
                                                                                                        tAGGCCTGATAAGACGCTCAAGCGTCGCATCAGGCAGTCGGCTCCTTTGCCGGAGTCGTCGTTAACGCCTTATCCGGCCTCCGGAGGGTGCGGGATTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGTCACCGTTGc                                            <  2:3964531/150‑1 (MQ=11)
                                                                                                           gCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGAcgc                                                                             >  1:2545174/1‑114 (MQ=32)
                                                                                                             cTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGAt                                       <  1:3020657/150‑1 (MQ=37)
                                                                                                              tGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCg                                                                                                      >  2:1248205/1‑86 (MQ=33)
                                                                                                              tGATAAGACGCGCAAGCGTCGCATCAGGCAGGCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGGGCGGGACGGGGAAGGCCAGATAACAAGCCCAAGCGGAGCATCTGGCAGGCGGCACCGTTGCCGGagg                                      >  1:459390/1‑148 (MQ=11)
                                                                                                                 tAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcgg                                   <  2:3667512/150‑1 (MQ=37)
                                                                                                                  aaGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcggc                                  <  2:1692432/150‑1 (MQ=37)
                                                                                                                  aaGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcggc                                  >  2:3897962/1‑150 (MQ=37)
                                                                                                                   aGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCTTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATgcggcg                                 <  1:406680/150‑1 (MQ=21)
                                                                                                                     aCGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGc                                                                                         >  1:423095/1‑92 (MQ=33)
                                                                                                                     aCGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCg                                                                   >  1:3156595/1‑114 (MQ=32)
                                                                                                                     aCGCGCAAGCGTCGCATCAGGCAGACGGCACCATTGCCGGATGCGGCGGAAACGCCTTATCCGGCCTACGGAGGGAGCGGGAATTTGTAGGCCGGATAAGACGCGCAAGCGTCGCAACAGGCAGTGGGCACCGTTGCCGGATgcggcggg                               >  1:203232/1‑148 (MQ=11)
                                                                                                                      cgcgCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTaa                              <  1:3797278/150‑1 (MQ=37)
                                                                                                                       gcgcAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAtt                                                                                                >  1:1768076/1‑83 (MQ=33)
                                                                                                                       gcgcAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAtt                                                                                                <  2:934347/83‑1 (MQ=33)
                                                                                                                       gcgcAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAtt                                                                                                >  1:934347/1‑83 (MQ=33)
                                                                                                                       gcgcAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAtt                                                                                                <  2:1767906/83‑1 (MQ=33)
                                                                                                                         gcAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa                                                              >  1:3925895/1‑115 (MQ=32)
                                                                                                                          ccaGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGc                          <  2:2092936/148‑1 (MQ=37)
                                                                                                                          cAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGc                          <  1:2243166/150‑1 (MQ=37)
                                                                                                                           aaGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGc                                                                       <  1:818040/104‑1 (MQ=32)
                                                                                                                           aaGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGc                                                                       >  2:818025/1‑104 (MQ=32)
                                                                                                                           aaGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACCGAGGGGGCGGGAATATGTAGGACTGATAATACTCGCACGCGTCGCATCAGGCAGTCGGCACCGAAGACGGATGCGGCGTAATCGcc                         >  1:2715667/1‑150 (MQ=11)
                                                                                                                                tCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGcc                                                                                                                             >  1:1969432/1‑45 (MQ=38)
                                                                                                                                tCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGcc                                                                                                                             <  2:1969432/45‑1 (MQ=38)
                                                                                                                                tCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAggg                                                                                                          <  2:3793116/64‑1 (MQ=35)
                                                                                                                                 cGCATCAGGCAGTCGGCACCCTTGCCGGATGCGGCGTACACGCCTTATCCGGCCTACGGAGGGTGCGTGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCTCATCAGGCAGTCGGCACCGTTGCCGGATGATGCTTAAACGCCTTATcc                   <  2:2353630/150‑1 (MQ=11)
                                                                                                                                 cGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATcc                   <  1:2898610/150‑1 (MQ=37)
                                                                                                                                    aTCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                <  1:2416509/150‑1 (MQ=37)
                                                                                                                                     tCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGACGTTGTAGGCCTGATAAGACGCGCACGCGTCGCATCAGGCCGGCGGCACCGGTGCCGGATGCGGCGGGAGCGCCGGATCCGGcc               >  1:2355395/1‑150 (MQ=11)
                                                                                                                                         gCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCat                                                         <  2:263080/104‑2 (MQ=14)
                                                                                                                                            gTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAt        <  1:344823/150‑1 (MQ=38)
                                                                                                                                               ggCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGACTTTGTAGGCCTGATAAGa                                                                                >  1:2292333/1‑75 (MQ=21)
                                                                                                                                                gCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAgg                                                                                                           >  2:3296661/1‑47 (MQ=38)
                                                                                                                                                gCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAgg                                                                                                           <  1:3296661/47‑1 (MQ=38)
                                                                                                                                                gCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTa                                                                                            <  1:3405091/62‑1 (MQ=35)
                                                                                                                                                gCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCTGTCGGCACCGGTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTg    >  2:3405059/1‑150 (MQ=255)
                                                                                                                                                gCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTg    >  2:70294/1‑150 (MQ=17)
                                                                                                                                                  aCCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGATTTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCg  <  2:2587468/150‑1 (MQ=17)
                                                                                                                                                  aCCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCg  <  1:2312960/150‑1 (MQ=35)
                                                                                                                                                  aCCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCg  <  1:65395/150‑1 (MQ=35)
                                                                                                                                                     |                                                                                                                                                  
AGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCATTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAGTCGGCACCGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCG  >  NC_000913/4325799‑4326094

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: