Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,246,967 G→T 6.3% intergenic (+34/‑76) dgkA → / → lexA diacylglycerol kinase/transcriptional repressor of SOS regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,246,9670GT 6.3% 1172.5 / 39.2 302intergenic (+34/‑76)dgkA/lexAdiacylglycerol kinase/transcriptional repressor of SOS regulon
Reads supporting (aligned to +/- strand):  ref base G (151/132);  new base T (10/9);  total (161/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 309 total aligned reads displayed.

TGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGC  >  CP009273/4246819‑4247115
                                                                                                                                                    |                                                                                                                                                    
tGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATtg                                                                                                                                                      >  1:197060/1‑149 (MQ=255)
tGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATtg                                                                                                                                                      <  2:2802006/149‑1 (MQ=255)
  aaTACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATtgtg                                                                                                                                                    >  1:4499079/1‑149 (MQ=255)
    tACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCa                                                                                                                                                  >  2:3078698/1‑149 (MQ=255)
      ccATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGt                                                                                                                                                <  2:3825406/149‑1 (MQ=255)
          gAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTAt                                                                                                                                            <  1:1526689/149‑1 (MQ=255)
            gCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATgg                                                                                                                                          >  1:4717894/1‑149 (MQ=255)
              tttCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGtt                                                                                                                                        <  1:4132693/149‑1 (MQ=255)
               ttCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTc                                                                                                                                       >  2:5188262/1‑149 (MQ=255)
                tCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTcc                                                                                                                                      <  2:1455693/149‑1 (MQ=255)
                     aCGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAAt                                                                                                                                 <  1:1878235/149‑1 (MQ=255)
                         gcAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGcc                                                                                                                             <  2:4947370/149‑1 (MQ=255)
                         gcAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGcc                                                                                                                             <  1:3078698/149‑1 (MQ=255)
                           aaaaGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCtt                                                                                                                           >  1:2420431/1‑149 (MQ=255)
                            aaaGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCttt                                                                                                                          <  1:3672700/149‑1 (MQ=255)
                            aaaGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCttt                                                                                                                          >  2:2483655/1‑149 (MQ=255)
                               gATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGc                                                                                                                       >  2:2247894/1‑149 (MQ=255)
                               gATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGc                                                                                                                       >  2:2263576/1‑149 (MQ=255)
                               gATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGc                                                                                                                       >  2:4516109/1‑149 (MQ=255)
                                atatGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTTTGCAGTTTATGGTTCCAAAATCGCCTTTTGCt                                                                                                                      <  2:1990022/149‑1 (MQ=255)
                                 tatGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTTTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTg                                                                                                                     <  2:1003415/149‑1 (MQ=255)
                                 tatGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTg                                                                                                                     <  2:4999004/149‑1 (MQ=255)
                                 tatGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTg                                                                                                                     <  1:2661567/149‑1 (MQ=255)
                                 tatGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTg                                                                                                                     <  2:3480542/149‑1 (MQ=255)
                                  atGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGt                                                                                                                    >  2:82989/1‑149 (MQ=255)
                                   tGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGta                                                                                                                   >  1:2853389/1‑149 (MQ=255)
                                    gggATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGtat                                                                                                                  <  2:3495124/149‑1 (MQ=255)
                                    gggATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGtat                                                                                                                  >  2:131005/1‑149 (MQ=255)
                                    gggATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGtat                                                                                                                  <  1:2434200/149‑1 (MQ=255)
                                      gATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGtatat                                                                                                                >  2:1855882/1‑149 (MQ=255)
                                      gATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGtatat                                                                                                                <  2:3965989/149‑1 (MQ=255)
                                       aTCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGtatata                                                                                                               <  2:2472466/149‑1 (MQ=255)
                                             tGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTcaca                                                                                                         <  2:3210019/149‑1 (MQ=255)
                                             tGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTcaca                                                                                                         <  1:4605228/149‑1 (MQ=255)
                                                ggTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCa                                                                                                      >  1:4041016/1‑149 (MQ=255)
                                                 gTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTTTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCTCAGCAt                                                                                                     >  2:2383511/1‑149 (MQ=255)
                                                 gTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCAt                                                                                                     >  2:589570/1‑149 (MQ=255)
                                                 gTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCAt                                                                                                     >  1:1489435/1‑149 (MQ=255)
                                                   gCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATaa                                                                                                   <  2:1898092/149‑1 (MQ=255)
                                                        ttGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGta                                                                                              >  1:2987893/1‑149 (MQ=255)
                                                         tGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGtat                                                                                             >  1:1283500/1‑149 (MQ=255)
                                                         tGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGtat                                                                                             >  1:3878562/1‑149 (MQ=255)
                                                          gCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTTTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGtata                                                                                            <  2:2877666/149‑1 (MQ=255)
                                                          gCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGtata                                                                                            <  2:173293/149‑1 (MQ=255)
                                                          gCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGtata                                                                                            <  1:5115063/149‑1 (MQ=255)
                                                            cATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGtatata                                                                                          >  1:6928/1‑149 (MQ=255)
                                                            cATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGtatata                                                                                          >  2:5149992/1‑149 (MQ=255)
                                                             aTTATCGTCGCCGTGATTACCTGGTGCAGTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATac                                                                                         <  2:4298549/149‑1 (MQ=255)
                                                              ttATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATaca                                                                                        >  2:1978538/1‑149 (MQ=255)
                                                                 tcgtcgCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACAccc                                                                                     >  1:2758017/1‑149 (MQ=255)
                                                                  cgtcgCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCa                                                                                    <  1:2461937/149‑1 (MQ=255)
                                                                   gtcgCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAg                                                                                   >  1:439585/1‑149 (MQ=255)
                                                                        cGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGc                                                                              <  2:3882419/149‑1 (MQ=255)
                                                                          tGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGTGGGCgg                                                                            >  2:3283171/1‑149 (MQ=255)
                                                                          tGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAAGTCGATAAATCTCTGGTTTCGTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCgg                                                                            >  1:2068188/1‑149 (MQ=255)
                                                                             ttACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTTTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAAt                                                                         <  2:4138176/149‑1 (MQ=255)
                                                                                 cTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATTaaa                                                                     >  2:1689760/1‑149 (MQ=255)
                                                                                   ggTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGc                                                                   <  2:4831500/149‑1 (MQ=255)
                                                                                    gTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCg                                                                  >  2:4564629/1‑149 (MQ=255)
                                                                                        aTTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTaa                                                              >  1:3795724/1‑149 (MQ=255)
                                                                                          tCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACg                                                            <  1:3618791/149‑1 (MQ=255)
                                                                                          tCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACg                                                            <  1:2160748/149‑1 (MQ=255)
                                                                                               tATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTTTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAg                                                       >  1:2692546/1‑149 (MQ=255)
                                                                                               tATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAg                                                       <  2:2678595/149‑1 (MQ=255)
                                                                                                   gTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGcac                                                   >  2:1337770/1‑148 (MQ=255)
                                                                                                   gTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGcaa                                                   >  1:3527866/1‑149 (MQ=255)
                                                                                                      gCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGcaacaa                                                <  2:483366/149‑1 (MQ=255)
                                                                                                      gCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGcaacaa                                                <  2:3768140/149‑1 (MQ=255)
                                                                                                        aTTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAaga                                              >  1:1094390/1‑149 (MQ=255)
                                                                                                         ttttGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAagag                                             <  1:1341672/149‑1 (MQ=255)
                                                                                                          tttGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAgg                                            <  2:3397824/149‑1 (MQ=255)
                                                                                                           ttGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGgt                                           <  2:5085780/149‑1 (MQ=255)
                                                                                                             ggATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGgtgt                                         <  1:926594/149‑1 (MQ=255)
                                                                                                               aTAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGttt                                       >  1:5006229/1‑149 (MQ=255)
                                                                                                                tAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTg                                      >  1:252391/1‑149 (MQ=255)
                                                                                                                   cccTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGAtc                                   >  1:3134836/1‑149 (MQ=255)
                                                                                                                   cccTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGAtc                                   >  1:3950655/1‑149 (MQ=255)
                                                                                                                     cTTCCAGAATTCGATAAATCTCTGGTTTATTTTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGAtctc                                 <  1:774392/149‑1 (MQ=255)
                                                                                                                       tCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCAt                               >  1:4887118/1‑149 (MQ=255)
                                                                                                                       tCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCAt                               <  2:753331/149‑1 (MQ=255)
                                                                                                                       tCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCAt                               <  1:4282816/149‑1 (MQ=255)
                                                                                                                       tCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCAt                               >  2:4287609/1‑149 (MQ=255)
                                                                                                                       tCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCAt                               <  2:2065662/149‑1 (MQ=255)
                                                                                                                          aGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCg                            >  1:2329732/1‑149 (MQ=255)
                                                                                                                               tCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATc                       <  1:317875/149‑1 (MQ=255)
                                                                                                                                 gATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATcac                     >  1:4355426/1‑149 (MQ=255)
                                                                                                                                 gATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATcac                     >  1:992066/1‑149 (MQ=255)
                                                                                                                                  aTAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATcaca                    >  2:4408102/1‑149 (MQ=255)
                                                                                                                                   tAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCCGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACAt                   >  1:3425880/1‑149 (MQ=255)
                                                                                                                                     aaTCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCa                 >  1:3265903/1‑149 (MQ=255)
                                                                                                                                            ggTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGAc          <  2:781524/149‑1 (MQ=255)
                                                                                                                                             gTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACa         <  2:473238/149‑1 (MQ=255)
                                                                                                                                             gTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACa         <  2:4623247/149‑1 (MQ=255)
                                                                                                                                               ttATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAgg       >  2:4434213/1‑149 (MQ=255)
                                                                                                                                                tATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGt      <  2:3658586/149‑1 (MQ=255)
                                                                                                                                                 aTTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTa     <  1:1235367/149‑1 (MQ=255)
                                                                                                                                                  ttGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTAt    >  1:4772447/1‑149 (MQ=255)
                                                                                                                                                  ttGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTAt    >  2:3548961/1‑149 (MQ=255)
                                                                                                                                                  ttGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTAt    >  2:783366/1‑149 (MQ=255)
                                                                                                                                                    gtgCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATgc  <  1:1004443/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGAATACCATGAGCTTTCCGGACGCGCAAAAGATATGGGATCCGCTGCGGTGCTGATTGCCATTATCGTCGCCGTGATTACCTGGTGCATTCTGTTATGGTCGCATTTTGGATAACCCTTCCAGAATTCGATAAATCTCTGGTTTATTGTGCAGTTTATGGTTCCAAAATCGCCTTTTGCTGTATATACTCACAGCATAACTGTATATACACCCAGGGGGCGGAATGAAAGCGTTAACGGCCAGGCAACAAGAGGTGTTTGATCTCATCCGTGATCACATCAGCCAGACAGGTATGC  >  CP009273/4246819‑4247115

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: