Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 739,056 G→T 15.4% R3L (CGA→CTA)  ybgJ → putative allophanate hydrolase, subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273739,0560GT15.4% 447.3 / 46.6 165R3L (CGA→CTA) ybgJputative allophanate hydrolase, subunit 1
Reads supporting (aligned to +/- strand):  ref base G (88/50);  new base T (19/7);  total (107/58)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 165 total aligned reads displayed.

TTTTATGCTGCAGGTCACCATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCGGCATGAACA  >  CP009273/738910‑739202
                                                                                                                                                  |                                                                                                                                                  
ttttATGCTGCAGGTCACCATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGa                                                                                                                                                   >  1:516645/1‑148 (MQ=255)
 tttATGCTGCAGGTCACCATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAg                                                                                                                                                  <  2:2893378/148‑1 (MQ=255)
     tgctgcAGGTCACCATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACTAGcgcg                                                                                                                                              <  2:12214/148‑1 (MQ=255)
     tgctgcAGGTCACCATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGcgcg                                                                                                                                              >  1:466381/1‑148 (MQ=255)
          cAGGTCACCATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCgttgtt                                                                                                                                         >  2:2114665/1‑148 (MQ=255)
             gTCACCATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATc                                                                                                                                      >  2:716388/1‑148 (MQ=255)
              tCACCATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACTAGCGCGTTGTTATCt                                                                                                                                     <  2:3069123/148‑1 (MQ=255)
                 ccATGCCACTGATCGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGAt                                                                                                                                  <  2:1959602/148‑1 (MQ=255)
                 ccATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGAt                                                                                                                                  >  2:2719877/1‑148 (MQ=255)
                 ccATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGAt                                                                                                                                  <  1:2293438/148‑1 (MQ=255)
                 ccATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGAt                                                                                                                                  <  1:2292926/148‑1 (MQ=255)
                  cATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATa                                                                                                                                 >  1:801285/1‑148 (MQ=255)
                  cATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATa                                                                                                                                 <  2:822404/148‑1 (MQ=255)
                         cTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGaaa                                                                                                                          <  1:1255307/148‑1 (MQ=255)
                          tGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACTAGCGCGTTGTTATCTGATAGGTGAAAc                                                                                                                         >  1:2794839/1‑148 (MQ=255)
                          tGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACTAGCGCGTTGTTATCTGATAGGTGAAAc                                                                                                                         >  1:1122976/1‑148 (MQ=255)
                            aaCGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAAcgg                                                                                                                       >  2:2315711/1‑148 (MQ=255)
                             aCGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAAcggc                                                                                                                      >  2:2303885/1‑148 (MQ=255)
                             aCGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAAcggc                                                                                                                      >  2:2969366/1‑148 (MQ=255)
                                 ggtggtATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACTAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTa                                                                                                                  >  2:276077/1‑148 (MQ=255)
                                 ggtggtATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTa                                                                                                                  >  2:1967997/1‑148 (MQ=255)
                                  gtggtATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAg                                                                                                                 >  1:3078665/1‑148 (MQ=255)
                                  gtggtATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAg                                                                                                                 >  1:3078029/1‑148 (MQ=255)
                                    ggtATTCGCGCTTTGAGCGAGTGGCTGAATTAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAATGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTg                                                                                                               <  1:1549306/148‑1 (MQ=255)
                                      tATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGct                                                                                                             >  1:1412626/1‑148 (MQ=255)
                                              cATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGctggaactgg                                                                                                     >  1:2585934/1‑148 (MQ=255)
                                               aTTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACTAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGctggaactgga                                                                                                    >  2:94012/1‑148 (MQ=255)
                                               aTTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGctggaactgga                                                                                                    >  2:1371831/1‑148 (MQ=255)
                                               aTTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGctggaactgga                                                                                                    <  2:154313/148‑1 (MQ=255)
                                                ttGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGctggaactggaa                                                                                                   >  1:1664918/1‑148 (MQ=255)
                                                 tGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCtggaactggaac                                                                                                  >  2:189683/1‑148 (MQ=255)
                                                  gAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAAcc                                                                                                 >  1:1892442/1‑148 (MQ=255)
                                                           gCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCAGGTGAc                                                                                        >  2:3192842/1‑148 (MQ=255)
                                                            cTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATTTGCCTAATCCTGCATAACCAATAATCTGTGGTATCTAAAGTGAAACGAGCGCGTTGTTATCTGATAGGTCAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACg                                                                                       <  2:1681913/148‑1 (MQ=255)
                                                              gaatgaaAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGGTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCt                                                                                     >  1:1070695/1‑148 (MQ=255)
                                                                atgaaAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTgg                                                                                   >  1:1832376/1‑148 (MQ=255)
                                                                atgaaAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTgg                                                                                   >  1:259679/1‑148 (MQ=255)
                                                                atgaaAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTgg                                                                                   >  1:1248879/1‑148 (MQ=255)
                                                                   aaaaTACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTa                                                                                >  1:223337/1‑148 (MQ=255)
                                                                     aaTACCGATCTTGATGTGACCTTTATTGATATTCCTACGCCTCCATAACGAATAATCAGAGGGATCGAAAGTGCAACTCGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGc                                                                              <  2:912103/148‑1 (MQ=255)
                                                                     aaTACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGc                                                                              >  2:658298/1‑148 (MQ=255)
                                                                      aTACCGATCTTGATGTGACCTTTATTGATTTTCCTAATCATGCATACCGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTGGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGcc                                                                             <  2:2463298/148‑1 (MQ=255)
                                                                      aTACCGATCTGGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGcc                                                                             <  2:2463286/148‑1 (MQ=255)
                                                                         ccGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGa                                                                          <  2:2014333/148‑1 (MQ=255)
                                                                          cGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGaa                                                                         >  1:222382/1‑148 (MQ=255)
                                                                           gATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGaaa                                                                        >  2:161681/1‑148 (MQ=255)
                                                                            aTCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACTAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAAc                                                                       >  2:2896533/1‑148 (MQ=255)
                                                                             tCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACTAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACg                                                                      >  1:1920746/1‑148 (MQ=255)
                                                                             tCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACTAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACg                                                                      >  1:612915/1‑148 (MQ=255)
                                                                             tCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACg                                                                      >  1:2822023/1‑148 (MQ=255)
                                                                             tCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACg                                                                      >  1:2822169/1‑148 (MQ=255)
                                                                              cTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACgg                                                                     >  2:1717257/1‑148 (MQ=255)
                                                                                tGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGAt                                                                   >  1:736949/1‑148 (MQ=255)
                                                                                  aTGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACTAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCt                                                                 >  1:3221316/1‑148 (MQ=255)
                                                                                  aTGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCt                                                                 >  2:3098081/1‑148 (MQ=255)
                                                                                   tgtgACCTTTATTTATATTCCTAATCCTGCATAACGAATAATCAGTGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATTGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTg                                                                <  2:1765925/148‑1 (MQ=255)
                                                                                    gtgACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTgg                                                               >  1:1160028/1‑148 (MQ=255)
                                                                                          tttatttaTATTCCTAATCCTGCATAACCAATAATCAGTGGTAGCTAAAGTGCAACGTGCGCTTTTTTATCTGATCGTTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTg                                                         <  1:3247440/148‑1 (MQ=255)
                                                                                            tATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGgc                                                       >  1:530431/1‑148 (MQ=255)
                                                                                              ttGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACTAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGgcgc                                                     >  2:920359/1‑148 (MQ=255)
                                                                                                  tatTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACTAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCg                                                 >  1:3189455/1‑148 (MQ=255)
                                                                                                  tatTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCg                                                 >  2:1855816/1‑148 (MQ=255)
                                                                                                    ttCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTc                                               >  1:3168551/1‑148 (MQ=255)
                                                                                                    ttCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTc                                               >  1:1934216/1‑148 (MQ=255)
                                                                                                        tAATCCTGAATAACGAATAATCAGCGGGATCGAAAGTGCACCGCGCGCGTTGTTATCTGATCGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCtggt                                           <  1:484429/148‑1 (MQ=255)
                                                                                                          aTCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCtggtgg                                         >  2:310751/1‑148 (MQ=255)
                                                                                                            gtttcATACCGAATAATCTGTGGGATCGAATGTGCAACGAGCGGGTTGTTATCTGATTGGTGAAAGGGCGGTAGTGGTGGAAATGGCACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGat                                       <  1:2218219/144‑1 (MQ=255)
                                                                                                            ccTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGat                                       >  1:913090/1‑148 (MQ=255)
                                                                                                            ccTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGat                                       <  2:516645/148‑1 (MQ=255)
                                                                                                               ccaaaccGAATAATCAGAGGGAGCGAATGTCCACCGCGCGCTTTGTAATCTGATCGGGGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATg                                    <  1:102705/142‑1 (MQ=255)
                                                                                                                   aaCGAATAATCAGAGGGATCGAAAGTGCACCGAGCGGGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGa                                <  1:3020415/148‑1 (MQ=255)
                                                                                                                   aaCGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGACCCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGa                                >  1:2836387/1‑148 (MQ=255)
                                                                                                                   aaCGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGa                                >  1:3252521/1‑148 (MQ=255)
                                                                                                                   aaCGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGa                                >  2:1063200/1‑148 (MQ=255)
                                                                                                                   aaCGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGa                                >  1:207322/1‑148 (MQ=255)
                                                                                                                   aaCGAATAAGCAGCGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTTATCGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGa                                <  1:1650915/148‑1 (MQ=255)
                                                                                                                    aCGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGaa                               >  1:1392041/1‑148 (MQ=255)
                                                                                                                     ctaaaaatCAGTGGTATCTAAAGTCAACCGAGGGGTTTTTTATATTATCGGTCAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAAt                              <  1:1536090/143‑1 (MQ=255)
                                                                                                                       aataatCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAAtgt                            >  2:267891/1‑148 (MQ=255)
                                                                                                                        ataatCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAAtgtg                           >  1:1758008/1‑148 (MQ=255)
                                                                                                                         taatCCGTGGTATCTAACGTGCAACCAGCGCGTTTTTATCTTATCGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTgg                          <  1:478864/148‑1 (MQ=255)
                                                                                                                             cAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGa                      >  2:3005965/1‑148 (MQ=255)
                                                                                                                             cAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGa                      >  2:301648/1‑148 (MQ=255)
                                                                                                                                agGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCTTCTGGTGGATATGCCGAATGTGGTTGAAGc                   >  1:2387851/1‑148 (MQ=255)
                                                                                                                                  ggATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCa                 <  2:447976/148‑1 (MQ=255)
                                                                                                                                   gATCGAAAGTGCAACGCGCGCTTTGTTATCTGATTGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCAt                <  1:2142109/148‑1 (MQ=255)
                                                                                                                                    aTCGAAAGTGCAACGAGCGCGTTGTTATCTGATCGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCAtt               <  1:2719877/148‑1 (MQ=255)
                                                                                                                                        aaaGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCg           >  1:272680/1‑148 (MQ=255)
                                                                                                                                        aaaGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCg           >  2:1421964/1‑148 (MQ=255)
                                                                                                                                         aaGTGCACCGCGCGCGTTGTTATGTGATCGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCgg          <  1:2969401/148‑1 (MQ=255)
                                                                                                                                         aaGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCgg          <  1:2969366/148‑1 (MQ=255)
                                                                                                                                         aaGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCgg          >  1:2537429/1‑148 (MQ=255)
                                                                                                                                          aGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCGGc         >  1:1173373/1‑148 (MQ=255)
                                                                                                                                            tGCAACTAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCGGCAt       >  2:1038771/1‑148 (MQ=255)
                                                                                                                                            tGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCGGCAt       >  1:512514/1‑148 (MQ=255)
                                                                                                                                             gCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCGGCATg      <  2:486590/148‑1 (MQ=255)
                                                                                                                                               aaCTAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCGGCATGaa    <  2:2327717/148‑1 (MQ=255)
                                                                                                                                                 cGAGCGGGTTGTTATCTGTTCGGTGAACCGGCGGTTGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCGGCATGAACa  <  2:1122976/148‑1 (MQ=255)
                                                                                                                                                 cGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTTGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCGGCATGAACa  <  1:1291417/148‑1 (MQ=255)
                                                                                                                                                 cGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCGGCATGAACa  >  1:1232501/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTTTATGCTGCAGGTCACCATGCCACTGAACGTGGTGGTATTCGCGCATTGAGCGAGTGGCTGAATGAAAATACCGATCTTGATGTGACCTTTATTGATATTCCTAATCCTGCATAACGAATAATCAGAGGGATCGAAAGTGCAACGAGCGCGTTGTTATCTGATAGGTGAAACGGCGGTAGTGCTGGAACTGGAACCGCCGGTGACGCTGGCTAGCCAGAAACGGATCTGGCGACTGGCGCAGCGTCTGGTGGATATGCCGAATGTGGTTGAAGCCATTCCCGGCATGAACA  >  CP009273/738910‑739202

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