Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,552,962 T→A 5.2% intergenic (+577/+84) 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,9620TA 5.2% 224.9 / 2.5 77intergenic (+577/+84)PP_4882/PP_5717iron ABC transporter permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/39);  new base A (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 = 9.33e-01

CGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTAC  >  NC_002947/5552814‑5553108
                                                                                                                                                    |                                                                                                                                                  
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                                                                                                                    <  2:2036546/127‑1 (MQ=255)
                                                      tgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtggt                                                                                                                    >  1:2036546/1‑127 (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                                                                                                     <  2:902169/126‑1 (MQ=255)
                                                                      cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGtt                                                                                                     >  1:902169/1‑126 (MQ=255)
                                                                       aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                >  1:1121544/1‑98 (MQ=255)
                                                                       aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                <  2:1121544/98‑1 (MQ=255)
                                                                         gaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTc                                                                               >  1:1251480/1‑145 (MQ=255)
                                                                         gaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTc                                                                               <  2:1251480/145‑1 (MQ=255)
                                                                         gaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGttt                                                                          <  2:876801/150‑1 (MQ=255)
                                                                             cAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGAcc                                                                      >  2:896488/1‑150 (MQ=255)
                                                                             cAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGAcc                                                                      >  2:2204683/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                                                                              >  2:1996127/1‑104 (MQ=255)
                                                                                                                   cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCt                                                                              <  1:1996127/104‑1 (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                                                                                                     <  2:1667762/78‑1 (MQ=255)
                                                                                                                     gAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTAGTGtt                                                                                                     >  1:1667762/1‑78 (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)
                                                                                                                                                    |                                                                                                                                                  
CGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTAC  >  NC_002947/5552814‑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: