Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,063,736 T→C 45.9% intergenic (+26/+231) ymdF → / ← rutG KGG family protein/pyrimidine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,063,7360TC45.9% 52.9 / inf 305intergenic (+26/+231)ymdF/rutGKGG family protein/pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base T (88/77);  new base C (82/58);  total (170/135)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 311 total aligned reads displayed.

AGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCG  >  CP009273/1063587‑1063881
                                                                                                                                                     |                                                                                                                                                 
 gCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAttt                                                                                                                                                  >  1:374783/1‑150 (MQ=255)
    tCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTg                                                                                                                                               >  2:1635374/1‑150 (MQ=255)
    tCAGAAGCAGGTAAAAAAGGTGGACACCACATCGGGGGTAATATCAAGAATCACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTCGTAAGAGCAGTCACCGCAAAACCGACAACTAGCCGGGCTAAACAATCACGAATGCATTATTg                                                                                                                                               >  2:1141352/1‑150 (MQ=255)
     cAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGATCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGt                                                                                                                                              <  2:1612742/150‑1 (MQ=255)
      aGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTc                                                                                                                                             >  1:1091438/1‑150 (MQ=255)
      aGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTc                                                                                                                                             >  1:946174/1‑150 (MQ=255)
      aGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTc                                                                                                                                             >  1:4473920/1‑150 (MQ=255)
       gAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCt                                                                                                                                            >  1:104978/1‑150 (MQ=255)
            aGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGc                                                                                                                                       <  1:4415573/150‑1 (MQ=255)
            aGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGc                                                                                                                                       <  2:323386/150‑1 (MQ=255)
            aGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGc                                                                                                                                       >  2:430153/1‑150 (MQ=255)
            aGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGc                                                                                                                                       >  2:4203519/1‑150 (MQ=255)
              gTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTc                                                                                                                                     <  2:520550/150‑1 (MQ=255)
              gTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTc                                                                                                                                     <  2:4251730/150‑1 (MQ=255)
               tAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCg                                                                                                                                    >  2:1544628/1‑150 (MQ=255)
                aaaaaaGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGt                                                                                                                                   >  1:555374/1‑150 (MQ=255)
                 aaaaaGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTc                                                                                                                                  >  1:4600968/1‑150 (MQ=255)
                 aaaaaGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTc                                                                                                                                  >  2:4367151/1‑150 (MQ=255)
                   aaaGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCaa                                                                                                                                <  1:2097135/150‑1 (MQ=255)
                    aaGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCaaa                                                                                                                               >  2:2662651/1‑150 (MQ=255)
                      ggTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCaaaaa                                                                                                                             <  2:4538766/150‑1 (MQ=255)
                      ggTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCaaaaa                                                                                                                             >  1:643970/1‑150 (MQ=255)
                       gTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAg                                                                                                                            <  1:2257034/150‑1 (MQ=255)
                        tGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGc                                                                                                                           >  1:4520201/1‑150 (MQ=255)
                        tGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGc                                                                                                                           >  1:1496790/1‑150 (MQ=255)
                        tGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGc                                                                                                                           >  1:3470932/1‑150 (MQ=255)
                          gACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCa                                                                                                                         >  2:615536/1‑150 (MQ=255)
                             agcacagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCa                                                                                                                      >  1:1806558/1‑150 (MQ=255)
                               cacagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCAcc                                                                                                                    <  1:2652652/150‑1 (MQ=255)
                                  agcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCAccgcc                                                                                                                 >  1:2811643/1‑150 (MQ=255)
                                  agcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCAccgcc                                                                                                                 <  2:3466379/150‑1 (MQ=255)
                                    cGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCgg                                                                                                               >  2:2308522/1‑150 (MQ=255)
                                     gggggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGt                                                                                                              <  1:4188001/150‑1 (MQ=255)
                                       gggTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTa                                                                                                            >  1:3290822/1‑150 (MQ=255)
                                         gTAATTTCAAGAATGACCCGCCGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                                                          >  1:1883514/1‑150 (MQ=255)
                                         gTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTAcc                                                                                                          >  1:1359968/1‑150 (MQ=255)
                                           aaTTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCg                                                                                                        <  1:4058528/150‑1 (MQ=255)
                                                cAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCggtggt                                                                                                   <  2:1794820/150‑1 (MQ=255)
                                                 aaGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGtt                                                                                                  <  2:1244888/150‑1 (MQ=255)
                                                   gAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGa                                                                                                >  2:352630/1‑150 (MQ=255)
                                                     aTGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATa                                                                                              >  1:442085/1‑150 (MQ=255)
                                                      tGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATAc                                                                                             >  2:405759/1‑150 (MQ=255)
                                                         cccGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtga                                                                                          <  2:4516733/150‑1 (MQ=255)
                                                           cGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACtgatg                                                                                        >  1:983316/1‑150 (MQ=255)
                                                                   gcATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGt                                                                                >  1:2346119/1‑150 (MQ=255)
                                                                     aTCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTaa                                                                              >  2:1503237/1‑150 (MQ=255)
                                                                     aTCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTaa                                                                              >  1:146744/1‑150 (MQ=255)
                                                                      tCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACGAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAg                                                                             <  2:3311/150‑1 (MQ=255)
                                                                       cTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGc                                                                            >  2:3589033/1‑150 (MQ=255)
                                                                          aaGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTg                                                                         >  1:1030552/1‑150 (MQ=255)
                                                                            gCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGcc                                                                       >  2:2713482/1‑150 (MQ=255)
                                                                             cAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCt                                                                      <  1:1911942/150‑1 (MQ=255)
                                                                             cAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCt                                                                      <  2:1026695/150‑1 (MQ=255)
                                                                              aGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTg                                                                     >  2:4097401/1‑150 (MQ=255)
                                                                               ggTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGa                                                                    <  1:2091453/150‑1 (MQ=255)
                                                                                 tAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATg                                                                  <  2:2549872/150‑1 (MQ=255)
                                                                                 tAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATg                                                                  <  1:4457200/150‑1 (MQ=255)
                                                                                   aaaaaGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcg                                                                >  2:4091858/1‑150 (MQ=255)
                                                                                    aaaaGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATgcgc                                                               >  1:3405493/1‑150 (MQ=255)
                                                                                     aaaGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGcgcg                                                              >  2:4588962/1‑150 (MQ=255)
                                                                                     aaaGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGcgcg                                                              >  1:3984261/1‑150 (MQ=255)
                                                                                       aGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGAt                                                            <  2:3069916/150‑1 (MQ=255)
                                                                                       aGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGAt                                                            <  1:1552830/150‑1 (MQ=255)
                                                                                        ggtggtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATg                                                           >  1:2470125/1‑150 (MQ=255)
                                                                                          tggtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCt                                                         >  1:1570502/1‑150 (MQ=255)
                                                                                          tggtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCt                                                         >  1:4302849/1‑150 (MQ=255)
                                                                                           ggtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCtt                                                        >  1:3574383/1‑150 (MQ=255)
                                                                                            gtAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATTTAAGCTTGCCTGATGCGCGATGCTTa                                                       <  1:1398593/150‑1 (MQ=255)
                                                                                               agagCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCa                                                    <  2:585913/150‑1 (MQ=255)
                                                                                                gagCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTAGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAg                                                   >  2:803025/1‑150 (MQ=255)
                                                                                                    aGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCt                                               >  1:500782/1‑150 (MQ=255)
                                                                                                        aCGGCACAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCa                                           <  2:3132109/150‑1 (MQ=255)
                                                                                                          ggCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCaga                                         >  1:436610/1‑150 (MQ=255)
                                                                                                           gCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCagaa                                        <  1:1198648/150‑1 (MQ=255)
                                                                                                            cAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCagaag                                       >  2:1934003/1‑150 (MQ=255)
                                                                                                             aaaaGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCagaaga                                      >  2:1360923/1‑150 (MQ=255)
                                                                                                             aaaaGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCagaaga                                      >  2:3681018/1‑150 (MQ=255)
                                                                                                             aaaaGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCagaaga                                      >  1:1166754/1‑150 (MQ=255)
                                                                                                              aaaGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGAt                                     >  2:2490079/1‑150 (MQ=255)
                                                                                                               aaGCGACAACTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGAtt                                    <  1:3161186/150‑1 (MQ=255)
                                                                                                                    aCAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAat                               <  2:3532310/150‑1 (MQ=255)
                                                                                                                        cTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAatatat                           >  2:219695/1‑150 (MQ=255)
                                                                                                                        cTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAatatat                           <  2:3953039/150‑1 (MQ=255)
                                                                                                                        cTAGCCGGGCTAATCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAatatat                           <  1:2612917/150‑1 (MQ=255)
                                                                                                                         tAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATAtt                          <  1:2750021/150‑1 (MQ=255)
                                                                                                                             cGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAAt                      <  1:717730/150‑1 (MQ=255)
                                                                                                                              gggCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAAtt                     <  1:1313823/150‑1 (MQ=255)
                                                                                                                               ggCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAAttt                    <  2:2136434/150‑1 (MQ=255)
                                                                                                                                gCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTg                   >  1:2880467/1‑150 (MQ=255)
                                                                                                                                   aaTCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCAc                >  2:3353146/1‑150 (MQ=255)
                                                                                                                                    aTCAATGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACt               <  2:3027078/150‑1 (MQ=255)
                                                                                                                                        aTGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGttt           >  1:2276176/1‑150 (MQ=255)
                                                                                                                                        aTGACGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGttt           <  2:2034579/150‑1 (MQ=255)
                                                                                                                                         tGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCAACACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACttattt          >  1:2133870/1‑145 (MQ=255)
                                                                                                                                           aCGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGt        >  1:2966799/1‑150 (MQ=255)
                                                                                                                                           aCGAATGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGt        >  2:744668/1‑150 (MQ=255)
                                                                                                                                            cGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTa       <  2:294573/150‑1 (MQ=255)
                                                                                                                                              aaTGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATCTTGAATTTGCACTGTTTTGTAgg     >  2:3786572/1‑150 (MQ=255)
                                                                                                                                              aaTGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAgg     >  2:921214/1‑150 (MQ=255)
                                                                                                                                               aTGCATCTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGc    >  1:1703894/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCATCACCGCCGGTTACCCGGTGGTTGATACTGATGACAAATGTAAGCTTGCCTGATGCGCGATGCTTATCAGGCCTACCAGAAGATTGCAATATATTGAATTTGCACTGTTTTGTAGGCCG  >  CP009273/1063587‑1063881

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: