Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,170,708 G→T 6.3% intergenic (‑25/+283) gatY ← / ← fbaB D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,170,7080GT 6.3% 697.2 / 11.8 242intergenic (‑25/+283)gatY/fbaBD‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base G (102/117);  new base T (4/19);  total (106/136)
Fisher's exact test for biased strand distribution p-value = 7.60e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 244 total aligned reads displayed.

CAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGG  >  CP009273/2170561‑2170856
                                                                                                                                                   |                                                                                                                                                     
cAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGa                                                                                                                                                      <  2:648525/149‑1 (MQ=255)
cAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGa                                                                                                                                                      <  2:2750325/149‑1 (MQ=255)
 aGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAt                                                                                                                                                     <  2:4179174/149‑1 (MQ=255)
   gTTGGCAGTGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAttt                                                                                                                                                   <  2:1708583/149‑1 (MQ=255)
   gTTGGCAGTGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAttt                                                                                                                                                   <  2:1708901/149‑1 (MQ=255)
   gTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAttt                                                                                                                                                   <  2:43514/149‑1 (MQ=255)
   gTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGAttt                                                                                                                                                   <  2:44123/149‑1 (MQ=255)
   gTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGtttt                                                                                                                                                   >  1:1390302/1‑149 (MQ=255)
    ttGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTc                                                                                                                                                  <  2:2506120/149‑1 (MQ=255)
    ttGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTc                                                                                                                                                  <  2:1878644/149‑1 (MQ=255)
         aGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAat                                                                                                                                             <  2:3335689/149‑1 (MQ=255)
          gCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAata                                                                                                                                            <  1:2802928/149‑1 (MQ=255)
           cGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAatat                                                                                                                                           <  2:1498096/149‑1 (MQ=255)
            ggTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAtata                                                                                                                                          <  2:2331418/149‑1 (MQ=255)
            ggTGTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAtata                                                                                                                                          <  2:2848397/149‑1 (MQ=255)
             gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATaa                                                                                                                                         <  2:4509139/149‑1 (MQ=255)
             gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATaa                                                                                                                                         <  2:584870/149‑1 (MQ=255)
             gTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTAGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGAATTCTAAaaatta                                                                                                                                         >  1:2491468/1‑147 (MQ=255)
              tttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAAt                                                                                                                                        <  2:2766079/149‑1 (MQ=255)
               ttCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATg                                                                                                                                       <  2:552169/149‑1 (MQ=255)
                  tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAAAATAAAGaaa                                                                                                                                    >  2:3055318/1‑149 (MQ=255)
                  tACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTATATTTAAAAAAATAATTaaa                                                                                                                                    >  1:2031671/1‑149 (MQ=255)
                   accaccACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAat                                                                                                                                   <  2:2465470/149‑1 (MQ=255)
                     caccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAAtttaaata                                                                                                                                 >  2:1952889/1‑142 (MQ=255)
                       ccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattatt                                                                                                                              <  1:2838293/149‑1 (MQ=255)
                       ccacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAattatt                                                                                                                              <  1:781453/149‑1 (MQ=255)
                        cacTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAAAAAAAAGAAAATAttt                                                                                                                             >  1:2580365/1‑149 (MQ=255)
                         acTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTg                                                                                                                            <  1:1491765/149‑1 (MQ=255)
                         acTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTg                                                                                                                            <  1:3881292/149‑1 (MQ=255)
                         acTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAAAATAATTAAAATAtttt                                                                                                                            >  2:1297921/1‑148 (MQ=255)
                           ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATTAAAATATTTGtt                                                                                                                          >  2:4007398/1‑149 (MQ=255)
                           ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATTAAATTAtttttt                                                                                                                          >  2:1834100/1‑149 (MQ=255)
                           ttGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTGTGTCGATCCCACGTACAGTTTCATCTCCTGTCGTGTGTGTTCGATTTGCAAATATAAGTGAATTAgttgtt                                                                                                                          >  2:4007450/1‑149 (MQ=255)
                             gCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAAAATAATGAAATTAtttttttt                                                                                                                        >  2:1046296/1‑149 (MQ=255)
                              cATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGttttt                                                                                                                       <  2:2044581/149‑1 (MQ=255)
                              cATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATAAAAAGAAATTATTTGttttt                                                                                                                       >  2:219859/1‑149 (MQ=255)
                                tCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTTTTTTAGATTTCAAAAAATAATTAAAATATTTTTTTTTaa                                                                                                                     >  1:1468457/1‑149 (MQ=255)
                                 cGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTaaa                                                                                                                    <  2:1553692/149‑1 (MQ=255)
                                    ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                 <  2:3976489/149‑1 (MQ=255)
                                    ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAatat                                                                                                                 <  2:4677824/149‑1 (MQ=255)
                                    ttCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTAAAAAAATAAATAAAATATTTTTTTTTAAatat                                                                                                                 >  2:4133829/1‑149 (MQ=255)
                                       gagaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCga                                                                                                              <  1:4600780/149‑1 (MQ=255)
                                        agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAATGTAATTATTTTTTTTTTAatatattg                                                                                                             >  2:695029/1‑145 (MQ=255)
                                        agaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTAGATTTAAAAAAATAAATAAAATAATTTTTTTTAAATATagag                                                                                                             >  2:1607336/1‑149 (MQ=255)
                                         gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCgaga                                                                                                            >  2:3765914/1‑149 (MQ=255)
                                         gaTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTCGATTTCAAAATATAATTAAAATATTTTTTTTTAAATATAgaga                                                                                                            >  2:252321/1‑149 (MQ=255)
                                         gaTTGTGAAGATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCgaga                                                                                                            <  2:1173017/149‑1 (MQ=255)
                                          aTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGAt                                                                                                           <  2:4112136/149‑1 (MQ=255)
                                          aTTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGAt                                                                                                           <  2:502323/149‑1 (MQ=255)
                                           ttGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATa                                                                                                          <  2:3894144/149‑1 (MQ=255)
                                           ttGTCAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACAGTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATAAAATGAAATTATTTGTTTTTAAAAATCGAtata                                                                                                          >  2:1796162/1‑149 (MQ=255)
                                             gtgAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAAc                                                                                                        <  2:3224485/149‑1 (MQ=255)
                                                 atatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCGTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTGGCTTTCCAAATATAAATTAATTATTTGTTTTTTTATATAGCGATAAccaaa                                                                                                    >  2:3977872/1‑145 (MQ=255)
                                                  tatTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAAAGAAATTATTTGTTTTTAAATATCGAGATAACGATca                                                                                                   >  2:1486363/1‑149 (MQ=255)
                                                    ttGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATAGAGATAAAGATcaca                                                                                                 >  1:1174937/1‑149 (MQ=255)
                                                    ttGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGGTTTGGATTTAAAAAAATAAAGTAATTAATTTTTTTTTAAAAAACGGATAACGaagaca                                                                                                 >  1:155343/1‑144 (MQ=255)
                                                    ttGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTGGATTTAAAAATATAATGAAATTATTTTTTTTTAAATATATAGATAAAGATaaca                                                                                                 >  1:3720565/1‑149 (MQ=255)
                                                        aTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaaa                                                                                             <  2:946641/149‑1 (MQ=255)
                                                        aTGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACaaaaa                                                                                             <  1:836885/149‑1 (MQ=255)
                                                         tGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAAc                                                                                            <  1:4179748/149‑1 (MQ=255)
                                                            cGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCTTTTGTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGAc                                                                                         <  2:3559937/149‑1 (MQ=255)
                                                            cGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTGTTTGATTTGAAAAAATAATGAAAATAATTTTTTTTTAATATAGAGATAAACATAACaaaaaacac                                                                                         >  1:610769/1‑145 (MQ=255)
                                                               aaCCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTTGTTTTTAAAAAAAAAATAAAATATTTTTTTTTTAATATAGAGTTAAAGATCACAAAAAAAACaaa                                                                                      >  1:1227980/1‑148 (MQ=255)
                                                                     aTACCCGCCGCGGTGTGCGTTGTTCAGCATCTGCTTTGTCGATCCCACTTACATTTTATTATCCTGTGGTTTTTTTTCTATTTAAAAATATAATTAAATTTTTTGTTTTTAAATATCGATATAACGATCACAAACACGACAATATGaaa                                                                                <  2:4336617/149‑1 (MQ=255)
                                                                     aTACCCGCCGCGCTGTGCGTTGTACAGAATCTGCTTTGTCGATACCACTTACTTTTTCATATCCTGGCGTTTTTTTTCGATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCTATATAACGAACACAAAAACGACAATATGaaa                                                                                <  2:855943/149‑1 (MQ=255)
                                                                     aTAACCGCGGCGCTGTGCGTTGTTCAGAATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCTTTTGTTTTCTATTTCAAAATATAATTAAATTATTTTTTTTTAATTATCGAGATAACGATCACAAAAACGACAATATGaaa                                                                                <  1:3021166/149‑1 (MQ=255)
                                                                      tAACCGCCGCGGTTTGCGTTTTTCTGAATCTGATTTGTCGATACCACTTACTTTTTATTATCCTGTCTTTTGTTTTCGTTTTCAAAATATAATTAAATTTTTTGTTTTTAAATATCTTTATAACGATCACAAAAACGACAATATGaaaa                                                                               <  1:2261419/149‑1 (MQ=255)
                                                                      tAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTTGGTTTTAAAAAATAATTTAATTATTTTTTTTTATAAATAGTTATAAAGATAAAAAAAACAAAAAAATTaaaa                                                                               >  1:2033270/1‑149 (MQ=255)
                                                                       aaCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAAAAAAAATAAATTATTTGTTTTTTAAAATAGATATAAAGCTAACAAAAACAACAAAATGTAAat                                                                              >  1:3438802/1‑149 (MQ=255)
                                                                        aCCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTGAAAATATAATGTAATTATTTGTTTTTAAATATAGAGATAACGATCACAAAAACGACAATATGtaaagt                                                                             >  2:2934031/1‑147 (MQ=255)
                                                                          cgccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATAAAATGAAAATATTTTTTTTTAAAAATAGAGATAACGATCACAAAAAACACAAAATTAAAaatat                                                                           >  2:2109211/1‑149 (MQ=255)
                                                                           gccgcGGTGTGGGTTGTTCAGCATCTGCTTTGTCTATACCCCGTACATTTTCATATCGTGTCGTTTGTTTTCGATTTCAAAATATAATGAATTTATTTGTTTTTAATTATCGAGATAACGATCACAAAAACGACAATATGAAAattatt                                                                          <  2:1226494/149‑1 (MQ=255)
                                                                           gccgcGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATAATCAAATTATTTTTTTTTAAATATCGAGATAACGATCAAAAAAACGACAATATGAAAattatt                                                                          <  2:3563079/149‑1 (MQ=255)
                                                                               cgcTGTGCGTTGTTCAGCATCTGCTTTGGCGATACCAGTTACATTTTATTATCCTGCCGTTTTTTTGCGTTTTCAAAATATATTGAAATTATTTGTTTTTAAATATCGATATACCGATCACAAACACGACAATATGAAAATTATTCgag                                                                      <  2:2765321/149‑1 (MQ=255)
                                                                                  tgtgGGTTGTTCTTGATGTGCTTTGTCGATACCACGTAATTTTTCATATCCTTTGGTTTGTTTTCTATTTCAAAATATTATGAAATTTTTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCgaggag                                                                   <  2:4235587/149‑1 (MQ=255)
                                                                                  tgtgCGTTGTACAGCATCTGCTTTGTCGATACCACGTACATTTTCATAGCCTGTCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTTTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCgaggag                                                                   <  1:3419646/149‑1 (MQ=255)
                                                                                              agaaTCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGTTATAACGATCACAACCACGACAATATGAAAATTATTCGAGGAGTGAAAGGCaaaa                                                       <  1:1927591/146‑1 (MQ=255)
                                                                                               gCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTAAAAATATAAATAAATTATTTTTTTTTAAATATAGAGATAACGATAAAAAAAAAGAAAAAATTAAAATTTTTCGAGGAGGGAAAGTCaaaaa                                                      >  2:2432224/1‑149 (MQ=255)
                                                                                                    tGCTTTGTCGATACCACGTACTTTTTCATATCCTGTCGTTTGTTTTCTATTTCAAAATATATTTAAATTTTTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACgg                                                 <  2:1765552/149‑1 (MQ=255)
                                                                                                      cTTTGTCGATACCACGTACATTTTCATATCCTGTCGTGTTTTTTGGAATTGAAAAAAAAAATAAAAAATTTTTTTTTTAAAATATAGATAAACATAAAAAAAAACAAAAAATGAAAAATATTAGAGGATTTTAATGCaaaaaaaaagac                                               >  1:1059375/1‑144 (MQ=255)
                                                                                                       tttGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATTAAATTTTTTTTTTTTATATATAGTGATTAAGATAACAAAAACGTCAAAATTAAAATTATTGGGGGTGTGGAATGCAAAAAATCGGCCt                                              >  2:3104953/1‑149 (MQ=255)
                                                                                                          gTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTAAAAATATAATGAAATTATTTTTTTTTTAATATCGAGATAAAGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAAGGCCTccc                                           >  2:378893/1‑149 (MQ=255)
                                                                                                           tCGATACCACGTACATTTTCATATCCTGTCGTGTGTTTTAGATTTCAAAATATAATGAAATTATTTTTTTTTTAATATCGAGTTAACGATCCAAAAAACGACAATATGAAAATTATTAGGGGTGGGGAAGGCCCAAAAAAGGccgcgcg                                          >  2:2422573/1‑144 (MQ=255)
                                                                                                                  cACGTACATTTTCATATCCTGTCGTTTGTTTTAGATTTCAAAATATAATTAAAATATTTTTTTTTAAATATGGAGATAAACATCACAAAAAAGACAATATTAAAATTATTAGATGGGTGAAAGGCAAAAAAAAGGCCTCCCGCTAGGGa                                   >  2:2909345/1‑149 (MQ=255)
                                                                                                                    cGTACATTTTAATTTCATTTGTTTTTTTTTCTTTTAAAAATTATAATTAATTTTTTTTTTTTTAATTATCGTGACATGGATAACAAACAAGAAAATATGAAAATTATTCGAGGAGTTAAAGTCAAAAAAACGGCCTCCCGATAGGGAAg                                 <  1:3727348/149‑1 (MQ=255)
                                                                                                                                  tCCTTTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgc                   <  2:1080438/149‑1 (MQ=255)
                                                                                                                                  tCCTGTCTTTTGTTTTCGATTACAAAATATAATTAAATTTTTTGTTTTTAAATATCGAGATAACTATCACAAACACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTgc                   <  1:3865640/149‑1 (MQ=255)
                                                                                                                                     tGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAATAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                >  2:736125/1‑149 (MQ=255)
                                                                                                                                     tGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGcgcg                >  2:738169/1‑149 (MQ=255)
                                                                                                                                       tCGTTTGTTTTCTATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCgtg              <  2:488543/149‑1 (MQ=255)
                                                                                                                                       tCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGATATAACGATCACAAAAAAGACAATATGAAAATTATTGGTGGAGTGAAAGGCAAAAATAAGGCCTCCCGATAGGGAAGCCGTAGCTATGTGCGCgtg              >  2:2575567/1‑149 (MQ=255)
                                                                                                                                         gtttgtttTCGATTACAAAATATAATTAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt            <  2:904190/149‑1 (MQ=255)
                                                                                                                                          tttgtttTCTATTACAAAATATAATGAAATTATTTGTTTTTAAATATCGATATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt           <  1:3112268/149‑1 (MQ=255)
                                                                                                                                          tttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGttt           >  2:524316/1‑149 (MQ=255)
                                                                                                                                              ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt       >  2:3268397/1‑149 (MQ=255)
                                                                                                                                              ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt       >  2:931805/1‑149 (MQ=255)
                                                                                                                                              ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt       >  2:2185520/1‑149 (MQ=255)
                                                                                                                                                ttCGATTTAAAAATATAATGAATTTATTTTTTTTTAAATATCGATATAACGATAACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc     <  1:2717597/149‑1 (MQ=255)
                                                                                                                                                   gATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGTCAATATGCAAATTATTCGAGGAGTGAAAGGCAAAAAATCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCgg  >  1:305600/1‑149 (MQ=255)
                                                                                                                                                   gATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCgg  >  2:2995167/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CAGGTTGGCAGCGGTTTCTACCACCACTTGCATCGTTTCGAGATTGTGAATATTGAATGCCGGAACCGCATAACCGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGG  >  CP009273/2170561‑2170856

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