Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,905,785 T→G 12.5% intergenic (+63/+44) PP_5176 → / ← spuH hypothetical protein/spermidine/putrescine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,905,7850TG12.5% 450.8 / 2.2 177intergenic (+63/+44)PP_5176/spuHhypothetical protein/spermidine/putrescine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (77/71);  new base G (26/2);  total (104/73)
Fisher's exact test for biased strand distribution p-value = 3.56e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

GCCCGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCT  >  NC_002947/5905636‑5905934
                                                                                                                                                     |                                                                                                                                                     
gCCCGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCt                                                                                                                                                       <  2:3358910/150‑1 (MQ=255)
   cGGTACATTGGTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCtgtg                                                                                                                                                    <  2:2021554/150‑1 (MQ=255)
    ggTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGa                                                                                                                                                   <  2:1279778/150‑1 (MQ=255)
    ggTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGa                                                                                                                                                   <  1:4182499/150‑1 (MQ=255)
    ggTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCGCACCACAACCCTGAAGGCTGTGa                                                                                                                                                   >  2:3164794/1‑150 (MQ=255)
    cgattcagcaggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGa                                                                                                                                                   >  1:2027651‑M1/110‑150 (MQ=255)
     gTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAg                                                                                                                                                  <  1:156149/150‑1 (MQ=255)
     gTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCGCCCACAGGGCTCACCACAACCCAGAAGGCGGTGAg                                                                                                                                                  >  1:1442175/1‑150 (MQ=255)
      tACATTGTTGCCAATCCCTTGAGGGCTGGTTTGGTCAGATGTATCAGGGACTTTCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAgg                                                                                                                                                 <  1:1295410/150‑1 (MQ=255)
       aCATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGGGGGCATGCCCGCTCCCACAGTGCTCACCACAACCCTGAAGGCTGTGAGGt                                                                                                                                                <  1:2850553/150‑1 (MQ=255)
       aCATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGt                                                                                                                                                <  1:954341/150‑1 (MQ=255)
       aCATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGt                                                                                                                                                <  2:2482118/150‑1 (MQ=255)
         aTTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTcc                                                                                                                                              <  2:1520666/150‑1 (MQ=255)
         aTTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCAGAAGGCGGTGAGGTcc                                                                                                                                              >  2:629517/1‑150 (MQ=255)
           gcaggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCt                                                                                                                                            >  2:4305086‑M1/103‑150 (MQ=255)
              tgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                         <  2:389217/150‑1 (MQ=255)
              ggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                         >  1:3526459‑M1/100‑150 (MQ=255)
              ggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGACCTCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                        >  1:2924921‑M1/100‑150 (MQ=255)
               gcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                        >  2:1147240‑M1/99‑150 (MQ=255)
               gcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccccttcgcgggcatGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                        >  2:3926219‑M1/99‑150 (MQ=255)
               gCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATTCCCGCACCCACACGGCGCACCACAAACCAGAAGGCGGGGGGGTGCCAGTgg                                                                                                                                        >  1:4224270/1‑150 (MQ=255)
               gCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCGGTGAGGTCCCGGTgg                                                                                                                                        >  2:1018958/1‑150 (MQ=255)
                 cAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATTCCCGCGCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCGGTGGGa                                                                                                                                      >  2:770834/1‑150 (MQ=255)
                  aaTCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAg                                                                                                                                     <  2:1384606/150‑1 (MQ=255)
                       cTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCCCAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGc                                                                                                                                <  2:4515577/150‑1 (MQ=255)
                       cTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGc                                                                                                                                <  2:1140024/150‑1 (MQ=255)
                       agagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGc                                                                                                                                >  2:2347289‑M1/91‑150 (MQ=255)
                             ggCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCAAGCCCGCGCCCACAGGGCTCAACACAACCCTGAAGGCGTTTTGGTCCCGGTTGGAGAGGGCGTGccc                                                                                                                          >  2:413748/1‑150 (MQ=255)
                                   tttGTTCATAAGTATCAGGGACTATCCACATTGTGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGaaa                                                                                                                    <  1:3576252/149‑1 (MQ=255)
                                       gTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAgg                                                                                                                                                 <  1:2099217/117‑1 (MQ=255)
                                        tCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCCGGCATGCCCGCTCCCACCGGGCGCACCACAACCCAGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAATGGcc                                                                                                               >  2:1961844/1‑150 (MQ=255)
                                            aaGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCCTGCCCCCTCCCCCAGGGCGCACCACAACCCCGTAGGCGGTGAGGCCCCCGTGGGAGCGGGCGTGCCCCCGCAAGGGCCCGGt                                                                                                           >  1:2671720/1‑150 (MQ=255)
                                              gTATCAGGGACTAGCCACATTGGGTTGCGGTATGGCTTTGAATCGGTGCTGGCGCCTTCGCGGGCAGGCCCGCTCCCACAGGGCTCACCACAACCCTGACGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCa                                                                                                         <  1:1451348/150‑1 (MQ=38)
                                                      gACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCCCGCCCACAGGGCTCACCACAACCCAGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACg                                                                                                 >  2:1648546/1‑150 (MQ=255)
                                                      gACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCCGGCATGCCCGCGCCCACAGGGCTCAACACAACCCAGAAGGCGGGGAGGTCCCGGTGGGAGCGGGCGTGCCCCCGAAAGGGCCAGGTCCGGCCGACg                                                                                                 >  2:1648544/1‑150 (MQ=255)
                                                       aCTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGAATTCCCGCAACCAAAGGGCTAACAACAACCAGGAAGGCGGTGAGGTCCCTGGGGGAGCGGGCGTGCCCGCGCAAGGGCCAGGTCAGGCCGACGc                                                                                                >  2:1594312/1‑150 (MQ=255)
                                                         tATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAg                                                                                              >  2:2866128/1‑150 (MQ=255)
                                                              acaTTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCCCGCCCACACCGCTCACCACAAACCAGAAAGCTGTTATGTCCCGGTGGGCGCGGGCGGGGCCCCGAAACGGCCCCGGCCGGCCCACCCAGCcccc                                                                                         >  1:805456/1‑150 (MQ=255)
                                                                  tGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCAATCCCGCTCCCACACGGCTAAACACAAACCAGAAGGCTGTGAGGTCCCGGTGGGAGCGGGCGTGCCCCCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCAgcg                                                                                     >  1:1851508/1‑150 (MQ=255)
                                                                  ccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGa                                                                                                                                      >  1:332951‑M1/48‑101 (MQ=255)
                                                                  ccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGa                                                                                                                                      <  2:332951‑M1/54‑1 (MQ=255)
                                                                    ggATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGACGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGcc                                                                                   <  1:4083073/150‑1 (MQ=255)
                                                                       tGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCGGAAGGCGGGGAGGTCCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGGCAGGCCGACGCAGCCACCGGCGCCggg                                                                                >  2:3091310/1‑150 (MQ=39)
                                                                            tATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCACTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCtt                                                                           <  1:2550458/150‑1 (MQ=255)
                                                                                gCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCCCAGGGCTCACCCCAACCCTGAACGCGGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGAcc                                                                       >  2:3496750/1‑150 (MQ=255)
                                                                                 cTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACCGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGGGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCCCCGGCGCCGGGCGCTTGACCt                                                                      >  2:3441925/1‑150 (MQ=255)
                                                                                 cTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACACGGCTCACCACAAACCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCACGTCAGGCCGACGCAGCCACCAGCGCCGGGCGCTTGTcca                                                                      >  2:936684/1‑149 (MQ=38)
                                                                                  tttGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCCCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTg                                                                     >  2:2850762/1‑150 (MQ=255)
                                                                                  tttGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCAGAAGGCTGTGAGGTCCCTGGGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTg                                                                     >  2:4350971/1‑150 (MQ=255)
                                                                                     gAAGCGGTGCTGGAGCCGTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCgg                                                                  <  2:4390264/150‑1 (MQ=255)
                                                                                         cGGTGCTGGTGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGc                                                              <  2:4279845/150‑1 (MQ=255)
                                                                                              cgggcgCCTTCGCGGGCAGGCCCGCTCCCACCGGGCTCACCACACCCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGc                                                         <  1:192104/145‑1 (MQ=255)
                                                                                               tGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCg                                                        <  2:4427254/150‑1 (MQ=255)
                                                                                               tGGAGCCTTCGCGGGCATGCCCGCGCCCACAGGGCTCACCACAACCCAGAAGGCGGGGAGGTCCCTGGGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCCCCAGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCg                                                        >  1:4125625/1‑150 (MQ=39)
                                                                                                ggAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGt                                                       >  2:4020806/1‑150 (MQ=255)
                                                                                                ggAGCCTTCGCGGGCATGCCCCCTCCCACAGGGCGCACCACAACCCAGAACGCGGGGTGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCACGCCGACGCAGCCACCGGCGCCGGGCGCGTGACCTGCGGCTGCGGAGCgg                                                       >  2:532188/1‑149 (MQ=255)
                                                                                                 gcgCCTGCGCGGGCCTGCCCGCGCCCACCGGGCACACCACACCCCTGAAGGCGGTGCGGGCCCGGGGGGCGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGtt                                                      <  1:2408448/148‑1 (MQ=255)
                                                                                                  aGCCTCCGCGGGCAGGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTg                                                     <  1:2482047/150‑1 (MQ=255)
                                                                                                   gCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACCACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTgg                                                    >  2:1428154/1‑150 (MQ=255)
                                                                                                    ccTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGc                                                   <  1:4124098/150‑1 (MQ=255)
                                                                                                     cTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCt                                                  >  2:1059998/1‑150 (MQ=255)
                                                                                                        cgcgGGCATGCCCGCTCCCCCAGGGCTCACCACAACCCTGAAGGCGGTGAGGTCCCTGAGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgca                                               <  2:4228683/150‑1 (MQ=255)
                                                                                                        cgcgGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgca                                               <  1:3164794/150‑1 (MQ=255)
                                                                                                        cgcgGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgca                                               <  2:2168204/150‑1 (MQ=255)
                                                                                                           gggCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgcagca                                            >  1:3192625/1‑150 (MQ=255)
                                                                                                             gCATGCCCGCTCCCCCCGGGCACACCACAACCCTGACGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCAcc                                          <  1:2055785/150‑1 (MQ=255)
                                                                                                                tGCCCGCTCCCACAGGGCTCACCACAACCCTGAACGCGGTGTGGTCCCTGTGGGAGCGGGCGTGCCCCCGAAAGGGCCAGGTCAGCCCGACGCAGCCACCAGCGCCGGGCGCGAGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCctc                                       >  2:2219116/1‑150 (MQ=255)
                                                                                                                 gctcGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCctct                                      <  1:2990376/147‑1 (MQ=255)
                                                                                                                 GCCCGCTCCCACAGGGCTCACCACAACCCTGACGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCctct                                      >  1:2184044‑M1/1‑150 (MQ=255)
                                                                                                                  CCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCtt                                     <  2:2911835‑M1/150‑1 (MQ=255)
                                                                                                                  CCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCtt                                     <  2:3328286‑M1/150‑1 (MQ=255)
                                                                                                                   CCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTc                                    <  1:3219604‑M1/150‑1 (MQ=255)
                                                                                                                   CCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTc                                    <  1:2289265‑M1/150‑1 (MQ=255)
                                                                                                                         CCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGTCCCCTGCGCCGGGCGCTTGACCTGCGGCTGCGACGCGTTGGCTGCCGCTCCCTCTTCGATCGc                              <  1:4328495‑M1/150‑1 (MQ=255)
                                                                                                                         CCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgcag                                              >  2:4254325‑M1/1‑134 (MQ=255)
                                                                                                                         CCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGc                              <  2:63966‑M1/150‑1 (MQ=255)
                                                                                                                         CCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGc                              <  2:2896201‑M1/150‑1 (MQ=255)
                                                                                                                          CACAGGGCTCACCACAACCCTGACGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGcc                             >  1:247592‑M1/1‑150 (MQ=255)
                                                                                                                          CACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGcc                             <  2:1248266‑M1/150‑1 (MQ=255)
                                                                                                                           cccGGGCCCACCAAACCCCTGCCGGCTGTTTGGTCCCGGGGGGAGGGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCct                            >  1:1501913‑M1/4‑150 (MQ=255)
                                                                                                                               GGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCGGCACCCTCTTCGATCGCctgctg                        <  1:3420845‑M1/150‑1 (MQ=255)
                                                                                                                               GGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCctgctg                        <  1:3420838‑M1/150‑1 (MQ=255)
                                                                                                                                GCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGGGGGTGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGa                       >  1:524144‑M1/1‑150 (MQ=255)
                                                                                                                                 cTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGaa                      >  2:2738629/1‑150 (MQ=255)
                                                                                                                                     ccacAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGc                  >  1:2978767/1‑150 (MQ=255)
                                                                                                                                      cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc                 >  2:2336416/1‑150 (MQ=255)
                                                                                                                                      cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc                 >  1:2370838/1‑150 (MQ=255)
                                                                                                                                      cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc                 >  1:1785716/1‑150 (MQ=255)
                                                                                                                                      cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc                 >  1:2754036/1‑150 (MQ=255)
                                                                                                                                        caACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGTCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCtt               <  2:3308221/150‑1 (MQ=255)
                                                                                                                                        caACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCtt               >  1:186812/1‑150 (MQ=255)
                                                                                                                                        caACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCtt               >  2:3143333/1‑150 (MQ=255)
                                                                                                                                        caACCCTGAAGGCTGTGAGGTCCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCtt               <  1:3731543/150‑1 (MQ=255)
                                                                                                                                         aaCCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTg              >  2:3896595/1‑150 (MQ=255)
                                                                                                                                            ccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGTACCCTCTTCGATCGCCTGCTGAATCGCCTTGcga           >  1:1772500/1‑150 (MQ=255)
                                                                                                                                            ccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcga           >  2:612741/1‑150 (MQ=255)
                                                                                                                                             cTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgac          >  2:1121914/1‑150 (MQ=255)
                                                                                                                                             cTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgac          >  1:1006860/1‑150 (MQ=255)
                                                                                                                                              tGAAGGCTGTGAGGTCCCTGTGGTAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgacg         >  2:337424/1‑150 (MQ=255)
                                                                                                                                                     tgtgAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGctct  >  2:3845513/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCCGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCT  >  NC_002947/5905636‑5905934

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