Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,228,623 T→A 6.6% intergenic (+15/+15) PP_4605 → / ← PP_4606 AraC family transcriptional regulator/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,228,6230TA 6.6% 635.8 / 22.7 194intergenic (+15/+15)PP_4605/PP_4606AraC family transcriptional regulator/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (83/98);  new base A (5/8);  total (88/106)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 199 total aligned reads displayed.

CGGCTGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTTGTGTTC  >  NC_002947/5228477‑5228772
                                                                                                                                                  |                                                                                                                                                     
    tGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGccc                                                                                                                                                >  2:3981810/1‑150 (MQ=255)
    tGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTGTCATCTGAAACCGGGGCCGCAGAGCGGccc                                                                                                                                                <  2:3218995/150‑1 (MQ=255)
     gggAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGcccc                                                                                                                                               >  1:787713/1‑150 (MQ=255)
     gggAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCAGAGCGGcccc                                                                                                                                               <  2:4027735/150‑1 (MQ=255)
      ggAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCa                                                                                                                                              >  1:2667150/1‑150 (MQ=255)
            aTGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGcc                                                                                                                                        >  2:4735222/1‑150 (MQ=255)
            aTGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCAGTTGcc                                                                                                                                        >  1:3247047/1‑150 (MQ=255)
            aTGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTAGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCAGAGCGGccccggtttca                                                                                                                                        <  2:4190516/150‑8 (MQ=255)
              gCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTa                                                                                                                                      >  1:222639/1‑150 (MQ=255)
               ccGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTaa                                                                                                                                     <  1:3493105/150‑1 (MQ=255)
                       agagaTTGCCTTTTTGTGCGGGTTTGACATTGCCAGCCATTTCATCAACCTGTTTATCCAAGCCCTGGGGGCCTCACCTTGCGAATACCGCAAGGCCTTTTTCATCTTAATGCGGGTCCTCTCTGCGGCCCCTTTTGCCTAAATCTCgta                                                                                                                             <  2:1110718/150‑3 (MQ=255)
                        gagaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAg                                                                                                                            <  1:361728/150‑1 (MQ=255)
                          gaTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATGTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGt                                                                                                                          <  1:2931635/150‑1 (MQ=255)
                           aTTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTGTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTg                                                                                                                         <  2:3128694/150‑1 (MQ=255)
                             tGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCt                                                                                                                       >  2:1447293/1‑150 (MQ=255)
                             tGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCt                                                                                                                       <  1:2980794/150‑1 (MQ=255)
                              gCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTc                                                                                                                      <  1:1969684/150‑1 (MQ=255)
                              gCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTc                                                                                                                      <  1:2416207/150‑1 (MQ=255)
                              gCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTc                                                                                                                      <  2:4961861/150‑1 (MQ=255)
                              gCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTc                                                                                                                      <  1:1413956/150‑1 (MQ=255)
                                cTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGa                                                                                                                    >  1:2187046/1‑150 (MQ=255)
                                  tATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACa                                                                                                                  >  2:1424522/1‑150 (MQ=255)
                                  tATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACa                                                                                                                  >  2:5450510/1‑150 (MQ=255)
                                     tgtgCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAAc                                                                                                               <  1:5410810/150‑1 (MQ=255)
                                     tgtgCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTGCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTGTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAAc                                                                                                               <  1:3018357/150‑1 (MQ=255)
                                               tgccagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCcg                                                                                                     <  2:2351288/150‑1 (MQ=255)
                                                gccagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCcgc                                                                                                    <  1:1239549/150‑1 (MQ=255)
                                                 ccagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGTGGCCCCATTTGCCTAAAACTCGAGGTGCTCGACAAACTGACCTGc                                                                                                       >  2:3232578/1‑145 (MQ=255)
                                                  cagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCAttt                                                                                                                                           <  2:3711653/109‑1 (MQ=255)
                                                  cagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCAttt                                                                                                                                           >  1:3711653/1‑109 (MQ=255)
                                                  cagtgccagCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCgcgc                                                                                                  <  2:746150/150‑1 (MQ=255)
                                                     tgccagCCATTTCAGCAACCGGTGTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatc                                                                                               <  2:1192211/150‑1 (MQ=255)
                                                      gccagccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcg                                                                                              <  1:1530107/150‑1 (MQ=255)
                                                       ccagccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgc                                                                                             >  1:4477956/1‑150 (MQ=255)
                                                        cagccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTc                                                                                                                      <  1:5392532/124‑1 (MQ=255)
                                                        cagccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTc                                                                                                                      >  2:5392532/1‑124 (MQ=255)
                                                         agccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgcca                                                                                           >  1:1003138/1‑150 (MQ=255)
                                                          gccaTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaa                                                                                          >  2:1716224/1‑150 (MQ=255)
                                                            caTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatc                                                                                        <  2:2736519/150‑1 (MQ=255)
                                                             aTTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCctcgccaatcg                                                                                       >  1:5083189/1‑150 (MQ=255)
                                                                     aaCCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAg                                                                               >  1:3786284/1‑150 (MQ=255)
                                                                       ccGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTa                                                                             >  2:4681510/1‑150 (MQ=255)
                                                                          gTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTg                                                                          >  1:1599342/1‑150 (MQ=255)
                                                                            ttCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGt                                                                        <  2:1844496/150‑1 (MQ=255)
                                                                               gCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCa                                                                     >  1:4276432/1‑150 (MQ=255)
                                                                                ccAAGCCCTGGGGGCCACACCTGGCGACTACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCAc                                                                    >  2:3355537/1‑150 (MQ=255)
                                                                                ccAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCAc                                                                    <  1:1799947/150‑1 (MQ=255)
                                                                                ccAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCAc                                                                    <  2:2660402/150‑1 (MQ=255)
                                                                                 cAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCAcc                                                                   >  1:5232748/1‑150 (MQ=255)
                                                                                      ccTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCACTCGCCACAAAGTACTGGTTCACCGCCGa                                                              >  1:1016519/1‑150 (MQ=255)
                                                                                      ccTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGa                                                              <  1:2779529/150‑1 (MQ=255)
                                                                                      ccTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGa                                                              <  1:4487095/150‑1 (MQ=255)
                                                                                      ccTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCAACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGa                                                              <  1:399339/150‑1 (MQ=255)
                                                                                        tGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGCCAAACTGACCTGCCGCGCATCGCCACTCGCCACACAGTACGGGTGCACCGCCGAGc                                                            >  2:4072029/1‑150 (MQ=255)
                                                                                        tGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGc                                                            >  1:3326614/1‑150 (MQ=255)
                                                                                        tGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGc                                                            >  2:1974957/1‑150 (MQ=255)
                                                                                          ggggCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTg                                                          <  2:3856790/150‑1 (MQ=255)
                                                                                             gCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTaa                                                       <  2:446897/150‑1 (MQ=255)
                                                                                             gCCACACCTGGCGAATACCGCAACGCCTATTTCAACTGAAACCGGGGCCGCTCTGCGGCCCCAATTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTCCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGACCTGTaa                                                       <  2:872786/150‑1 (MQ=255)
                                                                                                  acCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcgc                                                                                      >  1:3224691/1‑114 (MQ=255)
                                                                                                   ccTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcgc                                                                                      <  1:4607550/113‑1 (MQ=255)
                                                                                                   ccTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcgc                                                                                      >  2:4607550/1‑113 (MQ=255)
                                                                                                     tGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCtt                                               >  1:4398230/1‑150 (MQ=255)
                                                                                                       gCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCatcgccaatcgc                                                                                      >  2:2957835/1‑109 (MQ=255)
                                                                                                         naaTACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCg                                           >  1:2913061/2‑150 (MQ=255)
                                                                                                               cGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAgg                                     >  1:4046716/1‑150 (MQ=255)
                                                                                                               cGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAgg                                     >  1:4102959/1‑150 (MQ=255)
                                                                                                               cGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAgg                                     <  1:2674383/150‑1 (MQ=255)
                                                                                                               cGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAgg                                     >  1:2109970/1‑150 (MQ=255)
                                                                                                                  aaGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGttc                                  <  1:4947777/150‑1 (MQ=255)
                                                                                                                  aaGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGttc                                  <  1:4107417/150‑1 (MQ=255)
                                                                                                                    ggCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGttctt                                >  2:120591/1‑150 (MQ=255)
                                                                                                                     gacTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCcac                                                                                   >  1:3881932/3‑98 (MQ=255)
                                                                                                                     gacTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCcac                                                                                   <  2:3881932/96‑1 (MQ=255)
                                                                                                                     gCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGttcttc                               <  2:2679288/150‑1 (MQ=255)
                                                                                                                     gCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGttcttc                               <  2:1401869/150‑1 (MQ=255)
                                                                                                                     gCCTTGTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGttcttc                               <  2:3450769/150‑1 (MQ=255)
                                                                                                                       cTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCAc                             >  1:664809/1‑150 (MQ=255)
                                                                                                                       cTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCAc                             >  2:1317804/1‑150 (MQ=255)
                                                                                                                        tttttCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACg                            >  2:1271486/1‑150 (MQ=255)
                                                                                                                            tCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCa                        <  1:707440/150‑1 (MQ=255)
                                                                                                                            tCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCa                        <  2:5316431/150‑1 (MQ=255)
                                                                                                                              aTCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGt                      <  1:120591/150‑1 (MQ=255)
                                                                                                                               tCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAAACGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGtt                     >  1:25513/1‑150 (MQ=255)
                                                                                                                                 tGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACa                                       <  2:2926737/130‑1 (MQ=255)
                                                                                                                                 tGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGc                   <  2:1352957/150‑1 (MQ=255)
                                                                                                                                 tGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAAATGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGc                   >  1:3937933/1‑150 (MQ=255)
                                                                                                                                  gcaaatggggCCGCAGAGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCa                  >  2:404268/7‑150 (MQ=255)
                                                                                                                                  gAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCa                  <  2:5142969/150‑1 (MQ=255)
                                                                                                                                  gAAACCGGGGCCGCTCAGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCa                  <  1:3932992/150‑1 (MQ=255)
                                                                                                                                    aaatggggCCGCAGAGCGGCCCCGTTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTg                                                                          >  2:2021914/5‑92 (MQ=38)
                                                                                                                                    aaatggggCCGCAGAGCGGCCCCGTTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTg                                                                          <  1:2021914/88‑1 (MQ=38)
                                                                                                                                    aaatggggCCGCAGAGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAgg                >  1:4965952/5‑150 (MQ=255)
                                                                                                                                      atggggCCGCAGAGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGgcg              >  1:5135589/3‑150 (MQ=255)
                                                                                                                                         gggCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCAc           <  2:5059153/150‑1 (MQ=255)
                                                                                                                                              gCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACaaa                                                                                <  2:4263802/76‑1 (MQ=255)
                                                                                                                                              gCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTtgt      <  1:1447293/150‑1 (MQ=255)
                                                                                                                                                 ctGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTTGTGtt   >  2:3330210/1‑150 (MQ=255)
                                                                                                                                                  tGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTTGTGTtc  <  2:2187046/150‑1 (MQ=255)
                                                                                                                                                  tGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTTGTGTtc  <  2:2118045/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGGCTGGGAGTAATGCCGCTGGGAGAGATTGCCTTATTGTGCGGGTTTGCCAGTGCCAGCCATTTCAGCAACCGGTTTCGCCAAGCCCTGGGGGCCACACCTGGCGAATACCGCAAGGCCTTTTTCATCTGAAACCGGGGCCGCTCTGCGGCCCCATTTGCCTAAAACTCGAAGGTGCTCGACAAACTGACCTGCCGCGCATCGCCAATCGCCACAAAGTACTGGTTCACCGCCGAGCTGTAATACACCTTGTCGAACAGGTTCTTCACGTTCAGTTGCAGGCGCACGTTGTGTTC  >  NC_002947/5228477‑5228772

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