Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 677,557 C→A 9.5% K431N (AAG→AAT hscC ← Hsp70 family chaperone Hsc62, binds to RpoD and inhibits transcription

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273677,5570CA 9.5% 882.5 / 45.0 242K431N (AAG→AAThscCHsp70 family chaperone Hsc62, binds to RpoD and inhibits transcription
Reads supporting (aligned to +/- strand):  ref base C (114/105);  new base A (10/13);  total (124/118)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 250 total aligned reads displayed.

CTTCAATCTGCTGCGCGCTCAATGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGCACAGTGGTGTTTCGTT  >  CP009273/677409‑677705
                                                                                                                                                    |                                                                                                                                                    
ctTCAATCTGCTGCGCGCTCAATGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGANTATCAATCGACTGATAAGCCCCGGTTttct                                                                                                                                                     <  1:4277326/150‑1 (MQ=255)
  tCAATCTGCTGCGCGCTCAATGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTttcttc                                                                                                                                                   >  2:3055460/1‑150 (MQ=255)
     aTCTGCTGCGCGCTCAATGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAAc                                                                                                                                                <  1:4208133/150‑1 (MQ=255)
     aTCTGCTGCGCGCTCAATGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAAc                                                                                                                                                <  1:316842/150‑1 (MQ=255)
           tgcGCGCTCAATGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGcaca                                                                                                                                          >  1:3474481/1‑150 (MQ=255)
            gcgcgcTCTATGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTCCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGTATATCACTCGACTGATAAGCCCCGGTTTTCTTAAACGGCACAt                                                                                                                                         <  2:4582814/150‑1 (MQ=255)
              gcgcTCAATGTTACCGGGCTGTGGTTAATCACTCTGGACATAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCg                                                                                                                                       >  2:509042/1‑150 (MQ=255)
                gcTCAATGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGaa                                                                                                                                     <  2:3079133/150‑1 (MQ=255)
                 cTCAATGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAg                                                                                                                                    >  2:3638490/1‑150 (MQ=255)
                  tCAATGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAgg                                                                                                                                   <  1:2919689/150‑1 (MQ=255)
                  tCAATGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAgg                                                                                                                                   >  1:170756/1‑150 (MQ=255)
                    aaTGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGAt                                                                                                                                 >  1:1289477/1‑150 (MQ=255)
                    aaTGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGAt                                                                                                                                 >  1:3207567/1‑150 (MQ=255)
                    aaTGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGAt                                                                                                                                 >  2:4042576/1‑150 (MQ=255)
                      tGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTc                                                                                                                               <  2:873624/150‑1 (MQ=255)
                                   tgGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGACGCACGTCAACTTCAAGCAACCCGTTGATATAATAACTAAAGCGAATTTCAATAGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAAtgttg                                                                                                                  <  2:2199563/150‑1 (MQ=255)
                                   tgGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAAtgttg                                                                                                                  <  1:2438058/150‑1 (MQ=255)
                                         aTCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTa                                                                                                            <  1:3715425/150‑1 (MQ=255)
                                            aCTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACt                                                                                                         <  2:3195468/150‑1 (MQ=255)
                                                 ggACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTtgtg                                                                                                    >  2:262805/1‑150 (MQ=255)
                                                  gACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTgg                                                                                                   >  2:2662258/1‑150 (MQ=255)
                                                    cTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGtt                                                                                                 >  1:3967257/1‑150 (MQ=255)
                                                     ttAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGttt                                                                                                >  2:1858860/1‑150 (MQ=255)
                                                     ttAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGCTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGttt                                                                                                >  1:1982029/1‑150 (MQ=255)
                                                             gCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTg                                                                                        <  1:434024/150‑1 (MQ=255)
                                                             gCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTg                                                                                        >  2:4455151/1‑150 (MQ=255)
                                                             gCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTg                                                                                        >  2:1078978/1‑150 (MQ=255)
                                                             gCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTg                                                                                        >  2:4288086/1‑150 (MQ=255)
                                                              ccGTCTTCCTGAAGCACGTCCTCTTCAAGCAACCCTTTGATATCATAACTAAAGCGAATATCACTAGACTGATAAGCCCCGGTTTTCTTCTACTGCAGATCGAAGGATTCGACCAGAATGTTGTTTTTATCTTTGTGTTTTTCTCCCTGa                                                                                       <  2:1262215/150‑1 (MQ=255)
                                                              ccGTCTTCCAGAAGCACGTCAACTTCATGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGACGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGa                                                                                       <  2:542869/150‑1 (MQ=255)
                                                              ccGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGa                                                                                       <  2:388456/150‑1 (MQ=255)
                                                                   ttCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGAc                                                                                  >  1:657935/1‑150 (MQ=255)
                                                                      cAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGtt                                                                               <  1:1647832/150‑1 (MQ=255)
                                                                       aGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTa                                                                              <  1:914240/150‑1 (MQ=255)
                                                                          aGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCGCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTACCTGATAGACGTTAAcc                                                                           >  1:970180/1‑150 (MQ=255)
                                                                          aGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAAcc                                                                           <  2:4163931/150‑1 (MQ=255)
                                                                           gCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCg                                                                          >  2:2343546/1‑150 (MQ=255)
                                                                            cACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGt                                                                         <  1:2271593/150‑1 (MQ=255)
                                                                               gTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAAt                                                                      >  1:751462/1‑150 (MQ=255)
                                                                               gTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAAt                                                                      >  2:923739/1‑150 (MQ=255)
                                                                               gTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAAt                                                                      >  2:337244/1‑150 (MQ=255)
                                                                                tCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAAtt                                                                     <  1:231010/150‑1 (MQ=255)
                                                                                 cAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTg                                                                    >  1:2371756/1‑150 (MQ=255)
                                                                                 cAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTATTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTg                                                                    <  2:3718774/150‑1 (MQ=255)
                                                                                    cTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAAt                                                                 >  1:1201294/1‑150 (MQ=255)
                                                                                     ttCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATc                                                                >  1:233440/1‑150 (MQ=255)
                                                                                       cAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCt                                                              >  1:3779802/1‑150 (MQ=255)
                                                                                         aGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGt                                                            >  1:2908011/1‑150 (MQ=255)
                                                                                          gCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGtt                                                           <  1:2272627/150‑1 (MQ=255)
                                                                                          gCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTGTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGtt                                                           >  2:673354/1‑150 (MQ=255)
                                                                                           cAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTc                                                          <  1:577731/150‑1 (MQ=255)
                                                                                           cAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTc                                                          >  2:2007990/1‑150 (MQ=255)
                                                                                            aaCCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTcc                                                         <  1:4603881/150‑1 (MQ=255)
                                                                                            aaCCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTATTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTcc                                                         <  1:400905/150‑1 (MQ=255)
                                                                                             aCCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCg                                                        <  2:3005730/150‑1 (MQ=255)
                                                                                             aCCCGTTGATATCAGAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCg                                                        <  2:1948261/150‑1 (MQ=255)
                                                                                              cccGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTACTTGAATCCTGTTCCgg                                                       >  2:1936538/1‑150 (MQ=255)
                                                                                               ccGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCggg                                                      >  2:386395/1‑150 (MQ=255)
                                                                                               ccGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCggg                                                      <  1:3331920/150‑1 (MQ=255)
                                                                                                 gTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTg                                                    >  1:1048398/1‑150 (MQ=255)
                                                                                                   tGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCa                                                  >  2:1583043/1‑150 (MQ=255)
                                                                                                    gATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCAt                                                 >  1:1513167/1‑150 (MQ=255)
                                                                                                     atatCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATg                                                >  2:692393/1‑150 (MQ=255)
                                                                                                       atCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGt                                              <  1:3559218/150‑1 (MQ=255)
                                                                                                       atCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTATTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGt                                              >  2:775207/1‑150 (MQ=255)
                                                                                                       atCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTATTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGt                                              >  2:2932349/1‑150 (MQ=255)
                                                                                                       atCATAACTAAAGCGAATATCAATCGACCGATAAGCCCCGGTTTTCTTCAACTTCACATCGAAGGTTTCCACGAGTATGTTGTTTTTTACTTTGTGGTTTTCTCCCTGATAGCCGTTATCCGTAATTGATTCCTGTTCCGGGTGCATGGt                                              <  2:4399778/150‑1 (MQ=255)
                                                                                                        tCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGtt                                             <  2:645921/150‑1 (MQ=255)
                                                                                                         cATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTg                                            <  1:718928/150‑1 (MQ=255)
                                                                                                          ctAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCATCGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTACCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGa                                           <  2:1159454/149‑1 (MQ=255)
                                                                                                          atcaCTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGa                                           <  1:3819515/147‑1 (MQ=255)
                                                                                                          aTAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGa                                           >  1:3404529/1‑150 (MQ=255)
                                                                                                              cTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATaa                                       >  1:2488412/1‑150 (MQ=255)
                                                                                                                 aaGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGtt                                    >  2:1374123/1‑150 (MQ=255)
                                                                                                                 aaGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGtt                                    >  2:260643/1‑150 (MQ=255)
                                                                                                                 aaGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGtt                                    >  2:2401977/1‑150 (MQ=255)
                                                                                                                  aGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGttt                                   <  2:1531464/150‑1 (MQ=255)
                                                                                                                    cGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCt                                 >  1:2768157/1‑150 (MQ=255)
                                                                                                                     gAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTa                                <  1:47580/150‑1 (MQ=255)
                                                                                                                     gAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTa                                >  2:1676684/1‑150 (MQ=255)
                                                                                                                      aaTCTCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTac                               <  1:849922/150‑1 (MQ=255)
                                                                                                                        tatCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTacac                             <  1:1955772/150‑1 (MQ=255)
                                                                                                                        tatCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTacac                             >  2:2081894/1‑150 (MQ=255)
                                                                                                                           cAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACAcgcg                          >  2:893687/1‑150 (MQ=255)
                                                                                                                            aaTCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGa                         >  1:1277615/1‑150 (MQ=255)
                                                                                                                               cGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGacac                      >  2:3236977/1‑150 (MQ=255)
                                                                                                                                gACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACg                     <  1:1557808/150‑1 (MQ=255)
                                                                                                                                gACTGATAAGCCCCGGTTTTATTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACg                     <  1:2204430/150‑1 (MQ=255)
                                                                                                                                   tGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGc                  >  2:1615317/1‑150 (MQ=255)
                                                                                                                                   tGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGc                  <  1:3909124/150‑1 (MQ=255)
                                                                                                                                   tGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGc                  >  1:4485478/1‑150 (MQ=255)
                                                                                                                                    gATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGca                 <  2:198789/150‑1 (MQ=255)
                                                                                                                                    gATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGca                 >  1:3307335/1‑150 (MQ=255)
                                                                                                                                    gATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGca                 <  1:1188727/150‑1 (MQ=255)
                                                                                                                                      tAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGcaca               <  2:2476771/150‑1 (MQ=255)
                                                                                                                                       aaGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGCACAg              >  2:3497546/1‑150 (MQ=255)
                                                                                                                                          ccccGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGCACAgtgg           >  2:630233/1‑150 (MQ=255)
                                                                                                                                           cccGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGCACAgtggt          <  1:1014805/150‑1 (MQ=255)
                                                                                                                                               gTTTTATTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGCACAGTGGTGttt      <  2:1229634/150‑1 (MQ=255)
                                                                                                                                                ttttCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGCACAGTGGTGTTTc     <  1:1631283/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTCAATCTGCTGCGCGCTCAATGTTACCGGGCTGTGGTTAATCACTCTGGACTTAACGCTGCCGTCTTCCAGAAGCACGTCAACTTCAAGCAACCCGTTGATATCATAACTAAAGCGAATATCAATCGACTGATAAGCCCCGGTTTTCTTCAACGGCACATCGAAGGATTCCACCAGAATGTTGTTTTTAACTTTGTGGTTTTCTCCCTGATAGACGTTAACCGTAATTGAATCCTGTTCCGGGTGCATGGTTGAATAAGTTTCTACACGCGACACGGGCACAGTGGTGTTTCGTT  >  CP009273/677409‑677705

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