Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,238,283 C→A 100% intergenic (‑183/‑180) uxaC ← / → exuT uronate isomerase/hexuronate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,238,2830CA98.6% 940.9 / 2.5 290intergenic (‑183/‑180)uxaC/exuTuronate isomerase/hexuronate transporter
Reads supporting (aligned to +/- strand):  ref base C (0/4);  new base A (147/139);  total (147/143)
Fisher's exact test for biased strand distribution p-value = 5.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 291 total aligned reads displayed.

TGCGATGCTGACAAAGTTATCACACCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTCAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGATTCATTCATTTGTTAGATGAATCGGG  >  CP009273/3238141‑3238425
                                                                                                                                              |                                                                                                                                              
tGCGATGCTGACAAATTTATCACACCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTa                                                                                                                                           <  2:2498171/148‑1 (MQ=255)
      gCTGACAAAGTTATCACACCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTcc                                                                                                                                     >  1:2837128/1‑148 (MQ=255)
       cTGACAAAGTTATCACACCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTccc                                                                                                                                    <  2:1734583/148‑1 (MQ=255)
       cTGACAAAGTTATCACACCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTccc                                                                                                                                    >  1:1594538/1‑148 (MQ=255)
        tGACAAAGTTATCACACCAATATCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTCAGGTAATGTCCCt                                                                                                                                   <  1:32598/148‑1 (MQ=255)
         gACAAAGTTATCACACCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTa                                                                                                                                  >  1:3033767/1‑148 (MQ=255)
         gACAAAGTTATCACACCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTa                                                                                                                                  >  1:1176182/1‑148 (MQ=255)
              aGTTATCACACCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAat                                                                                                                             <  1:3135895/148‑1 (MQ=255)
                ttATCACACCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAatat                                                                                                                           >  1:2445208/1‑148 (MQ=255)
                ttATCACACCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAatat                                                                                                                           <  1:2659506/148‑1 (MQ=255)
                   tCACACCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTcc                                                                                                                        <  1:3236650/148‑1 (MQ=255)
                    cacaCCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTccc                                                                                                                       >  1:1577870/1‑148 (MQ=255)
                     acacCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCca                                                                                                                      >  2:2060028/1‑148 (MQ=255)
                      cacCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTCATGTCCCTACAAATATTCCcac                                                                                                                     >  2:3256588/1‑148 (MQ=255)
                      cacCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCcac                                                                                                                     <  2:43312/148‑1 (MQ=255)
                       acCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCcaca                                                                                                                    <  2:624902/148‑1 (MQ=255)
                         cAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACAtt                                                                                                                  <  1:828146/148‑1 (MQ=255)
                               ccAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGAt                                                                                                            <  2:542515/148‑1 (MQ=255)
                                cAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATg                                                                                                           <  1:2525992/148‑1 (MQ=255)
                                  gagTCAGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGc                                                                                                         >  2:981512/1‑148 (MQ=255)
                                    gTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGctc                                                                                                       <  2:76516/148‑1 (MQ=255)
                                     tCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGctct                                                                                                      <  1:1260256/148‑1 (MQ=255)
                                        gAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCAc                                                                                                   <  1:3729301/148‑1 (MQ=255)
                                         aaGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCAcc                                                                                                  >  1:1505532/1‑148 (MQ=255)
                                          aGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCt                                                                                                 >  2:2444712/1‑148 (MQ=255)
                                           gATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCtt                                                                                                >  2:1022813/1‑148 (MQ=255)
                                            atatTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCttt                                                                                               >  1:3683570/1‑148 (MQ=255)
                                            atatTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCttt                                                                                               <  1:1921134/148‑1 (MQ=255)
                                            atatTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCttt                                                                                               >  1:1793262/1‑148 (MQ=255)
                                              atTTTCGTGAGTTAGATCAATAAACGTAGTTATAAAAATGACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTa                                                                                             <  2:2612130/148‑1 (MQ=255)
                                              atTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTa                                                                                             >  2:2276912/1‑148 (MQ=255)
                                               ttttCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTaa                                                                                            >  1:1330918/1‑148 (MQ=255)
                                               ttttCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTaa                                                                                            <  1:3207281/148‑1 (MQ=255)
                                                tttCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTTAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTaaa                                                                                           <  1:1011778/148‑1 (MQ=255)
                                                tttCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTTAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTaaa                                                                                           <  1:3469809/148‑1 (MQ=255)
                                                tttCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTaaa                                                                                           >  1:2974438/1‑148 (MQ=255)
                                                tttCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTaaa                                                                                           <  2:1847834/148‑1 (MQ=255)
                                                     tGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGt                                                                                      <  2:2200387/148‑1 (MQ=255)
                                                         ttAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGa                                                                                  >  1:2878721/1‑148 (MQ=255)
                                                         ttAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATACTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGa                                                                                  >  1:1127653/1‑148 (MQ=255)
                                                          tAGATCAATAAACGTAGTTAAAAAAATTAATCTCAAAGTGGTAAATCTCGCTGCAGGCCGCTCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCATACAAATATTCACACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGAc                                                                                 >  1:69120/1‑148 (MQ=255)
                                                           aGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACa                                                                                <  1:2795140/148‑1 (MQ=255)
                                                            gATCAATAAACGTAGTTAAAAAAATTGCTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACaa                                                                               >  1:1995099/1‑148 (MQ=255)
                                                             aTCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAg                                                                              <  1:3104829/148‑1 (MQ=255)
                                                             aTCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAg                                                                              >  2:3176802/1‑148 (MQ=255)
                                                               cAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGtt                                                                            >  1:291236/1‑148 (MQ=255)
                                                                 aTAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTAt                                                                          >  2:3668207/1‑148 (MQ=255)
                                                                  tAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATc                                                                         >  2:1584865/1‑148 (MQ=255)
                                                                  tAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATc                                                                         >  1:1199618/1‑148 (MQ=255)
                                                                  tAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATc                                                                         >  1:2168258/1‑148 (MQ=255)
                                                                  tAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATc                                                                         >  1:1063275/1‑148 (MQ=255)
                                                                   aaaCGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCt                                                                        <  2:3033767/148‑1 (MQ=255)
                                                                     aCGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCttt                                                                      <  1:1417209/148‑1 (MQ=255)
                                                                      cGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTc                                                                     >  2:2006267/1‑148 (MQ=255)
                                                                      cGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTc                                                                     >  2:1848269/1‑148 (MQ=255)
                                                                           ttAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCATGCTGTCGTAATGTAATGTCACTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTTTATGACAAGTTATCTTTCTGCCg                                                                >  1:3410291/1‑148 (MQ=255)
                                                                               aaaaaTTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGc                                                            >  1:1647762/1‑148 (MQ=255)
                                                                               aaaaaTTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGc                                                            >  1:1138152/1‑148 (MQ=255)
                                                                                 aaaTTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATTTCCCTACAAATATTCCCACATATGTGATGGATCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCaa                                                          >  1:1946467/1‑148 (MQ=255)
                                                                                  aaTTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCaaa                                                         <  1:679410/148‑1 (MQ=255)
                                                                                   aTTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAAt                                                        >  2:2441829/1‑148 (MQ=255)
                                                                                    ttACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATc                                                       >  1:3495715/1‑148 (MQ=255)
                                                                                     tACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCa                                                      <  1:466652/148‑1 (MQ=255)
                                                                                       ctctCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATa                                                    >  1:3095030/1‑148 (MQ=255)
                                                                                        tctcAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCTTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATaa                                                   <  2:1404547/148‑1 (MQ=255)
                                                                                        tctcAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATaa                                                   >  2:1431226/1‑148 (MQ=255)
                                                                                        tctcAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATaa                                                   >  1:1485007/1‑148 (MQ=255)
                                                                                          tcAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGt                                                 >  2:907005/1‑148 (MQ=255)
                                                                                            aaaGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCg                                               <  2:3642075/148‑1 (MQ=255)
                                                                                             aaGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGa                                              >  2:415682/1‑148 (MQ=255)
                                                                                                 ggTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGa                                          >  2:3512598/1‑148 (MQ=255)
                                                                                                 ggTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGACGTCGCAAATCATAAGTCGACGGa                                          >  1:2515225/1‑148 (MQ=255)
                                                                                                  gTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGaa                                         <  2:254735/148‑1 (MQ=255)
                                                                                                   tAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATACGTCGACGGAAt                                        >  2:130988/1‑148 (MQ=255)
                                                                                                   tAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAAt                                        >  1:2867167/1‑148 (MQ=255)
                                                                                                   tAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTGTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATCAGTCGACGGAAt                                        >  2:1087355/1‑148 (MQ=255)
                                                                                                    aaaTCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATg                                       <  2:3192619/148‑1 (MQ=255)
                                                                                                       tctcGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCaa                                    >  1:3608523/1‑148 (MQ=255)
                                                                                                        ctcGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTCAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCaaa                                   <  1:3793085/148‑1 (MQ=255)
                                                                                                          cGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAAtt                                 >  2:3117971/1‑148 (MQ=255)
                                                                                                           gCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTg                                <  2:1615887/148‑1 (MQ=255)
                                                                                                            cTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGc                               <  2:1767009/148‑1 (MQ=255)
                                                                                                            cTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGc                               <  2:3659039/148‑1 (MQ=255)
                                                                                                               cAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGa                            <  1:2222878/148‑1 (MQ=255)
                                                                                                               cAGGCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGa                            >  1:3768928/1‑148 (MQ=255)
                                                                                                                  gCCGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGattc                         <  2:3363537/148‑1 (MQ=255)
                                                                                                                   ccGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGattca                        <  1:1865505/148‑1 (MQ=255)
                                                                                                                   ccGCGCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGattca                        <  1:2125834/148‑1 (MQ=255)
                                                                                                                    cgcgCCAGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGattcat                       <  1:2170857/148‑1 (MQ=255)
                                                                                                                          aGTACTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGATTCATTCAttt                 <  2:3743099/148‑1 (MQ=255)
                                                                                                                             acgggCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCCCCTTTGTGATGGCTCTCCCCTTTTAAAGTTGTCTGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGATTCATTCATTTGtt              <  2:2295346/145‑1 (MQ=255)
                                                                                                                             aCTGGCCTTGCTGTCTTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGATTCATTCATTTGtt              <  1:2603494/148‑1 (MQ=255)
                                                                                                                             aCTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGATTCATTCATTTGtt              <  1:1848262/148‑1 (MQ=255)
                                                                                                                             aCTGGCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGATTCATTCATTTGtt              <  1:3321393/148‑1 (MQ=255)
                                                                                                                                ggCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGATTCATTCATTTGTTAGa           <  1:2060028/148‑1 (MQ=255)
                                                                                                                                 gCCTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGATTCATTCATTTGTTAGAt          >  2:1702662/1‑148 (MQ=255)
                                                                                                                                   cTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGATTCATTCATTTGTTAGATGa        <  2:1946636/148‑1 (MQ=255)
                                                                                                                                   cTTGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAACCATAAGTCGACGGAATGCAAATTGCCGATTCATTCATTTGTTAGATGa        <  2:1159679/148‑1 (MQ=255)
                                                                                                                                    ttGCTGTCGTAAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGATTCATTCATTTGTTAGATGaa       <  1:1035706/148‑1 (MQ=255)
                                                                                                                                         gtcgtaAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGATTCATTCATTTGTTAGATGAATCggg  <  1:32771/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
TGCGATGCTGACAAAGTTATCACACCAATTTCCAGAGTCCGAAGATATTTTCGTGAGTTAGATCAATAAACGTAGTTAAAAAAATTACTCTCAAAGTGGTAAATCTCGCTGCAGGCCGCGCCAGTACTGGCCTTGCTGTCGTCAGGTAATGTCCCTACAAATATTCCCACATTTGTGATGGCTCTCACCTTTTAAAGTTGTATGACAAGTTATCTTTCTGCCGTCGCAAATCATAAGTCGACGGAATGCAAATTGCCGATTCATTCATTTGTTAGATGAATCGGG  >  CP009273/3238141‑3238425

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: