Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 899,961 T→A 8.4% intergenic (‑38/‑88) ybjP ← / → ybjQ lipoprotein/UPF0145 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273899,9610TA 8.4% 758.9 / 46.5 225intergenic (‑38/‑88)ybjP/ybjQlipoprotein/UPF0145 family protein
Reads supporting (aligned to +/- strand):  ref base T (85/121);  new base A (12/7);  total (97/128)
Fisher's exact test for biased strand distribution p-value = 8.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 232 total aligned reads displayed.

CAACGCAAGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGGTGTGGT  >  CP009273/899812‑900110
                                                                                                                                                     |                                                                                                                                                     
cAACGCAAGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCt                                                                                                                                                       >  1:200931/1‑150 (MQ=255)
cAACGCAAGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCt                                                                                                                                                       <  1:4330509/150‑1 (MQ=255)
 aaCGCAAGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTa                                                                                                                                                      <  2:1795818/150‑1 (MQ=255)
 aaCGCAAGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTa                                                                                                                                                      >  1:3515159/1‑150 (MQ=255)
      aaGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGata                                                                                                                                                 >  1:469596/1‑150 (MQ=255)
       aGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGataa                                                                                                                                                <  2:3841436/150‑1 (MQ=255)
        ggACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGataat                                                                                                                                               <  2:466130/150‑1 (MQ=255)
            cACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATg                                                                                                                                           >  1:2519661/1‑150 (MQ=255)
                gCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGatat                                                                                                                                       <  2:21983/150‑1 (MQ=255)
                 cGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATc                                                                                                                                      >  2:1484405/1‑150 (MQ=255)
                    gtgCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCtt                                                                                                                                   <  1:4136859/150‑1 (MQ=255)
                    gtgCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCtt                                                                                                                                   <  1:4135771/150‑1 (MQ=255)
                    gtgCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCtt                                                                                                                                   <  1:4378390/150‑1 (MQ=255)
                     tgacGTTATCTTTTTATGCTGGAGTGACTGTGGTGCAGGCACTGTGGAGCAGTGCGCAGGGGATAAGCATTGTCATTTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCttt                                                                                                                                  <  1:3589128/147‑1 (MQ=255)
                     tgCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCttt                                                                                                                                  <  2:781382/150‑1 (MQ=255)
                     tgCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCttt                                                                                                                                  <  2:3167241/150‑1 (MQ=255)
                         gTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCaat                                                                                                                              <  1:4427158/150‑1 (MQ=255)
                             tCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGTTAAGCATTGTCAATTTGCTGTAGCGCATTATGTTTTCCTTATAAGCGATCGATCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAataata                                                                                                                          <  1:1982971/150‑1 (MQ=255)
                             tCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAataata                                                                                                                          <  1:2160948/150‑1 (MQ=255)
                               tttATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGc                                                                                                                        <  2:3815624/150‑1 (MQ=255)
                               tttATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGc                                                                                                                        <  1:2341519/150‑1 (MQ=255)
                                 tataAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGTGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGt                                                                                                                      <  1:160626/150‑1 (MQ=255)
                                  ataAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTa                                                                                                                     <  2:3728899/150‑1 (MQ=255)
                                    aaGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATc                                                                                                                   >  2:51514/1‑150 (MQ=255)
                                     aGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCa                                                                                                                  <  2:2790341/150‑1 (MQ=255)
                                     aGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCa                                                                                                                  <  2:4005705/150‑1 (MQ=255)
                                          gAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGCATCAGTCTg                                                                                                             <  2:478270/150‑1 (MQ=255)
                                          gAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTg                                                                                                             <  2:4176716/150‑1 (MQ=255)
                                             tGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATa                                                                                                          >  2:1976298/1‑150 (MQ=255)
                                              gACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATaa                                                                                                         >  2:139952/1‑150 (MQ=255)
                                              gACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATaa                                                                                                         >  2:3413440/1‑150 (MQ=255)
                                               aCTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAAt                                                                                                        <  2:4402025/150‑1 (MQ=255)
                                               aCTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAAt                                                                                                        <  1:4176898/150‑1 (MQ=255)
                                               aCTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAAt                                                                                                        <  2:22923/150‑1 (MQ=255)
                                                 tgtgGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGc                                                                                                      >  2:1036941/1‑150 (MQ=255)
                                                    ggtgCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCttt                                                                                                   <  1:87084/150‑1 (MQ=255)
                                                    ggtgCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCttt                                                                                                   <  1:4510725/150‑1 (MQ=255)
                                                     gtgCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCtttt                                                                                                  <  2:750257/150‑1 (MQ=255)
                                                     gtgCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGAGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCtttt                                                                                                  <  1:698854/150‑1 (MQ=255)
                                                         aGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTgaga                                                                                              <  1:3330894/150‑1 (MQ=255)
                                                          ggCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGAt                                                                                             >  1:3561162/1‑150 (MQ=255)
                                                           gCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATc                                                                                            <  1:2219005/150‑1 (MQ=255)
                                                             aCTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATAGa                                                                                          >  2:2239305/1‑150 (MQ=255)
                                                             aCTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGa                                                                                          >  2:3171715/1‑150 (MQ=255)
                                                               tGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAg                                                                                        <  2:462503/150‑1 (MQ=255)
                                                               tGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAg                                                                                        >  2:3929014/1‑150 (MQ=255)
                                                                 aggagCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGc                                                                                      >  2:605215/1‑150 (MQ=255)
                                                                   gagCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCtt                                                                                    <  1:800520/150‑1 (MQ=255)
                                                                   gagCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCtt                                                                                    <  1:2297578/150‑1 (MQ=255)
                                                                     gCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAg                                                                                  >  2:3508673/1‑150 (MQ=255)
                                                                       aGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCa                                                                                >  1:681268/1‑150 (MQ=255)
                                                                        gTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCaa                                                                               >  1:2923097/1‑150 (MQ=255)
                                                                                  gATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgat                                                                     <  2:2095559/150‑1 (MQ=255)
                                                                                   aTAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgatc                                                                    <  1:1115517/150‑1 (MQ=255)
                                                                                   aTAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgatc                                                                    >  2:1254608/1‑150 (MQ=255)
                                                                                   aTAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgatc                                                                    >  2:1774681/1‑150 (MQ=255)
                                                                                    tAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgatcg                                                                   <  2:3592521/150‑1 (MQ=255)
                                                                                      aGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgatcgat                                                                 <  1:30630/150‑1 (MQ=255)
                                                                                      aGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAgatcgat                                                                 <  1:31000/150‑1 (MQ=255)
                                                                                        cATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCa                                                               <  1:756608/150‑1 (MQ=255)
                                                                                        cATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCa                                                               <  2:4755309/150‑1 (MQ=255)
                                                                                        cATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCa                                                               >  2:2617922/1‑150 (MQ=255)
                                                                                        cATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCa                                                               <  1:1482513/150‑1 (MQ=255)
                                                                                        cATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCa                                                               >  1:640943/1‑150 (MQ=255)
                                                                                            gTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCa                                                           <  1:3909849/150‑1 (MQ=255)
                                                                                            gTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCa                                                           <  2:1665649/150‑1 (MQ=255)
                                                                                                aTTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAAttt                                                       <  1:172396/150‑1 (MQ=255)
                                                                                                   tGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCt                                                    <  2:2472874/150‑1 (MQ=255)
                                                                                                   tGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCt                                                    <  2:469596/150‑1 (MQ=255)
                                                                                                      tGTAGCTCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACa                                                 <  1:2457037/150‑1 (MQ=255)
                                                                                                      tGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACa                                                 <  1:2722012/150‑1 (MQ=255)
                                                                                                      tGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACa                                                 <  2:1809309/150‑1 (MQ=255)
                                                                                                        tAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAAc                                               >  2:3025454/1‑150 (MQ=255)
                                                                                                           cgcATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAAcccc                                            >  2:3693713/1‑150 (MQ=255)
                                                                                                             cATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCaa                                          >  2:1959683/1‑150 (MQ=255)
                                                                                                              aTAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAAc                                         <  1:2201150/150‑1 (MQ=255)
                                                                                                              aTAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAAc                                         >  2:3790935/1‑150 (MQ=255)
                                                                                                              aTAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAAc                                         >  2:3793181/1‑150 (MQ=255)
                                                                                                                 aTGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAActct                                      <  1:2088520/150‑1 (MQ=255)
                                                                                                                   gATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTgg                                    <  2:569305/150‑1 (MQ=255)
                                                                                                                      ttCCTTATAAGCGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAg                                 >  2:2431578/1‑150 (MQ=255)
                                                                                                                         cTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGcc                              >  2:672207/1‑150 (MQ=255)
                                                                                                                          ttATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCa                             >  2:737115/1‑150 (MQ=255)
                                                                                                                           tataAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAg                            >  1:3741620/1‑150 (MQ=255)
                                                                                                                               aGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCa                        >  2:2983554/1‑150 (MQ=255)
                                                                                                                                 cGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATc                      >  2:1752942/1‑150 (MQ=255)
                                                                                                                                 cGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATc                      >  2:4500523/1‑150 (MQ=255)
                                                                                                                                 cGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATc                      >  1:1247432/1‑150 (MQ=255)
                                                                                                                                 cGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATc                      >  1:4020151/1‑150 (MQ=255)
                                                                                                                                 cGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATc                      >  1:3777816/1‑150 (MQ=255)
                                                                                                                                 cGATCGCTCTGAAAGCGTTCAACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATc                      >  1:3034629/1‑150 (MQ=255)
                                                                                                                                   aTCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGt                    >  2:4517575/1‑150 (MQ=255)
                                                                                                                                    tCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGtt                   >  1:2747101/1‑150 (MQ=255)
                                                                                                                                       ctctGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGaa                <  1:4011091/150‑1 (MQ=255)
                                                                                                                                         ctGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAata              <  1:1844150/150‑1 (MQ=255)
                                                                                                                                         ctGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAata              >  1:3697050/1‑150 (MQ=255)
                                                                                                                                          tGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAatat             >  1:4626000/1‑150 (MQ=255)
                                                                                                                                                cGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGgt       >  2:3296996/1‑150 (MQ=255)
                                                                                                                                                 gTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGgtg      <  2:4065298/150‑1 (MQ=255)
                                                                                                                                                 gTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGgtg      <  2:200931/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAACGCAAGGACCACTGCGTGTGCCGTTATCTTTATAAGCTGGAGTGACTGTGGTGCAGGCACTGAGGAGCAGTGCGCAGGGGATAAGCATTGTCAATTTGCTGTAGCGCATAATGATTTCCTTATAAGCGATCGCTCTGAAAGCGTTCTACGATAATAATGATATCCTTTCAATAATAGCGTATCAGTCTGATAATGCTTTTGAGATCGAAGGCTTAGCAAACAAGGAGATCGATCATGCAATTTTCTACAACCCCAACTCTGGAAGGCCAGACCATCGTTGAATATTGCGGTGTGGT  >  CP009273/899812‑900110

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: