Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,591,139 C→A 47.6% D66Y (GAT→TAT)  dnaT ← DNA biosynthesis protein (primosomal protein I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,591,1390CA47.6% 151.6 / inf 306D66Y (GAT→TAT) dnaTDNA biosynthesis protein (primosomal protein I)
Reads supporting (aligned to +/- strand):  ref base C (76/83);  new base A (75/71);  total (151/155)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

CGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTAC  >  CP009273/4590990‑4591286
                                                                                                                                                     |                                                                                                                                                   
cGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATc                                                                                                                                                     >  1:3627990/1‑150 (MQ=255)
 ggCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCg                                                                                                                                                    >  2:3772745/1‑150 (MQ=255)
  gCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGc                                                                                                                                                   <  1:4325050/150‑1 (MQ=255)
    tctctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCtt                                                                                                                                                 >  2:3655243/1‑150 (MQ=255)
       ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTccc                                                                                                                                              <  2:433473/150‑1 (MQ=255)
       ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTccc                                                                                                                                              <  1:1883937/150‑1 (MQ=255)
       ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTccc                                                                                                                                              >  1:1517026/1‑150 (MQ=255)
       ctTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTccc                                                                                                                                              >  1:2510632/1‑150 (MQ=255)
         tAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCgg                                                                                                                                            >  1:2576507/1‑150 (MQ=255)
         tAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCgg                                                                                                                                            >  1:2426589/1‑150 (MQ=255)
         tAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCgg                                                                                                                                            >  1:3459229/1‑150 (MQ=255)
           gcgcCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGAc                                                                                                                                          <  1:2580522/150‑1 (MQ=255)
                aCGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGcc                                                                                                                                     >  1:541444/1‑150 (MQ=255)
                 cGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCa                                                                                                                                    <  1:81252/150‑1 (MQ=255)
                  gCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAg                                                                                                                                   <  2:3192215/150‑1 (MQ=255)
                    cccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCt                                                                                                                                 >  1:2929515/1‑150 (MQ=255)
                    cccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCt                                                                                                                                 <  1:2369538/150‑1 (MQ=255)
                    cccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCt                                                                                                                                 >  1:2843321/1‑150 (MQ=255)
                    cccATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCt                                                                                                                                 <  1:196353/150‑1 (MQ=255)
                        tAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTttc                                                                                                                             >  2:3370565/1‑150 (MQ=255)
                            gccgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTcc                                                                                                                         <  2:3844217/150‑1 (MQ=255)
                              cgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAg                                                                                                                       <  2:1826617/150‑1 (MQ=255)
                              cgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAg                                                                                                                       <  2:344657/150‑1 (MQ=255)
                              cgccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAg                                                                                                                       <  1:3136712/150‑1 (MQ=255)
                                ccAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAgcg                                                                                                                     <  1:2253902/150‑1 (MQ=255)
                                     cGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCcagc                                                                                                                <  1:1099711/150‑1 (MQ=255)
                                       gATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCcagcag                                                                                                              >  2:4024811/1‑150 (MQ=255)
                                        aTAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGt                                                                                                             <  2:1287724/150‑1 (MQ=255)
                                         tAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGtt                                                                                                            <  2:605321/150‑1 (MQ=255)
                                          aaaaTCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTc                                                                                                           >  2:1080636/1‑150 (MQ=255)
                                               cGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCca                                                                                                      <  2:3128696/150‑1 (MQ=255)
                                               cGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCca                                                                                                      <  1:3031363/150‑1 (MQ=255)
                                               cGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCca                                                                                                      <  2:500618/150‑1 (MQ=255)
                                                 gCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcagg                                                                                                    <  1:1664302/150‑1 (MQ=255)
                                                  cATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggc                                                                                                   <  2:928509/150‑1 (MQ=255)
                                                   aTCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcg                                                                                                  <  2:1333004/150‑1 (MQ=255)
                                                   aTCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcg                                                                                                  >  2:2736947/1‑150 (MQ=255)
                                                    tCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgt                                                                                                 >  1:1815930/1‑150 (MQ=255)
                                                     cGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtg                                                                                                <  2:2459866/150‑1 (MQ=255)
                                                       ggTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgca                                                                                              <  2:2509246/150‑1 (MQ=255)
                                                          tGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcaggc                                                                                           >  2:2840091/1‑150 (MQ=255)
                                                           gCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcaggcg                                                                                          >  1:4042830/1‑150 (MQ=255)
                                                             cAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcaggcgtg                                                                                        <  1:1294005/150‑1 (MQ=255)
                                                             cAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCcaggcgtgcaggcgtg                                                                                        <  2:846241/150‑1 (MQ=255)
                                                                    gATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCAt                                                                                 <  1:409165/150‑1 (MQ=255)
                                                                     aTACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATa                                                                                <  2:1941967/150‑1 (MQ=255)
                                                                     aTACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATa                                                                                <  1:2529022/150‑1 (MQ=255)
                                                                       aCATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATaaa                                                                              >  1:682140/1‑150 (MQ=255)
                                                                           ggCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGc                                                                          >  2:483947/1‑150 (MQ=255)
                                                                           ggCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGc                                                                          <  2:1879391/150‑1 (MQ=255)
                                                                            gCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGcc                                                                         <  1:1616118/150‑1 (MQ=255)
                                                                             cGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCg                                                                        <  1:1047735/150‑1 (MQ=255)
                                                                             cGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCg                                                                        >  2:469954/1‑150 (MQ=255)
                                                                             cGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCg                                                                        <  1:889696/150‑1 (MQ=255)
                                                                              gAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCgg                                                                       <  1:3124505/150‑1 (MQ=255)
                                                                                 tttCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGc                                                                    >  1:608086/1‑150 (MQ=255)
                                                                                   tCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCAt                                                                  >  1:300937/1‑150 (MQ=255)
                                                                                    ccccATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCAtt                                                                 >  1:2185034/1‑150 (MQ=255)
                                                                                     cccATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTg                                                                <  2:3593004/150‑1 (MQ=255)
                                                                                      ccATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGt                                                               <  2:3159472/150‑1 (MQ=255)
                                                                                       cATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGtt                                                              >  2:418457/1‑150 (MQ=255)
                                                                                          gggTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGAATAAAACGCCGGGGCATTGTTAGc                                                           <  1:4183966/150‑1 (MQ=255)
                                                                                             tACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTTACGGCATAAAACGCCGGGGCATTGTTAGCaaa                                                        >  1:849379/1‑150 (MQ=255)
                                                                                             tACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCaaa                                                        <  1:1707823/150‑1 (MQ=255)
                                                                                             tACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCaaa                                                        <  2:1517026/150‑1 (MQ=255)
                                                                                              aCAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAAt                                                       <  1:3365023/150‑1 (MQ=255)
                                                                                              aCAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGGTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCACAt                                                       >  2:3873286/1‑150 (MQ=255)
                                                                                                aGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATAc                                                     >  2:2318338/1‑150 (MQ=255)
                                                                                                  cAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACgg                                                   >  2:1622603/1‑150 (MQ=255)
                                                                                                   aaCGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGc                                                  >  2:1130879/1‑150 (MQ=255)
                                                                                                   aaCGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGc                                                  >  1:1308155/1‑150 (MQ=255)
                                                                                                    aCGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCCGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCAGTGTTAGCAAATACGGCa                                                 >  2:2901221/1‑150 (MQ=255)
                                                                                                     cGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCaa                                                >  1:496503/1‑150 (MQ=255)
                                                                                                         aGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACcac                                            >  1:3915211/1‑150 (MQ=255)
                                                                                                            ggCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACAcc                                         >  1:323750/1‑150 (MQ=255)
                                                                                                             gCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCg                                        <  2:2360321/150‑1 (MQ=255)
                                                                                                                 gCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCtt                                    <  1:1188087/150‑1 (MQ=255)
                                                                                                                  cGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTc                                   <  1:1379783/150‑1 (MQ=255)
                                                                                                                   ggTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCa                                  <  2:855351/150‑1 (MQ=255)
                                                                                                                   ggTTCCTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCa                                  <  2:1240367/150‑1 (MQ=255)
                                                                                                                     ttCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGc                                >  2:622896/1‑150 (MQ=255)
                                                                                                                       ccTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCtt                              >  2:3497718/1‑150 (MQ=255)
                                                                                                                       ccTGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCtt                              >  2:2762869/1‑150 (MQ=255)
                                                                                                                        cTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCttt                             <  1:2770689/150‑1 (MQ=255)
                                                                                                                        cTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCttt                             >  1:2396684/1‑150 (MQ=255)
                                                                                                                         tGATAGAGTTGATCGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCtttt                            >  1:3766118/1‑150 (MQ=255)
                                                                                                                              gagTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAg                       <  2:669221/150‑1 (MQ=255)
                                                                                                                                  tGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACg                   <  2:647360/150‑1 (MQ=255)
                                                                                                                                  tGAGCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACg                   <  1:1905336/150‑1 (MQ=255)
                                                                                                                                      cGTCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGttt               <  2:2219335/150‑1 (MQ=255)
                                                                                                                                        tCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTgg             <  2:1960234/150‑1 (MQ=255)
                                                                                                                                        tCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTgg             >  1:3841216/1‑150 (MQ=255)
                                                                                                                                        tCCAGAGCGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTgg             >  2:3238289/1‑150 (MQ=255)
                                                                                                                                          cAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTg           >  1:2090964/1‑150 (MQ=255)
                                                                                                                                           agagCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGa          >  1:3943570/1‑150 (MQ=255)
                                                                                                                                           agagCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGa          >  2:2150686/1‑150 (MQ=255)
                                                                                                                                           agagCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGa          >  1:3975466/1‑150 (MQ=255)
                                                                                                                                               cGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgt      >  2:254502/1‑150 (MQ=255)
                                                                                                                                               cGACATAGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgt      >  2:813413/1‑150 (MQ=255)
                                                                                                                                                gACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCgtg     >  2:2676451/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGGCTCTCTTAGCGCCACGCCCCATAGCGCCGCCAGGCGGATAAAATCGGCATCGGGTTGCCAGTCCGGATACATGGCGAATTTCCCCATGGGTACAGCAACGGGAGCGGCTTGCGGTTCCTGATAGAGTTGATCGTCCAGAGCGACATCGCTTCCCGGACGCGCCAGCTTTTCTTCCAGCGCCAGCAGTTCAGCCAGGCGTGCAGGCGTGACGGCATAAAACGCCGGGGCATTGTTAGCAAATACGGCAACCACACCGCCTTCAGCTTTTGCCAGAACGGTTTGGTGATCGTGTAC  >  CP009273/4590990‑4591286

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: