Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,757,516 C→A 10.2% V333F (GTC→TTC)  yfjL ← CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,757,5160CA10.2% 730.5 / 45.8 205V333F (GTC→TTC) yfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base C (68/116);  new base A (11/10);  total (79/126)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 205 total aligned reads displayed.

TATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTACTGCT  >  CP009273/2757369‑2757662
                                                                                                                                                   |                                                                                                                                                  
tatGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACgg                                                                                                                                                  >  1:3088575/1‑150 (MQ=255)
tatGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACgg                                                                                                                                                  <  2:2760483/150‑1 (MQ=255)
 atGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGa                                                                                                                                                 >  2:737406/1‑150 (MQ=255)
 atGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGa                                                                                                                                                 >  2:1279327/1‑150 (MQ=255)
 atGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGa                                                                                                                                                 <  1:1552131/150‑1 (MQ=255)
  tGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGaa                                                                                                                                                <  1:3232781/150‑1 (MQ=255)
  tGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGaa                                                                                                                                                <  1:2270672/150‑1 (MQ=255)
  tGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTACTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGaa                                                                                                                                                <  1:1366709/150‑1 (MQ=255)
      cGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGGGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCt                                                                                                                                            <  1:1937/150‑1 (MQ=255)
           tCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGt                                                                                                                                       <  2:3126477/150‑1 (MQ=255)
             aaGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCg                                                                                                                                     >  1:721093/1‑150 (MQ=255)
             aaGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCg                                                                                                                                     >  2:1541278/1‑150 (MQ=255)
                aaTTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTtc                                                                                                                                  <  1:1837250/150‑1 (MQ=255)
                aaTTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTtc                                                                                                                                  <  1:2122579/150‑1 (MQ=255)
                aaTTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTtc                                                                                                                                  <  1:3257381/150‑1 (MQ=255)
                aaTTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTtc                                                                                                                                  <  1:1837049/150‑1 (MQ=255)
                 aTTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTtct                                                                                                                                 >  1:751140/1‑150 (MQ=255)
                 aTTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTtct                                                                                                                                 >  2:608969/1‑150 (MQ=255)
                     tttCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCt                                                                                                                             <  2:1780632/150‑1 (MQ=255)
                     tttCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCt                                                                                                                             <  1:703831/150‑1 (MQ=255)
                     tttCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGAAGGAATTCTTTAGTCGTTCTCGCt                                                                                                                             <  2:2227326/150‑1 (MQ=255)
                      ttCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCtt                                                                                                                            <  2:1305052/150‑1 (MQ=255)
                          gCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTct                                                                                                                        <  1:478490/150‑2 (MQ=255)
                          gCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCa                                                                                                                        <  1:1477791/150‑1 (MQ=255)
                          gCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCa                                                                                                                        >  1:417222/1‑150 (MQ=255)
                             ccTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCaatg                                                                                                                     <  2:4328678/150‑1 (MQ=255)
                                tcTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGAAGGAATTCTTTAGTCGTTCTCGCTTTTCaatgaat                                                                                                                  <  2:4020945/149‑1 (MQ=255)
                                gCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGAAGGAATTCTTTAGTCGTTCTCGCTTTTCaatgaat                                                                                                                  >  2:2818401/1‑150 (MQ=255)
                                 cTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCaatgaatg                                                                                                                 >  2:3581591/1‑150 (MQ=255)
                                   ttttCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTc                                                                                                               <  2:567898/150‑1 (MQ=255)
                                    tttCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCa                                                                                                              >  2:2232577/1‑150 (MQ=255)
                                    tttCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGAAGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCa                                                                                                              <  2:46467/150‑1 (MQ=255)
                                     ttCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAg                                                                                                             <  2:1644735/150‑1 (MQ=255)
                                       ccACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGtt                                                                                                           <  1:3048358/150‑1 (MQ=255)
                                           gTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATg                                                                                                       >  2:1927905/1‑150 (MQ=255)
                                            tataGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGACTGTCAGTTTATGa                                                                                                      >  2:2919355/1‑150 (MQ=255)
                                            tataGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGa                                                                                                      >  1:2526382/1‑150 (MQ=255)
                                            tataGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGa                                                                                                      >  2:2502179/1‑150 (MQ=255)
                                               aGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAata                                                                                                   <  2:3908772/150‑1 (MQ=255)
                                               aGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAata                                                                                                   <  1:1917320/150‑1 (MQ=255)
                                                gCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAatat                                                                                                  <  2:1862258/150‑1 (MQ=255)
                                                  aTGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCa                                                                                                >  1:1999034/1‑150 (MQ=255)
                                                  aTGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCa                                                                                                <  2:285324/150‑1 (MQ=255)
                                                    gCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGa                                                                                              <  1:544941/150‑1 (MQ=255)
                                                    gCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGa                                                                                              <  1:304037/150‑1 (MQ=255)
                                                    gCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGa                                                                                              <  1:1012245/150‑1 (MQ=255)
                                                      aaGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACa                                                                                            >  2:1097938/1‑150 (MQ=255)
                                                        ggAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACAtt                                                                                          >  1:4489228/1‑150 (MQ=255)
                                                         gAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACAttt                                                                                         <  1:3132025/150‑1 (MQ=255)
                                                           aTAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGAAGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCt                                                                                       >  2:4214415/1‑150 (MQ=255)
                                                            tAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCtt                                                                                      >  1:4417121/1‑150 (MQ=255)
                                                               gctggctTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTg                                                                                   <  1:42876/150‑1 (MQ=255)
                                                                ctggctTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTgg                                                                                  <  2:2578758/150‑1 (MQ=255)
                                                                ctggctTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTgg                                                                                  <  2:3270730/150‑1 (MQ=255)
                                                                 tggctTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGt                                                                                 >  2:776669/1‑150 (MQ=255)
                                                                  ggctTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCTATTTGAAGACGGAATTCTTTAGTCGTTCTCTCTTTTCAATGAATTTCATTTTATGAATATCAGACATTTCTTTTGTTc                                                                                <  2:2156339/150‑1 (MQ=255)
                                                                  ggctTTCGTGGTGCCGCCGGATCGTTGCGTGCACCAGCACTCGACCCAAAAAGTCGCTGACGTCATGGCGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTc                                                                                <  2:4726949/150‑1 (MQ=255)
                                                                    ctTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCgg                                                                              >  2:4527520/1‑150 (MQ=255)
                                                                      ttCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGgag                                                                            <  1:885376/150‑1 (MQ=255)
                                                                       tCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGgaga                                                                           <  1:1325577/150‑1 (MQ=255)
                                                                        cGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGaa                                                                          <  1:593078/150‑1 (MQ=255)
                                                                        cGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGaa                                                                          >  2:713782/1‑150 (MQ=255)
                                                                        cGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGaa                                                                          <  2:3721234/150‑1 (MQ=255)
                                                                        cGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGaa                                                                          >  2:1034122/1‑150 (MQ=255)
                                                                        cGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGaa                                                                          <  2:2338386/150‑1 (MQ=255)
                                                                         gtggtgCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGaaa                                                                         >  1:3397426/1‑150 (MQ=255)
                                                                         gtggtgCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGaaa                                                                         >  1:2412134/1‑150 (MQ=255)
                                                                          tggtgCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGaaaa                                                                        >  2:1872462/1‑150 (MQ=255)
                                                                               ctgcCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTc                                                                   >  1:4490368/1‑150 (MQ=255)
                                                                                 gcCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCAt                                                                 <  1:395285/150‑1 (MQ=255)
                                                                                       tCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTGTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAAcc                                                           >  1:3453271/1‑150 (MQ=255)
                                                                                             cGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGAAGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGAc                                                     <  2:655446/150‑1 (MQ=255)
                                                                                              gTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACa                                                    <  1:278959/150‑1 (MQ=255)
                                                                                                  aaCAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACg                                                >  1:278069/1‑150 (MQ=255)
                                                                                                       cACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATc                                           <  1:3601846/150‑1 (MQ=255)
                                                                                                        aCTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATcc                                          >  1:3865090/1‑150 (MQ=255)
                                                                                                            gACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGa                                      <  2:3330246/150‑1 (MQ=255)
                                                                                                              cccAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGAtg                                    >  2:1477601/1‑150 (MQ=255)
                                                                                                                   aaaTTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCtt                               <  2:751140/150‑1 (MQ=255)
                                                                                                                     aTTCGCTGACGTCATGGAGCGATTTGAAGAAGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGt                             >  1:101976/1‑150 (MQ=255)
                                                                                                                         gCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCg                         <  1:2852397/150‑1 (MQ=255)
                                                                                                                         gCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCg                         <  2:2070580/150‑1 (MQ=255)
                                                                                                                          cTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCgg                        <  1:4196235/150‑1 (MQ=255)
                                                                                                                           tGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGt                       >  1:3772894/1‑150 (MQ=255)
                                                                                                                            gACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTa                      <  1:256209/150‑1 (MQ=255)
                                                                                                                             aCGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTAt                     <  2:4489228/150‑1 (MQ=255)
                                                                                                                                 cATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTAAATGAAAGTCAGTTTATGAAAATAAGACATTTATATTGTTCGGAGAAAAACCTAATCCAACAATTTAAAGACGTTATCATCGATGTGCTTGTAGAGGTAtttat                 >  1:4041795/1‑148 (MQ=255)
                                                                                                                                    ggAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCATACATATCTTTTGGTCGGAGAAAAACCTCATCAAACCATTGACAGACGTTATACGCGATGTGCTTGTAGCGGTATTTCTTaa              >  1:4722761/1‑150 (MQ=255)
                                                                                                                                     gAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAg             <  1:879009/150‑1 (MQ=255)
                                                                                                                                       gCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTg           <  1:2117428/150‑1 (MQ=255)
                                                                                                                                       gCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTg           <  1:1596414/150‑1 (MQ=255)
                                                                                                                                        cGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGa          >  1:382661/1‑150 (MQ=255)
                                                                                                                                         gATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAg         <  2:3392429/150‑1 (MQ=255)
                                                                                                                                           tttGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTa       >  2:58442/1‑150 (MQ=255)
                                                                                                                                             tGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTAct     >  2:2217032/1‑150 (MQ=255)
                                                                                                                                             tGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTAct     >  1:870807/1‑150 (MQ=255)
                                                                                                                                              gAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTActg    <  1:2300386/150‑1 (MQ=255)
                                                                                                                                               aaGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTActgc   >  2:4327216/1‑150 (MQ=255)
                                                                                                                                                aGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTActgct  >  1:2184786/1‑150 (MQ=255)
                                                                                                                                                aGAAGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTActgct  >  1:2963902/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TATGCACGTTTTCAAGAATTTTTTCGGCCCCTGCTTTTTCCACGTATAGCATGCAAGGAATAGGCTGGCTTTCGTGGTGCTGCCGGATCGTTGCGTGCAACAGCACTCGACCCAAAAATTCGCTGACGTCATGGAGCGATTTGAAGACGGAATTCTTTAGTCGTTCTCGCTTTTCAATGAATGTCAGTTTATGAATATCAGACATTTCTTTTGGTCGGAGAAAAACCTCATCCAACCATTGACAGACGTTATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTACTGCT  >  CP009273/2757369‑2757662

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: