Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,552,963 C→A 12.0% intergenic (+578/+83) 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,9630CA12.0% 261.4 / 18.0 92intergenic (+578/+83)PP_4882/PP_5717iron ABC transporter permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (37/44);  new base A (5/6);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

GTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAGTC  >  NC_002947/5552816‑5553112
                                                                                                                                                   |                                                                                                                                                     
gTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGc                                                                                                                                                     <  1:2904268/150‑1 (MQ=255)
    tCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAg                                                                                                                                                 >  1:2352147/1‑150 (MQ=255)
    tCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAg                                                                                                                                                 >  1:513422/1‑150 (MQ=255)
     cATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGc                                                                                                                                                >  1:651752/1‑150 (MQ=255)
      aTTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGcc                                                                                                                                               <  1:855671/150‑1 (MQ=255)
      aTTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGcc                                                                                                                                               >  1:2653543/1‑150 (MQ=255)
       ttttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGccc                                                                                                                                              <  2:1237361/150‑1 (MQ=255)
       ttttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGccc                                                                                                                                              <  2:442207/150‑1 (MQ=255)
       ttttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGccc                                                                                                                                              <  2:2653543/150‑1 (MQ=255)
        tttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCg                                                                                                                                             <  1:1341032/150‑1 (MQ=255)
            gcgcTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGGTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCGGAACCAGAAGCCGACACCAGCAGGTCCGGGCTTTGCTTCGCAGAGCCCGGaga                                                                                                                                         >  2:243811/1‑148 (MQ=17)
            gcgcTCTCATAAATCACCGCGCCTTAGCTGGAAAGACTGGAGAATAGCCTGCGCACTAAGAAGCAGGCTTAGGCTGAACAAGGTCTTGAGGAATTCCGAACAAGAAGCCAAGACCGGCAGGTCCGGGCTTTTCTTCGCAAAACCCGGaaa                                                                                                                                         >  2:1571907/1‑150 (MQ=17)
               ctctCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTg                                                                                                                                      <  2:2352147/150‑1 (MQ=255)
                 ccgatctgTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                    <  1:2193410/143‑1 (MQ=255)
                      aaGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                    <  2:275183/145‑1 (MQ=255)
                      aaGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                    >  1:275183/1‑145 (MQ=255)
                         tCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCGTTTATTCTTAAAGCCCGGAAACTGGGGtgaacgg                                                                                                                            >  2:2193410/1‑144 (MQ=17)
                             tGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCATGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAg                                                                                                                                                 <  1:2289307/125‑1 (MQ=255)
                             tGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCATGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAg                                                                                                                                                 >  2:2289307/1‑125 (MQ=255)
                                  aTCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGccc                                                                                                                                              >  1:681556/1‑123 (MQ=255)
                                  aTCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGccc                                                                                                                                              <  2:681556/123‑1 (MQ=255)
                                   tCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgctgc                                                                                                                  <  1:1080150/150‑1 (MQ=255)
                                     aGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                                <  2:2104912/150‑1 (MQ=255)
                                       ttGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                    >  2:1715181/1‑128 (MQ=255)
                                       ttGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                    <  1:1715181/128‑1 (MQ=255)
                                       ttGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGATGGTGGCTGCTTGa                                                                                                              <  1:2818914/150‑1 (MQ=255)
                                       ttGGTACGTCTGGTGTACAGCCGGTTCAACAAGAAGCATGCGTATGCTGAACAACGGCCTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                              <  1:1571907/150‑1 (MQ=18)
                                         ggTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAgt                                                                                                            <  2:651752/150‑1 (MQ=255)
                                          gTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTc                                                                                                                                                       >  1:806545/1‑106 (MQ=255)
                                          gTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTc                                                                                                                                                       <  2:806545/106‑1 (MQ=255)
                                           tACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtggtgg                                                                                                                      >  2:376763/1‑138 (MQ=255)
                                           tACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtggtgg                                                                                                                      <  1:376763/138‑1 (MQ=255)
                                            aCGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGGTACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGtt                                                                                                         <  1:1780533/150‑1 (MQ=255)
                                                 tGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGaa                                                                                                    >  1:1554125/1‑150 (MQ=255)
                                                   gtgtACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGaaaa                                                                                                  <  1:553892/150‑1 (MQ=255)
                                                       aCAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAAt                                                                                              <  2:1409277/150‑1 (MQ=255)
                                                           ccGGTTCACAAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTGTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGt                                                                                          >  1:27875/1‑150 (MQ=255)
                                                              gTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                                <  2:2768547/125‑1 (MQ=255)
                                                              gTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                                >  1:2768547/1‑125 (MQ=255)
                                                                  aCCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGa                                                                                                                                           >  2:2611917/1‑94 (MQ=255)
                                                                  aCCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGa                                                                                                                                           <  1:2611917/94‑1 (MQ=255)
                                                                   ccAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCt                                                                                  >  1:2141037/1‑150 (MQ=255)
                                                                   ccAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGCCCGGacc                                                                                                                                         <  2:2724724/95‑3 (MQ=17)
                                                                   ccAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGCCCGGacc                                                                                                                                         >  1:2724724/1‑93 (MQ=17)
                                                                    cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCg                                                                                                                                                      >  2:1047644/1‑81 (MQ=255)
                                                                    cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCg                                                                                                                                                      <  1:1047644/81‑1 (MQ=255)
                                                                    cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCgg                                                                                                                                            <  2:2093504/91‑1 (MQ=255)
                                                                    cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCgg                                                                                                                                            >  1:2093504/1‑91 (MQ=255)
                                                                    cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAgt                                                                                                            >  2:1591043/1‑123 (MQ=255)
                                                                    cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAgt                                                                                                            <  1:1591043/123‑1 (MQ=255)
                                                                    cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATgg                                                                                            <  1:2699122/139‑1 (MQ=255)
                                                                    cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATgg                                                                                            >  2:2699122/1‑139 (MQ=255)
                                                                    cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTg                                                                                 >  1:295551/1‑150 (MQ=255)
                                                                    cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGAACCGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTg                                                                                 >  2:418211/1‑150 (MQ=255)
                                                                     aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGTTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGt                                                                                >  1:1411825/1‑150 (MQ=255)
                                                                     aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGt                                                                                <  2:1554125/150‑1 (MQ=255)
                                                                     aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGCCCGGa                                                                                                                                           >  1:2309687/1‑91 (MQ=17)
                                                                     aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGCCCGGa                                                                                                                                           <  2:2309687/91‑1 (MQ=17)
                                                                         agCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGa                                                                            <  2:513422/150‑1 (MQ=255)
                                                                                ttATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGtt                                                                                                         >  1:6526/1‑114 (MQ=255)
                                                                                ttATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGtt                                                                                                         <  2:6526/114‑1 (MQ=255)
                                                                                    gCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTAc                                                                 <  2:1681869/150‑1 (MQ=255)
                                                                                     cTGAACAAGGTCTTGAGGCCTTCAGAACCCGCTGCCAAGCCCGGCAGGCCCGGGCTTTGCTCCGCAAAGGCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTAcc                                                                <  2:1411825/150‑1 (MQ=17)
                                                                                     cTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCg                                                                                                                                             <  2:1189583/73‑1 (MQ=255)
                                                                                     cTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCg                                                                                                                                             >  1:1189583/1‑73 (MQ=255)
                                                                                     cTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGGAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGGTTAGAAAAAAATGGGTTGCATTCTGTTTGACCAGGGCTTTcc                                                                >  2:2992527/1‑150 (MQ=255)
                                                                                      tGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGttt                                                                              >  2:2124178/1‑135 (MQ=255)
                                                                                      tGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGtct                                                                              <  1:2124178/135‑3 (MQ=255)
                                                                                              gTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGccc                                                                                                                                              >  1:69497/1‑63 (MQ=255)
                                                                                              gTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGccc                                                                                                                                              <  2:69497/63‑1 (MQ=255)
                                                                                               tCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTActc                                                      >  1:1955400/1‑150 (MQ=255)
                                                                                                   gAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATc                                                  >  1:2058927/1‑150 (MQ=255)
                                                                                                   gAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGGGGCTGCTTGAGTGTTTAGAAAACAATGGGCTGAATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATc                                                  <  1:243811/150‑1 (MQ=255)
                                                                                                      gAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCaa                                               <  2:2058927/150‑1 (MQ=255)
                                                                                                         ttCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                              <  1:60260/84‑1 (MQ=255)
                                                                                                         ttCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGa                                                                                                              >  2:60260/1‑84 (MQ=255)
                                                                                                            cGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGt                                                                                >  1:2872025/1‑111 (MQ=255)
                                                                                                            cGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGt                                                                                <  2:2872025/111‑1 (MQ=255)
                                                                                                            cGAACCAGAAACCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAAt                                                                                              <  2:1474144/97‑1 (MQ=255)
                                                                                                            cGAACCAGAAACCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAAt                                                                                              >  1:1474144/1‑97 (MQ=255)
                                                                                                                   gAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATc                                                  <  2:748155/134‑1 (MQ=255)
                                                                                                                   gAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATc                                                  >  1:748155/1‑134 (MQ=255)
                                                                                                                          aGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTTAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCAAAAGAAAAAAAAGAACAATCAg                           <  2:2141037/150‑1 (MQ=255)
                                                                                                                                gCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATggg                                                                                           <  1:100490/80‑1 (MQ=255)
                                                                                                                                gCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATggg                                                                                           >  2:100490/1‑80 (MQ=255)
                                                                                                                                gCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAAt                     <  1:2992527/150‑1 (MQ=255)
                                                                                                                                       cGGGCTTTGTTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGa                                   >  2:898372/1‑129 (MQ=255)
                                                                                                                                       cGGGCTTTGTTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGa                                   <  1:898372/129‑1 (MQ=255)
                                                                                                                                       cGGGCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTc                                                              <  1:11514/102‑1 (MQ=21)
                                                                                                                                       cGGGCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTc                                                              >  2:11514/1‑102 (MQ=21)
                                                                                                                                            tttGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTg         >  2:2175739/1‑150 (MQ=255)
                                                                                                                                               gCTGCGCAAATCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATGCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTAc      <  1:684136/150‑1 (MQ=255)
                                                                                                                                                cTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCa                   >  1:1282239/1‑136 (MQ=255)
                                                                                                                                                cTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCa                   <  2:1282239/136‑1 (MQ=255)
                                                                                                                                                   cGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAGTc  <  2:1955400/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGGTGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAGTC  >  NC_002947/5552816‑5553112

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: