Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 426,157 G→T 100% E33* (GAG→TAG)  yajD → HNH nuclease family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273426,1570GT96.5% 902.2 / 11.4 352E33* (GAG→TAG) yajDHNH nuclease family protein
Reads supporting (aligned to +/- strand):  ref base G (14/4);  new base T (157/176);  total (172/180)
Fisher's exact test for biased strand distribution p-value = 1.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 352 total aligned reads displayed.

CGCATCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCA  >  CP009273/426010‑426302
                                                                                                                                                   |                                                                                                                                                 
cGCATCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGTGGTCGCTGTTCCCGCTAg                                                                                                                                                 <  1:2580055/150‑1 (MQ=255)
  cATCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTACCGCGAgtt                                                                                                                                               >  1:1817682/1‑150 (MQ=255)
         ggATACTGCGAGTAACCCTACATTTCATCCATGTAACTTTTCCTGGCTATCCTCCCAAAAATCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTAt                                                                                                                                        <  1:4710527/150‑1 (MQ=255)
          gAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAACATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTAtt                                                                                                                                       <  1:902868/150‑1 (MQ=255)
           aaaCTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTc                                                                                                                                      <  2:2942678/150‑1 (MQ=255)
             aCTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCa                                                                                                                                    <  2:583612/150‑1 (MQ=255)
             aCTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCaa                                                                                                                                   <  2:3635451/150‑1 (MQ=255)
                 cgcgTAACCCTACATTTCCTACAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCt                                                                                                                                <  1:914213/150‑1 (MQ=255)
                 cgcgTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCt                                                                                                                                <  2:4273767/150‑1 (MQ=255)
                      aaCCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGTTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTg                                                                                                                           >  2:4895818/1‑150 (MQ=255)
                      aaCCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTg                                                                                                                           <  2:4806606/150‑1 (MQ=255)
                        cccTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGaa                                                                                                                         <  1:2997517/150‑1 (MQ=255)
                         ccTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGacc                                                                                                                        >  2:3042919/1‑148 (MQ=255)
                         ccTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAAc                                                                                                                        <  1:4899156/150‑1 (MQ=255)
                           tACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTACAACCTGCGTGAACtt                                                                                                                      >  1:1220748/1‑150 (MQ=255)
                             cATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTaa                                                                                                                    >  2:4122583/1‑149 (MQ=255)
                              aTTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTAcc                                                                                                                   >  2:3832530/1‑150 (MQ=255)
                              aTTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTAcc                                                                                                                   >  2:5147982/1‑150 (MQ=255)
                               tttCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCg                                                                                                                  >  1:3904601/1‑150 (MQ=255)
                                 tCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGtt                                                                                                                >  1:5217310/1‑150 (MQ=255)
                                  cATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTc                                                                                                               <  1:292957/150‑1 (MQ=255)
                                  cATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTc                                                                                                               >  2:243353/1‑150 (MQ=255)
                                   aTCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAATAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTACAACCTGCGTTAACTTACCGTTca                                                                                                              >  2:2685360/1‑150 (MQ=255)
                                   aTCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTCGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTca                                                                                                              >  2:1863910/1‑150 (MQ=255)
                                   aTCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTca                                                                                                              <  2:3152576/150‑1 (MQ=255)
                                      cAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTcacca                                                                                                           <  2:3227971/150‑1 (MQ=255)
                                       aGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGGTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTcaccac                                                                                                          <  1:4749526/150‑1 (MQ=255)
                                         gTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACcacaa                                                                                                        >  1:2785293/1‑149 (MQ=255)
                                           aaCTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTCTCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTg                                                                                                      >  1:4547873/1‑150 (MQ=255)
                                           aaCTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTg                                                                                                      >  2:1926781/1‑150 (MQ=255)
                                            aCTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGGTTGTGTATTCCAACCTGCGTGTCCTTACCGGTCACCCCATTGa                                                                                                     >  2:5086663/1‑150 (MQ=255)
                                               tttCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCa                                                                                                  >  2:3266905/1‑150 (MQ=255)
                                               tttCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCa                                                                                                  <  1:3276109/150‑1 (MQ=255)
                                                  cATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACAGTGATCACGa                                                                                               >  2:4640950/1‑150 (MQ=255)
                                                       cTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATa                                                                                          >  2:5194397/1‑150 (MQ=255)
                                                         atcatcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATAcc                                                                                        <  2:4475722/150‑1 (MQ=255)
                                                            atcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAAt                                                                                     <  2:4394051/150‑1 (MQ=255)
                                                            atcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAAt                                                                                     <  2:200500/150‑1 (MQ=255)
                                                            atcCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAAt                                                                                     >  2:5016282/1‑150 (MQ=255)
                                                               ccAAAAAACTAGGCGCGGGTAGAAAGCGGCTATCGCGAAACAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATacc                                                                                  <  1:169042/150‑3 (MQ=255)
                                                                cAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAAcc                                                                                 <  2:100680/150‑1 (MQ=255)
                                                                 aaaaaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTACAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAAAAAccc                                                                                >  2:1834181/1‑150 (MQ=255)
                                                                 aaaaaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGGTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAAccc                                                                                >  1:788922/1‑150 (MQ=255)
                                                                    aaaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTTTTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAAAAAACCGGa                                                                             >  1:2702827/1‑150 (MQ=255)
                                                                     aaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGaa                                                                            <  2:222969/150‑1 (MQ=255)
                                                                     aaCTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTACCGCTAGTTTGTTTATTACAACCTGCGTGAACATACCCTTCACCACAATGATCACGACCATACCAAAAACCCCGaa                                                                            >  1:1601234/1‑150 (MQ=255)
                                                                     aaCTATGCGCGGTTAGAAAGCGGCTATCGAGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGaa                                                                            >  2:1732047/1‑150 (MQ=255)
                                                                        tATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGAt                                                                         <  1:3629635/150‑1 (MQ=255)
                                                                        tATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGAt                                                                         <  1:1926236/150‑1 (MQ=255)
                                                                         aTGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATg                                                                        >  2:1496957/1‑150 (MQ=255)
                                                                          tGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATgg                                                                       <  1:2098218/150‑1 (MQ=255)
                                                                            cgcgGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCa                                                                     <  1:1037906/150‑1 (MQ=255)
                                                                              cgGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTACCGCTAGTTTGTTTATTAAAAACTGCGTGAACTTACCGTTCAACACATTGATCAAGACCATACCAAAAACCCGGAAGATGGCAGt                                                                   >  1:65628/1‑150 (MQ=255)
                                                                                gTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTaa                                                                 >  2:5228584/1‑150 (MQ=255)
                                                                                  tAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACt                                                               <  1:2460836/150‑1 (MQ=255)
                                                                                   aGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCCGTAACTg                                                              >  2:3995703/1‑150 (MQ=255)
                                                                                    gAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTgg                                                             >  2:3596461/1‑150 (MQ=255)
                                                                                    gAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACAATGATCACGACCATACCCATAACCCGGAAGATGGCAGTAACTgg                                                             >  1:1963258/1‑150 (MQ=255)
                                                                                       aGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGaa                                                          <  1:5063379/150‑1 (MQ=255)
                                                                                            cTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGACttgtt                                                     >  1:2831394/1‑150 (MQ=255)
                                                                                            cTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAAttgtt                                                     <  2:3900371/150‑1 (MQ=255)
                                                                                             tATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAAttgttg                                                    >  1:1459330/1‑150 (MQ=255)
                                                                                              aTCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATtgttgt                                                   <  2:1345895/150‑1 (MQ=255)
                                                                                              aTCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATtgttgt                                                   <  1:3710428/150‑1 (MQ=255)
                                                                                              aTCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGATGTTCCCGCTGGTTTGTTTATTCAAAAATGCGTGAACATACCGTTAACCACATTGATCACGACCAAACCCAAAACCCCGAAGATGGCAGTAACAGGGTATtgttgt                                                   >  1:1169855/1‑150 (MQ=255)
                                                                                                cgcgAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTCTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTgtgt                                                 <  2:2156173/150‑1 (MQ=255)
                                                                                                     aaaaGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTa                                            >  2:143629/1‑150 (MQ=255)
                                                                                                     aaaaGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCAGCGTGAACTTTCCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGGTACTGGGAATTGTTGTGTCTCTa                                            >  1:4259037/1‑150 (MQ=255)
                                                                                                       aatcATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACAACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTAtt                                          <  2:3689264/147‑1 (MQ=255)
                                                                                                        aGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTg                                         <  1:4796755/150‑1 (MQ=255)
                                                                                                         gCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTAACCACATTTATCACGACCAAACCAAAAACCCGGAAGATGGCAGTAACAGGGGATTGTTGTGTCTCTATTGc                                        >  1:4152790/1‑150 (MQ=255)
                                                                                                           aTTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCa                                      >  2:652380/1‑150 (MQ=255)
                                                                                                           aTTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTACAACCAGCGTGAACTTACCGTTCACCACAATTATAACGACCAAACCAATAAACCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTTCCa                                      >  1:1558440/1‑150 (MQ=255)
                                                                                                             tAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACg                                    <  2:1132258/150‑1 (MQ=255)
                                                                                                              aaaaaTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGa                                   >  1:4312489/1‑150 (MQ=255)
                                                                                                                aaaTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATc                                 <  1:724220/150‑1 (MQ=255)
                                                                                                                 aaTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCa                                >  2:730488/1‑150 (MQ=255)
                                                                                                                 aaTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCGATTGCCACGATCa                                >  2:4567992/1‑150 (MQ=255)
                                                                                                                  aTCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCAt                               <  2:2602040/150‑1 (MQ=255)
                                                                                                                    cTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGa                             >  2:1246483/1‑150 (MQ=255)
                                                                                                                    cTATCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGa                             <  1:4639558/150‑1 (MQ=255)
                                                                                                                      aTCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGc                           >  2:1769714/1‑150 (MQ=255)
                                                                                                                       tCCGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCa                          <  2:4947939/150‑1 (MQ=255)
                                                                                                                         cGTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCTGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCAtt                        <  1:366390/150‑1 (MQ=255)
                                                                                                                          gTGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTc                       <  1:2182674/150‑1 (MQ=255)
                                                                                                                           tGGGTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCg                      <  1:1742263/150‑1 (MQ=255)
                                                                                                                            gggTCTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGa                     <  1:4426793/150‑1 (MQ=255)
                                                                                                                                cTGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAata                 <  1:243353/150‑1 (MQ=255)
                                                                                                                                 tGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAatat                >  2:345563/1‑150 (MQ=255)
                                                                                                                                 tGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAatat                <  2:4548059/150‑1 (MQ=255)
                                                                                                                                 tGCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAatat                <  1:3976394/150‑1 (MQ=255)
                                                                                                                                  gCGGTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTTAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGTCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAAtata               <  2:4836173/150‑1 (MQ=255)
                                                                                                                                    ggTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATAcc             >  2:3771666/1‑150 (MQ=255)
                                                                                                                                    ggTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATAcc             >  1:1281828/1‑150 (MQ=255)
                                                                                                                                     gTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTCACTGGGCATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGACATATACCg            >  2:227110/1‑150 (MQ=255)
                                                                                                                                     gTCGCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCg            <  2:4152625/150‑1 (MQ=255)
                                                                                                                                        gCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAg         <  1:559197/150‑1 (MQ=255)
                                                                                                                                        gCTGTTCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAg         <  1:1113250/150‑1 (MQ=255)
                                                                                                                                            ttCCCGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGa     >  1:1916914/1‑150 (MQ=255)
                                                                                                                                               ccGCTAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCa  <  2:2782025/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGCATCTTGGGAAACTGCGCGTAACCCTACATTTCATCCAGGTAACTTTTCATGGCTATCATCCCAAAAAACTATGCGCGGTTAGAAAGCGGCTATCGCGAAAAAGCATTAAAAATCTATCCGTGGGTCTGCGGTCGCTGTTCCCGCGAGTTTGTTTATTCCAACCTGCGTGAACTTACCGTTCACCACATTGATCACGACCATACCAATAACCCGGAAGATGGCAGTAACTGGGAATTGTTGTGTCTCTATTGCCACGATCATGAGCATTCGAAATATACCGAAGCGGATCA  >  CP009273/426010‑426302

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: