Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,247,427 G→T 44.5% D129Y (GAC→TAC)  lexA → transcriptional repressor of SOS regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,247,4270GT44.5% 243.8 / inf 283D129Y (GAC→TAC) lexAtranscriptional repressor of SOS regulon
Reads supporting (aligned to +/- strand):  ref base G (76/81);  new base T (49/77);  total (125/158)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 284 total aligned reads displayed.

TAGGTCGTGTGGCTGCCGGTGAACCACTTCTGGCGCAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGT  >  CP009273/4247278‑4247575
                                                                                                                                                     |                                                                                                                                                    
 aGGTCGTGTGGCTGCCGGTGAACCACTTCTGGCGCAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGa                                                                                                                                                     >  2:981916/1‑150 (MQ=255)
   gTCGTGTGGCTGCCGGTGAACCACTTCTGGCGCAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACt                                                                                                                                                   <  1:4564566/150‑1 (MQ=255)
   gTCGTGTGGCTGCCGGTGAACCACTTCTGGCGCAACAGCATATAGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGTGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACt                                                                                                                                                   <  2:2981676/150‑1 (MQ=255)
                cGGTGAACCACTTCTGGCGCAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCa                                                                                                                                      <  2:2700436/150‑1 (MQ=255)
                cGGTGAACCACTTCTGGCGCAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCa                                                                                                                                      <  2:2700458/150‑1 (MQ=255)
                cGGTGAACCACTTCTGGCGCAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCa                                                                                                                                      <  2:1882035/150‑1 (MQ=255)
                   tGAACCACTTCTGGCGCAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGTATGGTGACTTGCTGGCAGTGCATaa                                                                                                                                   >  1:4209279/1‑150 (MQ=255)
                     aaCCACTTCTGGCGCAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTGATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATaaaa                                                                                                                                 <  2:3885734/150‑1 (MQ=255)
                      aCCACTTCTGGCGCAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAAc                                                                                                                                >  2:2309846/1‑150 (MQ=255)
                              tGGCGCAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTTATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATg                                                                                                                        <  2:1654748/150‑1 (MQ=255)
                              tGGCGCAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATg                                                                                                                        >  1:100024/1‑150 (MQ=255)
                                 cgcAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTAc                                                                                                                     >  2:4611482/1‑150 (MQ=255)
                                 cgcAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTAc                                                                                                                     <  1:4187671/150‑1 (MQ=255)
                                      cAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAAc                                                                                                                <  2:3782895/150‑1 (MQ=255)
                                       aGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACg                                                                                                               >  1:187762/1‑150 (MQ=255)
                                           tatTGAAGGTCATTATCAGGTCGATCCGTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCa                                                                                                           <  1:2695803/150‑1 (MQ=255)
                                              tGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGt                                                                                                        >  1:1700624/1‑150 (MQ=255)
                                                aaGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCg                                                                                                      <  1:1424798/150‑1 (MQ=255)
                                                 aGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGt                                                                                                     <  2:1214805/150‑1 (MQ=255)
                                                  ggTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGtt                                                                                                    <  1:537216/150‑1 (MQ=255)
                                                   gTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTg                                                                                                   >  1:1569958/1‑150 (MQ=255)
                                                   gTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTg                                                                                                   >  2:2976185/1‑150 (MQ=255)
                                                   gTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTg                                                                                                   <  1:4175115/150‑1 (MQ=255)
                                                      aTTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCAGTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCg                                                                                                >  1:1312350/1‑150 (MQ=255)
                                                        tATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCa                                                                                              >  2:3680803/1‑150 (MQ=255)
                                                        tATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCa                                                                                              >  2:1361639/1‑150 (MQ=255)
                                                         aTCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCAc                                                                                             <  1:558143/150‑1 (MQ=255)
                                                         aTCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCAc                                                                                             <  1:1144263/150‑1 (MQ=255)
                                                           cAGGTCGATCCTTNCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGt                                                                                           >  1:1409025/1‑150 (MQ=255)
                                                            aGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTa                                                                                          >  2:3648614/1‑150 (MQ=255)
                                                            aGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTa                                                                                          >  2:4112464/1‑150 (MQ=255)
                                                                cGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATtga                                                                                      >  2:49822/1‑150 (MQ=255)
                                                                 gATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATtgat                                                                                     <  1:1561425/150‑1 (MQ=255)
                                                                  aTCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATtgatg                                                                                    <  2:2526300/150‑1 (MQ=255)
                                                                   tccttcctTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATtgatga                                                                                   <  2:3915099/150‑1 (MQ=255)
                                                                   tccttcctTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATtgatga                                                                                   >  2:1447513/1‑150 (MQ=255)
                                                                    ccttccttATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGAc                                                                                  >  2:1900153/1‑150 (MQ=255)
                                                                        ccttATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAg                                                                              <  2:74480/150‑1 (MQ=255)
                                                                        ccttATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAg                                                                              <  2:3595282/150‑1 (MQ=255)
                                                                          ttATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGtt                                                                            <  2:1400574/150‑1 (MQ=255)
                                                                          ttATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGtt                                                                            <  2:2082492/150‑1 (MQ=255)
                                                                             ttCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTAcc                                                                         <  1:3721013/150‑1 (MQ=255)
                                                                             ttCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTAcc                                                                         <  1:3880372/150‑1 (MQ=255)
                                                                             ttCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTTGTTGTCGCACGTATTGATGACGAAGTTAcc                                                                         <  1:2378221/150‑1 (MQ=255)
                                                                                aaGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGtt                                                                      <  1:1829611/150‑1 (MQ=255)
                                                                                 aGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTa                                                                     <  1:148940/150‑1 (MQ=255)
                                                                                   ccGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAg                                                                   <  2:1330820/150‑1 (MQ=255)
                                                                                    cGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAgc                                                                  >  2:4120638/1‑150 (MQ=255)
                                                                                       aTGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGcc                                                               <  1:4144650/150‑1 (MQ=255)
                                                                                       aTGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGcc                                                               >  2:904977/1‑150 (MQ=255)
                                                                                         gCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTg                                                             <  1:2488906/150‑1 (MQ=255)
                                                                                          cTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGa                                                            >  1:3174457/1‑150 (MQ=255)
                                                                                          cTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGa                                                            <  2:608335/150‑1 (MQ=255)
                                                                                          cTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCACAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGa                                                            >  2:1744049/1‑150 (MQ=255)
                                                                                          cTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCACAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGa                                                            >  2:455514/1‑150 (MQ=255)
                                                                                           tGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGaa                                                           >  2:1005821/1‑150 (MQ=255)
                                                                                            gATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGaaa                                                          <  2:475852/150‑1 (MQ=255)
                                                                                            gATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGaaa                                                          <  2:3472646/150‑1 (MQ=255)
                                                                                               ttCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGaaaaaa                                                       >  2:2026119/1‑150 (MQ=255)
                                                                                                tCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAAc                                                      >  1:3761553/1‑150 (MQ=255)
                                                                                                  ctgctgCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAg                                                    <  1:2549118/150‑1 (MQ=255)
                                                                                                    gctgcGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAggg                                                  <  2:291808/150‑1 (MQ=255)
                                                                                                      tgcGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCa                                                <  1:4296005/150‑1 (MQ=255)
                                                                                                          cgTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATaa                                            <  2:3169319/150‑1 (MQ=255)
                                                                                                          cgTAAGCGGGATGTCGATGAAAGATATCGGAATGTTGGTTGGTTACTTGCTTTCAGTGCATAAAACTCTGGATGTACGTAACGGTCAGGTAGTTGTATCACTTATTGATGACGAAGTTACCGTTAAGCGCCTGATATATCAGGGCAATaa                                            <  2:734801/150‑1 (MQ=255)
                                                                                                           gTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATaaa                                           >  1:3941779/1‑150 (MQ=255)
                                                                                                              aGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTc                                        <  1:497984/150‑1 (MQ=255)
                                                                                                              aGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTc                                        <  1:301689/150‑1 (MQ=255)
                                                                                                              aGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTc                                        >  1:3296848/1‑150 (MQ=255)
                                                                                                               gCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCg                                       <  2:4011376/150‑1 (MQ=255)
                                                                                                                cGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGa                                      >  1:44520/1‑150 (MQ=255)
                                                                                                                cGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGa                                      >  2:2573280/1‑150 (MQ=255)
                                                                                                                 gtgATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGaa                                     <  2:3403203/148‑1 (MQ=255)
                                                                                                                 gggATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGaa                                     <  2:3402993/150‑1 (MQ=255)
                                                                                                                 gggATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGGATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGaa                                     >  2:3485856/1‑150 (MQ=255)
                                                                                                                  ggATGTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAAc                                    <  1:3387742/150‑1 (MQ=255)
                                                                                                                      gTCGATGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGtt                                <  2:830191/150‑1 (MQ=255)
                                                                                                                       tCGATGTTAGATATCGGCTTTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGATAATACAGGGCAATAATGTCGAACTGTTg                               <  2:183999/150‑1 (MQ=255)
                                                                                                                           tGAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAg                           <  2:2473781/150‑1 (MQ=255)
                                                                                                                           tGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAg                           >  2:2156374/1‑150 (MQ=255)
                                                                                                                            gAAAGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGa                          >  2:3246716/1‑150 (MQ=255)
                                                                                                                             aaaGATATCGGCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGaa                         >  2:2310341/1‑150 (MQ=255)
                                                                                                                                gATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAAt                      >  1:290040/1‑150 (MQ=255)
                                                                                                                                   atCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGc                   <  2:886100/150‑1 (MQ=255)
                                                                                                                                    tCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCg                  >  1:2028238/1‑150 (MQ=255)
                                                                                                                                     cGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGa                 >  2:2311844/1‑150 (MQ=255)
                                                                                                                                     cGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGa                 <  2:1107186/150‑1 (MQ=255)
                                                                                                                                      ggCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATagagag                >  2:831935/1‑150 (MQ=255)
                                                                                                                                       gCATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGt               >  2:2812300/1‑150 (MQ=255)
                                                                                                                                       gCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGt               <  2:4326722/150‑1 (MQ=255)
                                                                                                                                        cATTATGGATGGTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGtt              <  2:1669234/150‑1 (MQ=255)
                                                                                                                                        cATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGtt              >  1:3147790/1‑150 (MQ=255)
                                                                                                                                              ggatggTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAAcc        <  2:450668/150‑1 (MQ=255)
                                                                                                                                                atggTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCaa      >  1:3725369/1‑150 (MQ=255)
                                                                                                                                                atggTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCaa      <  1:2539990/150‑1 (MQ=255)
                                                                                                                                                atggTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCaa      <  1:2777323/150‑1 (MQ=255)
                                                                                                                                                 tggTTACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAAt     >  2:4373955/1‑150 (MQ=255)
                                                                                                                                                   gTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTg   <  2:4168703/150‑1 (MQ=255)
                                                                                                                                                    tGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGt  >  2:4132911/1‑150 (MQ=255)
                                                                                                                                                    tGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGt  <  1:733283/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TAGGTCGTGTGGCTGCCGGTGAACCACTTCTGGCGCAACAGCATATTGAAGGTCATTATCAGGTCGATCCTTCCTTATTCAAGCCGAATGCTGATTTCCTGCTGCGCGTCAGCGGGATGTCGATGAAAGATATCGGCATTATGGATGGTGACTTGCTGGCAGTGCATAAAACTCAGGATGTACGTAACGGTCAGGTCGTTGTCGCACGTATTGATGACGAAGTTACCGTTAAGCGCCTGAAAAAACAGGGCAATAAAGTCGAACTGTTGCCAGAAAATAGCGAGTTTAAACCAATTGT  >  CP009273/4247278‑4247575

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: