Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,552,961 T→G 5.2% intergenic (+576/+85) PP_4882 → / ← PP_5717 iron ABC transporter permease/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,552,9610TG 5.2% 212.8 / 3.1 77intergenic (+576/+85)PP_4882/PP_5717iron ABC transporter permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/39);  new base G (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

TTCGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTAC  >  NC_002947/5552812‑5553108
                                                                                                                                                     |                                                                                                                                                   
ttCGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCt                                                                                                                                                     <  2:1704129/150‑1 (MQ=255)
  cGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACACGGTCTTGAGGAATTCCGACCCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTc                                                                                                                                                   <  1:1370531/150‑1 (MQ=255)
  cGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTc                                                                                                                                                   <  1:2225608/150‑1 (MQ=255)
  cGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTc                                                                                                                                                   <  1:1410579/150‑1 (MQ=255)
       cTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCaaa                                                                                                                                              <  2:1573713/150‑1 (MQ=255)
           ttttCGCGCTCTCATAAGTCACTGCGCATCTGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGCCCGGCAGGTCCGGGCTTTGCTTCGCGAAGccc                                                                                                                                          <  2:412540/150‑1 (MQ=255)
            tttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCg                                                                                                                                         <  2:612791/150‑1 (MQ=255)
            tttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCg                                                                                                                                         <  2:2255415/150‑1 (MQ=255)
             ttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCgg                                                                                                                                        >  2:1787309/1‑150 (MQ=255)
                gcgcTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGaaa                                                                                                                                     >  2:292543/1‑150 (MQ=255)
                      tcATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACANCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgggg                                                                                                                               >  1:1954921/1‑150 (MQ=255)
                        aTAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                <  2:1555714/147‑1 (MQ=255)
                        aTAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                >  1:1555714/1‑147 (MQ=255)
                             tCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGCAACTGGGGTTTCAGc                                                                                                                        >  2:795703/1‑150 (MQ=255)
                               aCTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCtt                                                                                                                                                    <  1:100271/120‑1 (MQ=255)
                               aCTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCtt                                                                                                                                                    >  2:100271/1‑120 (MQ=255)
                                cTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtgg                                                                                                                     <  2:1478580/150‑1 (MQ=255)
                                     cATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgct                                                                                                                >  1:2160925/1‑150 (MQ=255)
                                      aTCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGccc                                                                                                                                          >  2:2230318/1‑123 (MQ=255)
                                      aTCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGccc                                                                                                                                          <  1:2230318/123‑1 (MQ=255)
                                          gTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGCCCg                                                                                                                                         <  2:434979/120‑1 (MQ=18)
                                          gTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGCCCg                                                                                                                                         >  1:434979/1‑120 (MQ=18)
                                             ggTACGCGTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGTTCTTGAGGAATCCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAgt                                                                                                        <  1:951601/150‑1 (MQ=255)
                                               tACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAgtgt                                                                                                      <  2:2160925/150‑1 (MQ=255)
                                               tACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAgtgt                                                                                                      <  1:292543/150‑1 (MQ=255)
                                                  gTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTa                                                                                                   <  2:451037/150‑1 (MQ=255)
                                                       gtgtCCAGCCGGGTCACCTAGAGGAAGGCTTTGGCAGAACAAGGTCTTACGCAATTCCGAACCAGATGCCAAGACCGGCGGGTCCGCGCTTTTCTTCGCAAAGCCCGGTAACTGGGGTTTCACCTAGTGGCTGCTTGAGTGTTTAGaaaa                                                                                              <  2:1122340/150‑1 (MQ=255)
                                                       gtgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGaaaa                                                                                              <  1:1569109/150‑1 (MQ=255)
                                                       gtgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTGTAGaaaa                                                                                              >  1:1122340/1‑150 (MQ=255)
                                                        tgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtggt                                                                                                                    >  1:2036546/1‑127 (MQ=255)
                                                        tgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtggt                                                                                                                    <  2:2036546/127‑1 (MQ=255)
                                                        tgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAAc                                                                                             <  1:2285501/150‑1 (MQ=255)
                                                         gtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                          <  1:2344497/136‑1 (MQ=255)
                                                         gtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                          >  2:2344497/1‑136 (MQ=255)
                                                             aGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAAATCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATgg                                                                                        >  2:1142556/1‑150 (MQ=255)
                                                                      aCCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTc                                                                               >  1:1424805/1‑150 (MQ=255)
                                                                      aCCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGATTGGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTc                                                                               >  1:1222411/1‑150 (MQ=255)
                                                                        cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGtt                                                                                                     >  1:902169/1‑126 (MQ=255)
                                                                        cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGtt                                                                                                     <  2:902169/126‑1 (MQ=255)
                                                                         aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                >  1:1121544/1‑98 (MQ=255)
                                                                         aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                <  2:1121544/98‑1 (MQ=255)
                                                                           gaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTc                                                                               <  2:1251480/145‑1 (MQ=255)
                                                                           gaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTc                                                                               >  1:1251480/1‑145 (MQ=255)
                                                                           gaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGttt                                                                          <  2:876801/150‑1 (MQ=255)
                                                                               cAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGAcc                                                                      >  2:2204683/1‑150 (MQ=255)
                                                                               cAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGAcc                                                                      >  2:896488/1‑150 (MQ=255)
                                                                                 ggCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAg                                                                    <  1:1098911/150‑1 (MQ=255)
                                                                                   cTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCATTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGCGGGTTTAGAAAACAATGGGTTGGATTCTGTTTGACCAggg                                                                  >  1:256882/1‑150 (MQ=255)
                                                                                         cTGAACGAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTGAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTTCAATCTGATTGACCAGGGCTTAcc                                                            >  2:1215083/1‑150 (MQ=18)
                                                                                          tGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCt                                                           >  2:2046229/1‑150 (MQ=255)
                                                                                          tGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGAGCTTGTCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCt                                                           >  1:299621/1‑150 (MQ=255)
                                                                                              cAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGCGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAg                                                       >  1:438661/1‑150 (MQ=255)
                                                                                              cAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAg                                                       <  1:2093683/150‑1 (MQ=255)
                                                                                                 ggTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTAc                                                    <  2:1424805/150‑1 (MQ=255)
                                                                                                      tGAGGAACTCCGAACCTGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTCACCTCAAGTACTCTAt                                               <  2:1222411/150‑1 (MQ=255)
                                                                                                       gcggAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATgg                                                                                        >  1:1046841/3‑108 (MQ=255)
                                                                                                       gAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCtt                                                                                                                                                    <  1:425572/48‑1 (MQ=255)
                                                                                                       gAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCtt                                                                                                                                                    >  2:425572/1‑48 (MQ=255)
                                                                                                       gAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATgg                                                                                        <  2:1046841/108‑1 (MQ=255)
                                                                                                          gAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCaa                                           <  1:896488/150‑1 (MQ=255)
                                                                                                           aaTTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCaaa                                          <  1:2191862/150‑1 (MQ=255)
                                                                                                                  aaCCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAgtgt                                                                                                      >  2:848620/1‑83 (MQ=255)
                                                                                                                  aaCCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAgtgt                                                                                                      <  1:848620/83‑1 (MQ=255)
                                                                                                                   aCCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaa                                  <  1:1215083/150‑1 (MQ=255)
                                                                                                                    ccAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaaaaaa                                 <  1:1787309/150‑1 (MQ=255)
                                                                                                                     cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCt                                                                              <  1:1996127/104‑1 (MQ=255)
                                                                                                                     cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCt                                                                              >  2:1996127/1‑104 (MQ=255)
                                                                                                                     cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGt                                                                            >  2:1547741/1‑106 (MQ=255)
                                                                                                                     cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGt                                                                            <  1:1547741/106‑1 (MQ=255)
                                                                                                                     cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGa                                                                        <  2:1723985/110‑1 (MQ=255)
                                                                                                                     cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGa                                                                        >  1:1723985/1‑110 (MQ=255)
                                                                                                                       gAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGa                                                                        >  2:1439168/1‑108 (MQ=255)
                                                                                                                       gAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGa                                                                        <  1:1439168/108‑1 (MQ=255)
                                                                                                                       gAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTAGTGtt                                                                                                     >  1:1667762/1‑78 (MQ=255)
                                                                                                                       gAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTAGTGtt                                                                                                     <  2:1667762/78‑1 (MQ=255)
                                                                                                                             aaGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGGTGCATTCTGTTTGACCAGGGCTTACCTCACGTACTCTATCCAACGAAAAAAAAGAACAATCa                        >  1:1213467/1‑150 (MQ=255)
                                                                                                                                            gggCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCaaaa         >  1:2050363/1‑150 (MQ=255)
                                                                                                                                              gCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGt                      <  2:994042/135‑1 (MQ=255)
                                                                                                                                              gCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGt                      >  1:994042/1‑135 (MQ=255)
                                                                                                                                               cTTTGCTTCGCAGAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATggg                                                                                       <  2:1256498/69‑1 (MQ=255)
                                                                                                                                               cTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATggg                                                                                       >  1:1256498/1‑69 (MQ=255)
                                                                                                                                                   gCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTAc  >  1:1032971/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTCGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTAC  >  NC_002947/5552812‑5553108

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: