Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,764,641 C→A 7.9% S45R (AGC→AGA yfjT → CP4‑57 prophage; putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,764,6410CA 7.9% 774.2 / 27.3 214S45R (AGC→AGAyfjTCP4‑57 prophage; putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (106/91);  new base A (7/10);  total (113/101)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 217 total aligned reads displayed.

AAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGATGGACCATAAGCAA  >  CP009273/2764494‑2764788
                                                                                                                                                   |                                                                                                                                                   
         aTAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATg                                                                                                                                            >  2:881424/1‑148 (MQ=255)
           aaaTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGtt                                                                                                                                          <  1:1123233/148‑1 (MQ=255)
           aaaTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGtt                                                                                                                                          >  1:245533/1‑148 (MQ=255)
           aaaTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGtt                                                                                                                                          >  1:245335/1‑148 (MQ=255)
             aTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTa                                                                                                                                        >  2:1062376/1‑148 (MQ=255)
             aTGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTa                                                                                                                                        >  2:2922812/1‑148 (MQ=255)
              tGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTAc                                                                                                                                       >  1:239441/1‑148 (MQ=255)
               gAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTAcc                                                                                                                                      <  2:923292/148‑1 (MQ=255)
                    tACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGAGGAAAGATGTTTACCATCaa                                                                                                                                 <  1:540215/148‑1 (MQ=255)
                      cGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACt                                                                                                                               >  2:1237197/1‑148 (MQ=255)
                        ttCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGc                                                                                                                             >  2:2417615/1‑148 (MQ=255)
                        ttCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGc                                                                                                                             >  2:1565740/1‑148 (MQ=255)
                         tCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCa                                                                                                                            <  1:268716/148‑1 (MQ=255)
                           aCTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACg                                                                                                                          <  2:490556/148‑1 (MQ=255)
                            cTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACgg                                                                                                                         >  2:2697495/1‑148 (MQ=255)
                              aaGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGa                                                                                                                       >  2:2129282/1‑148 (MQ=255)
                              aaGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGa                                                                                                                       >  2:2128589/1‑148 (MQ=255)
                               aGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGaa                                                                                                                      <  1:2504103/148‑1 (MQ=255)
                                gCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGAGGAAAGATGTTTACCATCAACTGCACGGGaaa                                                                                                                     >  2:1922368/1‑148 (MQ=255)
                                 ccGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAAt                                                                                                                    >  1:2153742/1‑148 (MQ=255)
                                 ccGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAAt                                                                                                                    <  2:466826/148‑1 (MQ=255)
                                 ccGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAAt                                                                                                                    <  2:779075/148‑1 (MQ=255)
                                    gTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGAGGAAAGATGTTTACCATCAACTGCACGGGAAATCCg                                                                                                                 <  1:3179333/148‑1 (MQ=255)
                                    gTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGAGGAAAGATGTTTACCATCAACTGCACGGGAAATCCg                                                                                                                 <  1:1067497/148‑1 (MQ=255)
                                     tttGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCgg                                                                                                                >  2:1268845/1‑148 (MQ=255)
                                      ttGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGa                                                                                                               >  1:2733020/1‑148 (MQ=255)
                                      ttGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGa                                                                                                               >  2:95907/1‑148 (MQ=255)
                                          aCTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAg                                                                                                           <  1:2268797/148‑1 (MQ=255)
                                           cTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGa                                                                                                          <  1:2995654/148‑1 (MQ=255)
                                            tGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGaa                                                                                                         <  1:2460612/148‑1 (MQ=255)
                                            tGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGaa                                                                                                         <  1:2341023/148‑1 (MQ=255)
                                                  gTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTAt                                                                                                   <  2:1438608/148‑1 (MQ=255)
                                                       aTGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTAcc                                                                                              <  2:173256/148‑1 (MQ=255)
                                                        tGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCa                                                                                             >  1:1559447/1‑148 (MQ=255)
                                                        tGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCa                                                                                             >  1:827622/1‑148 (MQ=255)
                                                         gTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGAGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAg                                                                                            >  1:1823746/1‑148 (MQ=255)
                                                           tCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCa                                                                                          <  2:52328/148‑1 (MQ=255)
                                                           tCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCa                                                                                          <  1:3067142/148‑1 (MQ=255)
                                                             gCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACt                                                                                        <  2:2583742/148‑1 (MQ=255)
                                                              ccAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGAGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTc                                                                                       <  2:883034/148‑1 (MQ=255)
                                                               cAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCg                                                                                      >  1:852179/1‑148 (MQ=255)
                                                                 gCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGt                                                                                    <  1:1932179/148‑1 (MQ=255)
                                                                 gCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGt                                                                                    <  1:262527/148‑1 (MQ=255)
                                                                    ttACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTAt                                                                                                   <  2:676569/130‑1 (MQ=255)
                                                                    ttACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTAt                                                                                                   >  1:676558/1‑130 (MQ=255)
                                                                    ttACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGtttt                                                                                 >  1:1299418/1‑148 (MQ=255)
                                                                     tACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGttttt                                                                                <  1:845801/148‑1 (MQ=255)
                                                                       cAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCt                                                                              >  2:602486/1‑148 (MQ=255)
                                                                          ccTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGAt                                                                           <  2:2478945/148‑1 (MQ=255)
                                                                          ccTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGAt                                                                           <  1:57159/148‑1 (MQ=255)
                                                                            tctGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGACAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATcc                                                                         >  2:2227198/1‑148 (MQ=255)
                                                                            tctGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATcc                                                                         >  2:3024758/1‑148 (MQ=255)
                                                                             ctGAGGTCCCCGTACTCTGGCAACAGGGTTACGGTCAGGGCAATGCGGAATACAGTGTGAAGGATGCCAGCGGAAAGATGATTACCAAAAACTGCACGGGAAATCCGGACCAGAACGGTATATAAAAGCACTAGGTTTATCTGATcat                                                                        >  2:1356929/1‑146 (MQ=21)
                                                                                cGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGc                                                                     >  2:2529390/1‑148 (MQ=255)
                                                                                  gTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCg                                                                   <  1:2031226/148‑1 (MQ=255)
                                                                                   tCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCgg                                                                  >  1:1357599/1‑148 (MQ=255)
                                                                                    ccccGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCggg                                                                 >  2:2171053/1‑148 (MQ=255)
                                                                                    ccccGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCggg                                                                 >  1:1398089/1‑148 (MQ=255)
                                                                                    ccccGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCggg                                                                 <  1:1705278/148‑1 (MQ=255)
                                                                                     cccGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCgggg                                                                <  1:2543446/148‑1 (MQ=255)
                                                                                     cccGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCgggg                                                                <  2:1953084/148‑1 (MQ=255)
                                                                                      ccGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGa                                                               <  1:1536145/148‑1 (MQ=255)
                                                                                       cGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGAt                                                              <  1:764232/148‑1 (MQ=255)
                                                                                             ctGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACg                                                        <  2:1722515/148‑1 (MQ=255)
                                                                                             ctGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGAGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACg                                                        <  1:396246/148‑1 (MQ=255)
                                                                                              tGGCAACAGGGGTACGGTCAGGGCAATGCGGAATATAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACgg                                                       <  1:3277716/148‑1 (MQ=255)
                                                                                              tGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACgg                                                       <  1:136139/148‑1 (MQ=255)
                                                                                              tGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACgg                                                       <  1:2156921/148‑1 (MQ=255)
                                                                                                  aaCAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAg                                                   >  2:66331/1‑148 (MQ=255)
                                                                                                  aaCAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAg                                                   <  1:3197080/148‑1 (MQ=255)
                                                                                                   aCAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGc                                                  >  1:884964/1‑148 (MQ=255)
                                                                                                   aCAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGc                                                  >  1:640041/1‑148 (MQ=255)
                                                                                                   aCAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGc                                                  >  1:2540749/1‑148 (MQ=255)
                                                                                                    cAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCt                                                 <  1:2845921/148‑1 (MQ=255)
                                                                                                       gggTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGc                                              >  1:2658918/1‑148 (MQ=255)
                                                                                                          tACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCatg                                           >  1:1965456/1‑148 (MQ=255)
                                                                                                          tACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCatg                                           >  1:2094131/1‑148 (MQ=255)
                                                                                                                  gggCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGt                                   >  2:113319/1‑148 (MQ=255)
                                                                                                                    gCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTAc                                 <  2:2745100/148‑1 (MQ=255)
                                                                                                                        tGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGc                             >  2:2825126/1‑148 (MQ=255)
                                                                                                                        tGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGc                             >  1:1219390/1‑148 (MQ=255)
                                                                                                                        tGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGc                             >  2:1166760/1‑148 (MQ=255)
                                                                                                                          cGGAATACAGTGTGACGGATGCCAGAGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCAt                           <  2:1323442/148‑1 (MQ=255)
                                                                                                                                tctgtgtGACGGATGCGTGCTGAAATATTTTTAACATCATGTGCACTGGAAATGAGTACATGTTTGGTATTTACGATGACTGGGTTTTTCTTATCCTCGCAGGTGTTAAAACGGTCAGCTCGCATGATGTCAGTAGCGGCTTcacttt                     <  2:3085555/145‑4 (MQ=255)
                                                                                                                                aCAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAgt                     >  2:2151068/1‑148 (MQ=255)
                                                                                                                                  aGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAgtgg                   <  2:1833271/148‑1 (MQ=255)
                                                                                                                                   gtgtGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAgtggt                  <  1:406340/148‑1 (MQ=255)
                                                                                                                                   gtgtGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAgtggt                  >  2:2318200/1‑148 (MQ=255)
                                                                                                                                    tgtgACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAgtggtg                 <  1:431869/148‑1 (MQ=255)
                                                                                                                                    tgtgACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAgtggtg                 >  2:1861290/1‑148 (MQ=255)
                                                                                                                                    tgtgACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAgtggtg                 >  2:1790031/1‑148 (MQ=255)
                                                                                                                                      tgACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGAt               <  2:149203/148‑1 (MQ=255)
                                                                                                                                      tgACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGAt               <  1:2052023/148‑1 (MQ=255)
                                                                                                                                      tgACGGATGCCAGAGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGAt               >  1:187900/1‑148 (MQ=255)
                                                                                                                                           gATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGATGGAcc          >  2:1286009/1‑148 (MQ=255)
                                                                                                                                              gCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGATGGACCATa       <  2:2050408/148‑1 (MQ=255)
                                                                                                                                               ccAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGATGGACCATaa      >  1:2276291/1‑148 (MQ=255)
                                                                                                                                                 aGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGATGGACCATAAGc    <  1:1453565/148‑1 (MQ=255)
                                                                                                                                                  gCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGATGGACCATAAGCa   >  1:257933/1‑148 (MQ=255)
                                                                                                                                                   cGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGATGGACCATAAGCaa  <  2:1096186/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAAGGAGACATAAATGAAAATACGTTCACTAAGCCGGTTTGTACTGGCAAGTACAATGTTCGCCAGCTTTACAGCCTCTGCGGTCCCCGGACTCTGGCAACAGGGGTACGGTCAGGGCAATGCGGAATACAGTGTGACGGATGCCAGCGGAAAGATGTTTACCATCAACTGCACGGGAAATCCGGACCAGAACGGTATTTACCAGCACTCGGTTTTTCTGATCCTCGCCGGGGATAAAACGGTCAGCTCGCATGATGACAGTACCGGCATCACAGTGGTGATGGACCATAAGCAA  >  CP009273/2764494‑2764788

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: