Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,909,967 C→G 7.7% intergenic (+315/+164) PP_4318 → / ← PP_4319 transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,909,9670CG 7.7% 231.2 / 10.0 78intergenic (+315/+164)PP_4318/PP_4319transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/39);  new base G (2/4);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

CCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAG  >  NC_002947/4909818‑4910116
                                                                                                                                                       |                                                                                                                                                      
ccTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTACGGGGTTCCCTCGATACGGAGTTCGGCATGGCGTGACGCGCGGCAGGCCccatcc                                                                                                                                                        >  2:155135/1‑150 (MQ=11)
 cTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCccatcca                                                                                                                                                       >  2:118607/1‑150 (MQ=37)
 cTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCccatcca                                                                                                                                                       >  1:866588/1‑150 (MQ=38)
        cccTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTTCGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCAAGGAGca                                                                                                                                                <  2:634020/150‑2 (MQ=12)
         ccTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCt                                                                                                                                               >  1:691779/1‑150 (MQ=32)
           ttACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGAGGCCCCATGGATGGGGCCtgca                                                                                                                                            <  2:980001/150‑4 (MQ=11)
           cgttctgcttctgcgtctgcgtctgcttctgcttttgctgttgcttctaagcgcgcgatagttcaggcgacgcagattgcgacttcaggaggccgagcgtagggcttgcTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCa                                                                                                                                            <  2:1408857‑M2/41‑1 (MQ=255)
               ggcggGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGAGGCCCCATGGATGGGGCCTGCAgcgc                                                                                                                                        <  2:1079063/150‑1 (MQ=11)
                gcggGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATGGATGGGGCCTGCAgcgc                                                                                                                                        <  1:155560/150‑1 (MQ=255)
                 cggGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGTGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCtt                                                                                                                                      <  2:691779/150‑1 (MQ=17)
                      cACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACgg                                                                                                                                 <  2:1289336/150‑1 (MQ=32)
                         tttttGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCa                                                                                                                              >  2:768166/1‑150 (MQ=32)
                         tttttGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCa                                                                                                                              <  1:488650/150‑1 (MQ=32)
                               tCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATg                                                                                                                        >  2:1261902/1‑150 (MQ=32)
                                 aaaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGa                                                                                                                                    <  1:597672/136‑1 (MQ=33)
                                 aaaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGa                                                                                                                                    >  2:597672/1‑136 (MQ=33)
                                  aaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGccc                                                                                                                     >  2:465798/1‑150 (MQ=32)
                                  aaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGccc                                                                                                                     <  2:39070/150‑1 (MQ=32)
                                  aaCGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGccc                                                                                                                     >  1:504451/1‑150 (MQ=32)
                                    cgcgACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGt                                                                                                                   >  1:208368/1‑150 (MQ=32)
                                     ccaccAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCTCTCCTTTACGCTCCGGGGTTCCCTCGCTCAACTTTTGCTCCGGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTc                                                                                                                  <  1:319111/146‑1 (MQ=11)
                                       gACAAAAAGTAACCAAAACACGCTAGATGTTATCACCTTGCCCTGGCTCGACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCaa                                                                                                                <  2:866588/150‑2 (MQ=11)
                                         cAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATGGATGCCCGTCaagc                                                                                                              <  1:322299/150‑4 (MQ=255)
                                            aaaGTAACCAACAACCGCTGGATCCTCTCGGCCGGCCGTTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTccc                                                                                                           <  2:1224628/150‑1 (MQ=11)
                                              tgcttttgcttgtaagcgcgcggtagttcaggcgacgcagattgcgacttcagaaggccgagcgcagggcttgcGAAGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTc                                                                                                        <  2:242732‑M2/76‑1 (MQ=255)
                                               gTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTc                                                                                                        <  2:195798/150‑1 (MQ=34)
                                                 aaCCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTc                                                                                                        <  2:1038498/148‑1 (MQ=33)
                                                 aaCCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTc                                                                                                        >  1:1038498/1‑148 (MQ=34)
                                                     aaaaaaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATggg                                                                                                                                                   >  1:1130204/1‑102 (MQ=33)
                                                     aaaaaaCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATggg                                                                                                                                                   <  2:1130204/102‑1 (MQ=33)
                                                         ctaagtgcgcggtagttcaggcgacgcagattgcgacttcaggaggccgagcgcagggcttgcTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAAGCCCTgc                                                                                             <  1:399384‑M2/87‑1 (MQ=255)
                                                         ctaagcgcgcgatagcccaggcgccgcggattgcgacttcaggaggccgagcgcagggcttgcTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTgc                                                                                             <  2:98245‑M2/87‑1 (MQ=255)
                                                          aCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAgc                                                                                                                                          >  2:1146106/1‑105 (MQ=33)
                                                          aCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAgc                                                                                                                                          <  1:1146106/105‑1 (MQ=33)
                                                          aCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTgc                                                                                             <  1:155135/150‑1 (MQ=38)
                                                          aCGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTgc                                                                                             >  1:109604/1‑150 (MQ=38)
                                                                  tCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGcc                                                                                     >  2:716317/1‑150 (MQ=38)
                                                                  tCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCACAACCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGcc                                                                                     >  1:956539/1‑150 (MQ=255)
                                                                    cTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                  <  2:204891/119‑1 (MQ=31)
                                                                    cTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTc                                                                                                                  >  1:204891/1‑119 (MQ=31)
                                                                      atcatcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCt                                                                                 <  2:113932/150‑1 (MQ=37)
                                                                       tcatcCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTg                                                                                <  1:584024/150‑1 (MQ=38)
                                                                        ttcaggcgacgcagattgcgacttcaggaggccgagcgtagggcttgcTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGccc                                                                                                                     <  1:946522‑M2/64‑1 (MQ=255)
                                                                        ttcaggcgacgcagattgcgacttcaggaggccgagcgtagggcttgcTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGccc                                                                                                                     >  2:946522‑M2/49‑112 (MQ=255)
                                                                                cccTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCaa                                                                       >  1:1081860/1‑150 (MQ=38)
                                                                                  cTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATc                                                                     <  2:956278/150‑1 (MQ=39)
                                                                                     cgcTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGc                                                                  >  2:878866/1‑150 (MQ=39)
                                                                                          cGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGc                                                             <  1:878866/150‑1 (MQ=39)
                                                                                                       ccTCGCTACGGTCTTGCTCCCGGGAGTACTCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACcgcg                                                >  1:1138935/1‑150 (MQ=255)
                                                                                                       ccTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACcgcg                                                >  2:1426356/1‑150 (MQ=255)
                                                                                                        cTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACcgcgc                                               >  2:712108/1‑150 (MQ=255)
                                                                                                         tCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCACCATCCATAGGGCCTCAGCGCTTTAGCGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAAATACcgcgcg                                              >  2:712315/1‑150 (MQ=14)
                                                                                                             tACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAGCTTGCAATCTGCGTCGCCTCAACTACCGCGCGCTTa                                          >  2:1142891/1‑150 (MQ=255)
                                                                                                                ggTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                       <  1:465798/150‑1 (MQ=255)
                                                                                                                ggTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                       <  2:1341071/150‑1 (MQ=255)
                                                                                                                 gTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCtt                                                                                                         <  1:836625/83‑1 (MQ=31)
                                                                                                                 gTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCtt                                                                                                         >  2:836625/1‑83 (MQ=31)
                                                                                                                        tCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGAACGGGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTacaagcacaagc                               >  1:331549/1‑150 (MQ=14)
                                                                                                                          ccGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                       <  1:1185433/140‑1 (MQ=255)
                                                                                                                          ccGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                       >  2:1185433/1‑140 (MQ=255)
                                                                                                                            gggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaa                           <  1:1267838/150‑1 (MQ=255)
                                                                                                                             ggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGc                                                                                      <  2:589348/90‑1 (MQ=35)
                                                                                                                             ggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGc                                                                                      >  1:589348/1‑90 (MQ=35)
                                                                                                                             ggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGc                                                                  >  2:1158628/1‑110 (MQ=38)
                                                                                                                             ggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGc                                                                  <  1:1158628/110‑1 (MQ=38)
                                                                                                                             ggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                       >  1:720847/1‑137 (MQ=255)
                                                                                                                             ggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                       <  2:720847/137‑1 (MQ=255)
                                                                                                                             ggAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaag                          <  2:504451/150‑1 (MQ=255)
                                                                                                                                  tcaagcgcTGAGGCCCCATGGATGGGGCCTGCAGCGCGTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTc                                                                                         >  1:1417160/5‑84 (MQ=11)
                                                                                                                                  tcaagcgcTGAGGCCCCATGGATGGGGCCTGCAGCGCGTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTc                                                                                         <  2:1417160/80‑1 (MQ=11)
                                                                                                                                  tcaagcgcTGAGGCCCCATGGATGGGGCCT‑‑GCGCTTGACGGGCATCCATGCCCGTCAACTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaag                    >  2:1406651/5‑150 (MQ=17)
                                                                                                                                     cgcgCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaag                    >  2:721534/1‑150 (MQ=255)
                                                                                                                                        gcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                       >  2:746623/1‑128 (MQ=255)
                                                                                                                                        gcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACaa                                       <  1:746623/128‑1 (MQ=255)
                                                                                                                                        gcTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcag                 <  2:1406906/150‑1 (MQ=255)
                                                                                                                                             aGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaagca            <  1:721534/150‑1 (MQ=255)
                                                                                                                                              ggCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAaagcagaagcagaagcagaagcag           <  1:1426356/150‑1 (MQ=255)
                                                                                                                                                ccccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagca         >  2:830274/1‑150 (MQ=255)
                                                                                                                                                ccccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagca         >  2:612357/1‑150 (MQ=255)
                                                                                                                                                 cccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagcag        <  1:716317/150‑1 (MQ=255)
                                                                                                                                                 cccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagcag        >  2:1110745/1‑150 (MQ=255)
                                                                                                                                                  ccatccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagcaga       >  1:638916/1‑150 (MQ=255)
                                                                                                                                                    atccatGGGGCCTCAGCGCTTGACGGGCATCCATGTCCTGCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAagcagaagcagcagaag     <  1:1142891/150‑1 (MQ=255)
                                                                                                                                                       catGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGcagaagcagcagaagcag  >  1:1283448/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
CCTTTCCGCCCTTACGGCGGGTCACTTTTTGTCAAACGCGACAAAAAGTAACCAAAAAACGCTGGCTCCTATCATCCGGCCCCTGCGCTACGCTCCGGGGTTCCCTCGCTACGGTCTTGCTCCCGGGAGTACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTTCGCAAGCCCTGCGCTCGGCCTTCTGAAGTCGCAATCTGCGTCGCCTGAACTACCGCGCGCTTACAAGCAAAAGCAGAAGCAGAAGCAGAAGCAGCAGAAGCAG  >  NC_002947/4909818‑4910116

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