Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,905,785 T→G 11.9% 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,7850TG11.9% 363.9 / 2.2 145intergenic (+63/+44)PP_5176/spuHhypothetical protein/spermidine/putrescine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (61/58);  new base G (18/4);  total (83/62)
Fisher's exact test for biased strand distribution p-value = 9.52e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

CCCGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCT  >  NC_002947/5905637‑5905934
                                                                                                                                                    |                                                                                                                                                     
cccGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGACCCCGCCCCCACCGCTAACCACAAACCACGAGGCtg                                                                                                                                                      >  1:4164086/1‑150 (MQ=255)
 ccGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCtgt                                                                                                                                                     <  1:1028601/150‑1 (MQ=255)
  cGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCtgtg                                                                                                                                                    <  1:798730/150‑1 (MQ=255)
      ttcagcaggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGACCCCACAACCCTGAAGGCTGTGAGGTc                                                                                                                                               >  2:2054217‑M1/107‑150 (MQ=255)
        aTTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTcc                                                                                                                                              >  2:4628297/1‑150 (MQ=255)
          tgttgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGTACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCt                                                                                                                                            <  1:1456747/150‑1 (MQ=255)
          gcaggcgatcgaagagggtgctgcagccaacgattcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCt                                                                                                                                            >  2:5011858‑M1/103‑150 (MQ=255)
           gttgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtg                                                                                                                                           <  1:1311041/150‑1 (MQ=255)
            ttgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgt                                                                                                                                          <  1:4628297/150‑1 (MQ=255)
            ttgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCCCGCCCACAGGGCGCACCACAACCCCGAAGGCGGTGGGGTgccggg                                                                                                                                          >  1:236163/1‑147 (MQ=255)
            aggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgt                                                                                                                                          >  2:3714688‑M1/101‑150 (MQ=255)
            aggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgt                                                                                                                                          >  1:2171422‑M1/101‑150 (MQ=255)
            aggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgt                                                                                                                                          >  1:2809046‑M1/101‑150 (MQ=255)
             tgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                         <  2:4046308/150‑1 (MQ=255)
             tgCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                         <  2:4046085/150‑1 (MQ=255)
             ggcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                         >  1:228699‑M1/100‑150 (MQ=255)
              gcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                        >  2:2601196‑M1/99‑150 (MQ=255)
              gcgatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcgggcacGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                        >  1:4875090‑M1/99‑150 (MQ=255)
              gCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCGCAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                        <  2:1963717/150‑1 (MQ=255)
              gCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                        <  1:2669564/150‑1 (MQ=255)
              gCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTgg                                                                                                                                        <  2:1877873/150‑1 (MQ=255)
                gatcgaagagggtgctgcagccaacgcttcgcagccgcaggtcaagcgcccggcgcaggtggctgcgtcggcctgacctggccctttcgcggtcacGCCCGCTCCCACAGTGCTCACCACAACCCTGAAGGCTGTGAGTTCCCTGTGGGa                                                                                                                                      >  2:886484‑M1/97‑150 (MQ=255)
                cAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCCCGCCCCCACGGCTCCCCACCACCCAGAAGGCGGTGAGGGCCCCGGGGGa                                                                                                                                      >  1:541379/1‑150 (MQ=255)
                 aaTCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGGGGGAg                                                                                                                                     >  2:4358792/1‑150 (MQ=255)
                  aTCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCCTGCCCCCGCCCCCACGGCTCACCCCAACCCAGAAGGCGGTGAGGTGCCAGTGGGAGc                                                                                                                                    >  1:4639656/1‑150 (MQ=255)
                     ccTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCAAGCCCGCGCCCACACCGCTCAACACAACCCAGAAGGCTGTGAGGTCCCAGTGGGGGCggg                                                                                                                                 >  1:4729779/1‑150 (MQ=255)
                      cTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCCAGCCCCCTCCCCCAGGGCTCACCACAAACCAGAAGGCGGTGGGGTCCCGGTGGGAGCCGGc                                                                                                                                >  1:1900209/1‑150 (MQ=255)
                      cTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCtgtg                                                                                                                                                    <  1:4449709/130‑1 (MQ=255)
                       ttGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCGCCCACAGGGCACACCACAACCCAGAAGGCAGTGTGGTCCCTGTGGGAGCGGGCg                                                                                                                               >  1:5182281/1‑150 (MQ=255)
                            ggCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCGGTGAGGTCCCGGTGGGAGCGGGCGTGccc                                                                                                                          >  1:3598263/1‑150 (MQ=255)
                            ggCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCCCGCCCACACGGCTCACCACAACCCAGAAGGCTGTGAGGTACCAGGGGGAGCGGGCGTGccc                                                                                                                          >  1:528197/1‑150 (MQ=255)
                              cTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCCCCGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCcgc                                                                                                                        <  2:1292182/150‑1 (MQ=255)
                               tGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGATGTGAGGTCCCTGTGGGAGCGGGCGTGCCcgcg                                                                                                                       >  1:4044110/1‑150 (MQ=255)
                                      gTCAGATGTATCAGGGACTATCCACATTGGGATGCGGTAGGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACCGGGCTCACCACACCCCTGAAGGCTGTGAGGTCCCGGGGGGAGCGGGGGTGCCCGCGAAAGGGc                                                                                                                <  2:1880585/150‑1 (MQ=255)
                                       tCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGACAGGGcc                                                                                                               >  1:5218567/1‑150 (MQ=255)
                                       gcAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATTCCCCCGGCCACAAGGCGCACCACAACACACAAAGCAGTGAGGTGCCGGTGGGAGCGGGCGGGGCCGCCCAAAGGcc                                                                                                               >  1:5007741/2‑150 (MQ=255)
                                         aGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCGCCCACAGGGCTCACCACAACCCGGAAGGCGGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAg                                                                                                             >  2:4641739/1‑150 (MQ=255)
                                         aGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCCGGCAAGCCCGCTACCACCGGGCGCACCCCCCCCCAGACGGCGCTGTGGGCCCGGTGGGAGCGGGCGTGCCCGCGGAAGGGcccg                                                                                                             >  1:1423976/1‑148 (MQ=255)
                                                tCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCGCAACACAAACCGGAATGCTGTGAGGTCCCGGGGGGAGCGGGCGTGCCCCCGAAAGGGCCAGGTCAGGc                                                                                                      >  1:1790139/1‑150 (MQ=255)
                                                  aGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCg                                                                                                    >  2:828969/1‑150 (MQ=255)
                                                   gggACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCAACACAACCCTGAAGGCGGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGa                                                                                                   >  2:4264195/1‑150 (MQ=255)
                                                    ggACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGAc                                                                                                  >  2:4837502/1‑150 (MQ=255)
                                                        tATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCtgtg                                                                                                                                         >  1:803115/1‑107 (MQ=255)
                                                        tATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTCAAGGCTGTGAGGTCCCGGTGGGCGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAg                                                                                              >  1:3592197/1‑150 (MQ=255)
                                                          tCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCCCTCCCCCAGGGCGCACCCCAACCCAGAAGGCGGTGAGGTCCCGGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTAAGGCCGACGCAGcc                                                                                            >  1:4314307/1‑150 (MQ=255)
                                                                  gggATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTgcgc                                                                                    <  2:2066913/150‑1 (MQ=255)
                                                                    gATGCGTTATGGCTTTTATTCGGTGCGGGCGCCTTCGGGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGCGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCg                                                                                  <  1:3032162/150‑1 (MQ=38)
                                                                           tttGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACCGGGCTCACCACACCCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCtt                                                                           <  1:2092029/148‑1 (MQ=255)
                                                                           tATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCAGGCCCGCTCCCACAGGGCTCACCACCACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCtt                                                                           <  2:4092450/150‑1 (MQ=255)
                                                                                 tttGACTGGGTGGGGGAGCCTTCGCGGGCAGGCCCGCGCCCACCGGGCTCCCCACACCCCTGAAGGCGGGGCGGGCCCTGTGGGCGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTg                                                                     <  1:5098501/150‑1 (MQ=255)
                                                                                 tttGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCGGAAGGCGGTGAGGTCCCGGTGGGAGCGGGCGTGCCCGCCAAAGGGCCAGGTCAGGCCGACGCAGCCCCCAGCGCCGGGCGCGTGACCTg                                                                     >  2:2340797/1‑150 (MQ=39)
                                                                                 tttGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACACGGCTCACCACAACCCTGAAGGCGGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCGTGACCTg                                                                     >  2:4814918/1‑150 (MQ=255)
                                                                                     aaTCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGt                                                                                                                              <  2:1619852/89‑1 (MQ=255)
                                                                                     aaTCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGt                                                                                                                              >  1:1619852/1‑89 (MQ=255)
                                                                                           tGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACCGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGaa                                                           <  2:4727987/150‑1 (MQ=255)
                                                                                           tGCTGGAGCCTTCGCGGGCATGCCCCCGCCCACAGGGCTCACCACAACCCGGAACGCGGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGAAGCACCCCCCATCGCCGGGCGCTTGAACGGCGGCTGCGaa                                                           >  2:1011969/1‑150 (MQ=255)
                                                                                             cTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGc                                                         <  2:2506331/150‑1 (MQ=255)
                                                                                                gcgCCTGCGCGGGCCTGCCCGCTCCCACCGGGCACACCACAACCCTGATGGCTGGGAGGTCCCTGGGGGCGGGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCATCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGtt                                                      <  1:553349/148‑1 (MQ=255)
                                                                                                gcgCCTGCGCGGGCATGCCCGCTCCCACCGGGCACACCACACCCCTGCCGGCTGTGAGGTCCCTGGGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGtt                                                      <  1:1479717/148‑1 (MQ=255)
                                                                                                     ttCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTg                                                 <  1:4805099/150‑1 (MQ=255)
                                                                                                     ttCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCGTGACCTGCGGCTGCGAAGCGTTGGCTg                                                 >  1:1591973/1‑150 (MQ=255)
                                                                                                        gcgGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgcag                                              <  2:567783/150‑1 (MQ=255)
                                                                                                         cgGGCATGCCCGCTCCCACCGGGCTCACCACAACCCGGAAGGCTGTGAGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgcagc                                             <  2:541379/150‑1 (MQ=255)
                                                                                                              aTGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCct                                        >  1:2208065/1‑150 (MQ=255)
                                                                                                                  CCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTc                                    <  2:5058431‑M1/150‑1 (MQ=255)
                                                                                                                     CTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCATGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGAt                                 <  2:5358770‑M1/150‑1 (MQ=255)
                                                                                                                      TCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATc                                <  1:5444092‑M1/150‑1 (MQ=255)
                                                                                                                      TCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATc                                <  1:5444087‑M1/150‑1 (MQ=255)
                                                                                                                         CACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGcc                             >  2:3089987‑M1/1‑150 (MQ=255)
                                                                                                                          ACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCct                            <  1:4638506‑M1/150‑1 (MQ=255)
                                                                                                                          ACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCct                            <  1:4817626‑M1/150‑1 (MQ=255)
                                                                                                                                 tCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCcgc                                                                                                                        <  1:3849198/51‑1 (MQ=255)
                                                                                                                                 tCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCcgc                                                                                                                        >  2:3849198/1‑51 (MQ=255)
                                                                                                                                 tCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTg                                                               >  2:1514456/1‑108 (MQ=255)
                                                                                                                                 tCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAAt                     >  1:5066756/1‑150 (MQ=255)
                                                                                                                                 tCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAAt                     >  1:4002230/1‑150 (MQ=255)
                                                                                                                                 tCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAAt                     >  1:962405/1‑150 (MQ=255)
                                                                                                                                  caccacAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAACACCCTCTTCGATCGCCTGCTGAATc                    >  1:3161133/1‑150 (MQ=255)
                                                                                                                                   accacAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCg                   >  1:4577171/1‑150 (MQ=255)
                                                                                                                                   accacAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCg                   >  2:2214981/1‑150 (MQ=255)
                                                                                                                                     cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCctgc                          >  2:2116100/1‑141 (MQ=255)
                                                                                                                                     cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCctgc                          <  1:2116100/141‑1 (MQ=255)
                                                                                                                                     cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc                 >  2:1775945/1‑150 (MQ=255)
                                                                                                                                     cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc                 >  1:3884852/1‑150 (MQ=255)
                                                                                                                                     cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc                 >  1:2837408/1‑150 (MQ=255)
                                                                                                                                     cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc                 >  1:2853755/1‑150 (MQ=255)
                                                                                                                                       ccaCCCTGAAGGCGGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCtt               <  2:1962986/148‑1 (MQ=255)
                                                                                                                                       caACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCtt               >  1:68366/1‑150 (MQ=255)
                                                                                                                                       caACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCtg               >  2:1596386/1‑149 (MQ=255)
                                                                                                                                          cccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATc                                <  1:1106049/130‑1 (MQ=255)
                                                                                                                                           ccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGaa                                                           >  2:2963702/1‑102 (MQ=255)
                                                                                                                                           ccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGaa                                                           <  1:2963702/102‑1 (MQ=255)
                                                                                                                                           ccTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcga           >  1:4703096/1‑150 (MQ=255)
                                                                                                                                            cTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgac          >  1:5308353/1‑150 (MQ=255)
                                                                                                                                             tGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgacg         >  2:649467/1‑150 (MQ=255)
                                                                                                                                             tGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgacg         >  2:811424/1‑150 (MQ=255)
                                                                                                                                              gAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgacga        >  2:2321792/1‑150 (MQ=255)
                                                                                                                                              gAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgacga        <  1:4837406/150‑1 (MQ=255)
                                                                                                                                               aaGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCgacgac       <  1:4666834/150‑1 (MQ=255)
                                                                                                                                                    tgtgAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGctct  >  2:5187824/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCCGGTACATTGTTGCCAATCCCTTGAGGGCTGGATTGGTCAGAAGTATCAGGGACTATCCACATTGGGATGCGGTATGGCTTTGAATCGGTGCTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCT  >  NC_002947/5905637‑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: