Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,024,003 C→A 15.7% D59Y (GAT→TAT)  hspQ ← heat shock protein involved in degradation of mutant DnaA; hemimethylated oriC DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,024,0030CA15.7% 801.3 / 88.1 278D59Y (GAT→TAT) hspQheat shock protein involved in degradation of mutant DnaA; hemimethylated oriC DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (137/97);  new base A (23/21);  total (160/118)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 281 total aligned reads displayed.

CAGGGATTAGTTACGCAGACGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTG  >  CP009273/1023854‑1024150
                                                                                                                                                     |                                                                                                                                                   
cAGGGATTAGTTACGCAGACGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATc                                                                                                                                                     <  2:738116/150‑1 (MQ=255)
 aGGGATTAGTTACGCAGACGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCg                                                                                                                                                    >  2:3076877/1‑150 (MQ=255)
 aGGGATTAGTTACGCAGACGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCg                                                                                                                                                    >  2:1764204/1‑150 (MQ=255)
  gggATTAGTTACGCAGACGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGTACCGGTAGGCCGTTATCGt                                                                                                                                                   >  1:4484469/1‑150 (MQ=255)
   ggATTAGTTACGCAGACGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATAGtc                                                                                                                                                  <  2:1139869/150‑1 (MQ=255)
    gATTAGTTACTCAGACGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCAACCATGGATGGCTGTTCAGGATGTACATCCTGCAAATCGCTACTTAGCTGGGCTACAGCCAGGTAGGTATGAACCGGTAGGCCGATATCGtcc                                                                                                                                                 >  2:713744/1‑150 (MQ=255)
    gATTAGTTACGCAGACGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGtcc                                                                                                                                                 >  1:4261763/1‑150 (MQ=255)
       tAGTTACGCAGACGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGtcctcc                                                                                                                                              <  1:921255/150‑1 (MQ=255)
        aGTTACGCAGACGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCa                                                                                                                                             <  2:332993/150‑1 (MQ=255)
            aCGCAGACGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGAAGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCAc                                                                                                                                         >  1:392198/1‑150 (MQ=255)
                  aCGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTg                                                                                                                                   <  2:1366266/150‑1 (MQ=255)
                  aCGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTg                                                                                                                                   <  1:3879262/150‑1 (MQ=255)
                  aCGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTg                                                                                                                                   <  1:3591526/150‑1 (MQ=255)
                      ggggCCAGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCAACCATGGATGGCTGTACAGGATGTTCAACCTGCAATACGCTAATTAGCTGGGCTACAGCCAGGTAGGTATGAAACGGAAGGCCGTTATCGACCACCAACACTAAGTGATaa                                                                                                                               >  1:2963999/1‑149 (MQ=255)
                       gggCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGCTGGCTGTTCAGGATGTTCATCCTGCAAGTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATAcc                                                                                                                              >  2:3250712/1‑150 (MQ=255)
                       gggCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATAcc                                                                                                                              <  1:3053999/150‑1 (MQ=255)
                       gggCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATAGTCCTCCATCACTACGTGATAcc                                                                                                                              >  1:3884124/1‑150 (MQ=255)
                         gCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATTGCTGTTCATGATGTTCATACTTCAATTCGCTACTTATCTGGGCTTAAGCCAGGTAGGTATGAAACGGTAGGCCGTTATAGTCCTCCATCACTAAGTGATAacac                                                                                                                            >  1:1409250/1‑150 (MQ=255)
                            tGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGa                                                                                                                         >  2:2979142/1‑150 (MQ=255)
                            tGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATACATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGa                                                                                                                         >  2:1903105/1‑150 (MQ=255)
                            tGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCAAACCAGGATGGCTGTTCAGGATGTTAATCCTGCAATTAGCTACATAGCTGGGCTTAAGCCACGTAGGTATGAACCGGTAGGCCGTTATAGTCCTCCAACACAACCTTATACCACGGa                                                                                                                         >  1:4946741/1‑150 (MQ=255)
                            tGGAGTTGTTTGCGGATGGTCTGCGCCAGCACATCCAAGGATGGCTGTTAAGGATGTTCATACTGAAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCAACACAACGTTATACCACGGa                                                                                                                         >  1:4052670/1‑150 (MQ=255)
                             ggAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAg                                                                                                                        >  2:5111609/1‑150 (MQ=255)
                             ggAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAg                                                                                                                        <  2:440138/150‑1 (MQ=255)
                               aGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAgcg                                                                                                                      >  1:3105680/1‑150 (MQ=255)
                                gttgttTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAgcgg                                                                                                                     >  1:1043515/1‑150 (MQ=255)
                                    ttttctgATGGTCTGCGCCAGATCATAAATGGATGGCTGTTCAGGATGTTCATACTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGGAGTCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGcgcg                                                                                                                 <  2:5238295/144‑1 (MQ=255)
                                    tttGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATAGTCCTCCATCACTACGTGATACCACGGAGCGGcgcg                                                                                                                 >  2:167947/1‑150 (MQ=255)
                                       gCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAg                                                                                                              >  2:3519857/1‑150 (MQ=255)
                                       gCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAg                                                                                                              <  2:4525384/150‑1 (MQ=255)
                                         ggATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCt                                                                                                            >  2:1865118/1‑150 (MQ=255)
                                         ggATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCt                                                                                                            >  1:4874236/1‑150 (MQ=255)
                                          gATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTc                                                                                                           >  2:3894819/1‑150 (MQ=255)
                                             ggTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTc                                                                                                        >  2:2308347/1‑150 (MQ=255)
                                              gTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACTGAGCGGCGCGAAGCTCGTCa                                                                                                       >  1:2727184/1‑150 (MQ=255)
                                                 tGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATAGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTc                                                                                                    >  1:1898625/1‑150 (MQ=255)
                                                   cgcCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCAc                                                                                                  <  2:4356829/150‑1 (MQ=255)
                                                      cAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGc                                                                                               >  2:3609158/1‑150 (MQ=255)
                                                       aGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGcc                                                                                              >  1:3655344/1‑150 (MQ=255)
                                                        gCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAGCCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCa                                                                                             <  1:3503453/150‑1 (MQ=255)
                                                          tCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGc                                                                                           >  1:2503408/1‑150 (MQ=255)
                                                           cATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCt                                                                                          <  1:4709833/150‑1 (MQ=255)
                                                            aTCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATAGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCtc                                                                                         >  1:112318/1‑150 (MQ=255)
                                                                 tggatggCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAg                                                                                    <  1:225398/150‑1 (MQ=255)
                                                                    atggCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCg                                                                                 >  1:1185673/1‑150 (MQ=255)
                                                                     tggCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGa                                                                                >  1:4313078/1‑150 (MQ=255)
                                                                          gTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGtt                                                                           <  2:2145834/150‑1 (MQ=255)
                                                                           ttCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTc                                                                          <  1:2023157/150‑1 (MQ=255)
                                                                            tCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCa                                                                         >  1:3072275/1‑150 (MQ=255)
                                                                               ggATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCagag                                                                      <  2:2241720/150‑1 (MQ=255)
                                                                               ggATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCagag                                                                      >  2:4725745/1‑150 (MQ=255)
                                                                               ggATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCagag                                                                      <  1:846056/150‑1 (MQ=255)
                                                                               ggATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCagag                                                                      <  1:813217/150‑1 (MQ=255)
                                                                                gATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCagaga                                                                     >  1:3176925/1‑150 (MQ=255)
                                                                                   gTTCATCCTGCTATTCGCTACTGAGCTGGGCTTCAGCCAGGTAGGTATGAATCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCg                                                                  <  2:4933565/150‑1 (MQ=255)
                                                                                   gTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGGATGAACCGGGAGGCCGTTATAGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCg                                                                  <  2:4834702/150‑1 (MQ=255)
                                                                                     tCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGaa                                                                <  2:1277880/150‑1 (MQ=255)
                                                                                        tCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATaa                                                             <  2:3596197/150‑1 (MQ=255)
                                                                                         ccTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATAGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATaaa                                                            <  1:1432773/150‑1 (MQ=255)
                                                                                            gCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATAGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCg                                                         >  1:2534686/1‑150 (MQ=255)
                                                                                             cAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCgg                                                        >  2:632307/1‑150 (MQ=255)
                                                                                               aTTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGt                                                      <  1:4069977/150‑1 (MQ=255)
                                                                                                gtCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTc                                                     <  2:4161084/149‑1 (MQ=255)
                                                                                                    cTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGata                                                 >  2:3430978/1‑150 (MQ=255)
                                                                                                    cTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGata                                                 >  1:523143/1‑150 (MQ=255)
                                                                                                     tACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGatat                                                <  2:1028711/150‑1 (MQ=255)
                                                                                                      aCTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATc                                               <  1:1082490/150‑1 (MQ=255)
                                                                                                      aCTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATc                                               >  1:210524/1‑150 (MQ=255)
                                                                                                      aCTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATc                                               >  1:2896462/1‑150 (MQ=255)
                                                                                                       cTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATcc                                              >  1:669104/1‑150 (MQ=255)
                                                                                                        ttAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCa                                             >  2:1527812/1‑150 (MQ=255)
                                                                                                         tAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATAGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCAc                                            <  2:1940922/150‑1 (MQ=255)
                                                                                                          aGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACg                                           <  1:1874327/150‑1 (MQ=255)
                                                                                                                  ttCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTcc                                   <  2:4807057/150‑1 (MQ=255)
                                                                                                                   tCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATAGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCg                                  <  2:2457397/150‑1 (MQ=255)
                                                                                                                       ccAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGAt                              >  2:4188195/1‑150 (MQ=255)
                                                                                                                       ccAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGAt                              <  2:3803921/150‑1 (MQ=255)
                                                                                                                        cAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGATGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATa                             >  1:1663465/1‑150 (MQ=255)
                                                                                                                        cAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATa                             >  1:188424/1‑150 (MQ=255)
                                                                                                                        cAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATa                             >  1:3323758/1‑150 (MQ=255)
                                                                                                                        cAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATa                             >  2:4108864/1‑150 (MQ=255)
                                                                                                                         aggtaggtATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATaa                            <  1:1472787/150‑1 (MQ=255)
                                                                                                                          ggtaggtaTGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAAc                           <  1:4633880/150‑1 (MQ=255)
                                                                                                                          ggtaggtaTGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAAc                           <  2:1669677/150‑1 (MQ=255)
                                                                                                                              ggtaTGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATCTCCACGACCACTCCGAGATAACCTaa                       >  2:2580631/1‑150 (MQ=255)
                                                                                                                              ggtaTGAACCGGTAGGCCGTTATAGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTaa                       >  1:300274/1‑150 (MQ=255)
                                                                                                                                taTGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACa                     >  1:559794/1‑150 (MQ=255)
                                                                                                                                 aTGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAg                    <  1:1764204/150‑1 (MQ=255)
                                                                                                                                   gAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAggg                  <  2:5317139/150‑1 (MQ=255)
                                                                                                                                    aaCCGGTAGGCCGTTATAGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGa                 >  1:3862840/1‑150 (MQ=255)
                                                                                                                                       cGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATg              >  2:2106131/1‑150 (MQ=255)
                                                                                                                                       cGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATg              >  2:2813999/1‑150 (MQ=255)
                                                                                                                                         gTAGGCCGTTATAGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGc            <  2:4871739/150‑1 (MQ=255)
                                                                                                                                          tAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCg           >  2:2730221/1‑150 (MQ=255)
                                                                                                                                           aGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCgg          >  2:405416/1‑150 (MQ=255)
                                                                                                                                            ggCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGa         <  2:1043515/150‑1 (MQ=255)
                                                                                                                                               cGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACct      <  2:5027202/150‑1 (MQ=255)
                                                                                                                                                gTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACctg     >  2:852992/1‑150 (MQ=255)
                                                                                                                                                 ttATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACctgc    >  1:3454200/1‑150 (MQ=255)
                                                                                                                                                  tataGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACctgct   >  1:3390111/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAGGGATTAGTTACGCAGACGCGGGGCCTGGAGTTGTTTGCGGATGGTCTGCGCCAGCTCATCCATGGATGGCTGTTCAGGATGTTCATCCTGCAATTCGCTACTTAGCTGGGCTTCAGCCAGGTAGGTATGAACCGGTAGGCCGTTATCGTCCTCCATCACTACGTGATACCACGGAGCGGCGCGAAGCTCGTCATTCACCGCCAGCTCATCAGGCGACGGTTCAGAGAGCGAATAAACCGGGTCGATATCCACGACCACTCCGAGATAACCTAACAGGGAATGGCGGACCTGCTG  >  CP009273/1023854‑1024150

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