Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,960,118:1 +C 100% coding (1676/2247 nt) ptsP ← fused PTS enzyme: PEP‑protein phosphotransferase (enzyme I)/GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,960,1181.C96.8% 228.1 / 0.8 62coding (1676/2247 nt)ptsPfused PTS enzyme: PEP‑protein phosphotransferase (enzyme I)/GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (26/34);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGG  >  CP009273/2959970‑2960260
                                                                                                                                                     |                                                                                                                                              
cAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGAc                                                                                                                                               <  2:359753/150‑1 (MQ=255)
cAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGAc                                                                                                                                               <  1:472539/150‑1 (MQ=255)
  gATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTTCTTCATCAACTTCATCGAGGCTTGTGAccc                                                                                                                                             <  2:701198/150‑2 (MQ=255)
       ggCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCgg                                                                                                                                         >  2:240316/1‑150 (MQ=255)
       ggCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCgg                                                                                                                                         >  2:856637/1‑150 (MQ=255)
            caTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACa                                                                                                                                    <  2:861056/150‑1 (MQ=255)
            caTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACa                                                                                                                                    <  2:444032/150‑1 (MQ=255)
                 aTACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAata                                                                                                                               <  1:954705/150‑1 (MQ=255)
                 aTACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAata                                                                                                                               <  2:200290/150‑1 (MQ=255)
                 aTACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAata                                                                                                                               <  2:913716/150‑1 (MQ=255)
                  tACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAatat                                                                                                                              <  2:471577/150‑1 (MQ=255)
                     cATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCa                                                                                                                           <  2:362401/150‑1 (MQ=255)
                       ttGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAgg                                                                                                                         >  1:286212/1‑150 (MQ=255)
                         gACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGtt                                                                                                                       >  2:522596/1‑150 (MQ=255)
                           cGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGc                                                                                                                     <  1:16490/150‑1 (MQ=255)
                            ggCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGcc                                                                                                                    >  2:683103/1‑150 (MQ=255)
                             gCACTTCCAGCATGATGCCGGTACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGccc                                                                                                                   <  1:522683/150‑1 (MQ=255)
                             gCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGccc                                                                                                                   <  1:699804/150‑1 (MQ=255)
                                   ccAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGc                                                                                                              <  2:149940/149‑1 (MQ=255)
                                   ccAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCg                                                                                                             >  2:144601/1‑150 (MQ=255)
                                    cAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCg                                                                                                             <  2:813838/149‑1 (MQ=255)
                                       cATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCAt                                                                                                         >  2:789171/1‑150 (MQ=255)
                                         tgatgCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTa                                                                                                       <  1:686367/150‑1 (MQ=255)
                                           atgCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAg                                                                                                      <  2:560127/149‑1 (MQ=255)
                                           atgCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGc                                                                                                     <  2:310597/150‑1 (MQ=255)
                                              ccGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACg                                                                                                  >  2:179355/1‑150 (MQ=255)
                                                  tACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACgcagc                                                                                              <  2:212978/150‑1 (MQ=255)
                                                    cGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCAt                                                                                            >  2:257092/1‑150 (MQ=255)
                                                     gTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATc                                                                                           >  2:689411/1‑150 (MQ=255)
                                                         tttGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGccc                                                                                       >  1:930413/1‑150 (MQ=255)
                                                             ggAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCa                                                                                    <  1:977574/149‑1 (MQ=255)
                                                                   tCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGa                                                                              <  1:247922/149‑1 (MQ=255)
                                                                    cGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGAt                                                                             <  2:204416/149‑1 (MQ=255)
                                                                     gTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATc                                                                            >  2:793730/1‑149 (MQ=255)
                                                                     gTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCa                                                                           <  1:522596/150‑1 (MQ=255)
                                                                           gATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCaag                                                                         >  2:40827/1‑146 (MQ=255)
                                                                             tcatcTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGAtc                                                                   <  1:769084/150‑1 (MQ=255)
                                                                               atcTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGAtctc                                                                 >  1:864226/1‑150 (MQ=255)
                                                                               atcTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGAtctc                                                                 >  1:579631/1‑150 (MQ=255)
                                                                                    ctcGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCt                                                            <  2:259033/150‑1 (MQ=255)
                                                                                    ctcGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCg                                                               <  2:773078/146‑1 (MQ=255)
                                                                                        aCTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATc                                                        <  1:789056/150‑1 (MQ=255)
                                                                                          ttCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATc                                                        >  2:141280/1‑148 (MQ=255)
                                                                                          ttCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATc                                                        >  2:833207/1‑148 (MQ=255)
                                                                                                    gCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCg                                            <  1:240316/150‑1 (MQ=255)
                                                                                                     cACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCg                                            >  2:1015025/1‑149 (MQ=255)
                                                                                                        gTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATc                                        >  1:957352/1‑150 (MQ=255)
                                                                                                           cAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATccc                                      <  1:257092/149‑1 (MQ=255)
                                                                                                           cAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCa                                     >  2:207084/1‑150 (MQ=255)
                                                                                                            aaTCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCAc                                    <  1:940344/150‑1 (MQ=255)
                                                                                                             aTCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCAAGCGGCATGGCACGCGCCTGGATCAAGAAGATCTCCGGCTGCTCGGGCGTAATGGGAATCCCACg                                   <  2:408727/150‑1 (MQ=255)
                                                                                                              tCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGc                                  >  1:159565/1‑150 (MQ=255)
                                                                                                                aGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGcc                                 >  1:131105/1‑149 (MQ=255)
                                                                                                                 ggcggcGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCaa                               <  1:587945/150‑1 (MQ=255)
                                                                                                                       gTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAgg                         <  1:844570/150‑1 (MQ=255)
                                                                                                                           ttCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATg                     >  1:658304/1‑150 (MQ=255)
                                                                                                                           ttCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATg                     >  1:991241/1‑150 (MQ=255)
                                                                                                                             catcaACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATgg                    <  1:856637/149‑1 (MQ=255)
                                                                                                                              atcaACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGAtt                 <  2:705844/150‑1 (MQ=255)
                                                                                                                               tcaACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGAtt                 <  2:145139/150‑1 (MQ=255)
                                                                                                                                 aaCTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATtct               <  2:427567/150‑1 (MQ=255)
                                                                                                                                  aCTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATtctc              >  2:276073/1‑150 (MQ=255)
                                                                                                                                   cTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTctct             >  2:817691/1‑150 (MQ=255)
                                                                                                                                      cATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTc           <  1:179355/149‑1 (MQ=255)
                                                                                                                                              cTTGTGACCCATCAGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCgg  <  1:123075/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
CAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGG  >  CP009273/2959970‑2960260

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