Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,531,143 A→T 100% I124F (ATT→TTT)  fimB → tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,531,1430AT99.6% 889.3 / ‑4.6 278I124F (ATT→TTT) fimBtyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (124/153);  minor base C (1/0);  total (125/153)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 283 total aligned reads displayed.

CAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTC  >  CP009273/4531000‑4531287
                                                                                                                                               |                                                                                                                                                
caacaaCGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCAttttat                                                                                                                                              <  2:1016456/148‑1 (MQ=255)
 aacaacGCACCCGCTATTGAATAAAGAGGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtc                                                                                                                                             <  2:1761272/148‑1 (MQ=255)
 aacaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtc                                                                                                                                             <  2:674633/148‑1 (MQ=255)
  acaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtct                                                                                                                                            >  2:1510749/1‑148 (MQ=255)
  acaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtct                                                                                                                                            <  1:54361/148‑1 (MQ=255)
  acaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTAtct                                                                                                                                            <  2:2931406/148‑1 (MQ=255)
    aacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCg                                                                                                                                          <  1:1458496/148‑1 (MQ=255)
     acGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGa                                                                                                                                         <  2:2614902/148‑1 (MQ=255)
     acGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGa                                                                                                                                         >  1:3276834/1‑148 (MQ=255)
      cGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGAc                                                                                                                                        >  2:2307545/1‑148 (MQ=255)
        cACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACtt                                                                                                                                      <  2:1698938/148‑1 (MQ=255)
        cACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACtt                                                                                                                                      >  1:836429/1‑148 (MQ=255)
         aCCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTc                                                                                                                                     >  2:2205792/1‑148 (MQ=255)
         aCCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTc                                                                                                                                     <  1:3426351/148‑1 (MQ=255)
             gCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCggt                                                                                                                                 >  2:1842108/1‑148 (MQ=255)
               tATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCggtgg                                                                                                                               <  2:3253412/148‑1 (MQ=255)
                  tGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTaa                                                                                                                            <  1:354825/148‑1 (MQ=255)
                    aaTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATg                                                                                                                          <  1:1706212/148‑1 (MQ=255)
                    aaTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATg                                                                                                                          <  1:2475883/148‑1 (MQ=255)
                    aaTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATg                                                                                                                          >  1:1179373/1‑148 (MQ=255)
                    aaTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATg                                                                                                                          >  2:3114254/1‑148 (MQ=255)
                     aTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGc                                                                                                                         <  1:1774850/148‑1 (MQ=255)
                         agaagTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCggg                                                                                                                     <  2:1376163/148‑1 (MQ=255)
                             gTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGt                                                                                                                 >  1:998522/1‑148 (MQ=255)
                             gTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGt                                                                                                                 <  1:2239914/148‑1 (MQ=255)
                              ttCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTc                                                                                                                >  2:3012301/1‑148 (MQ=255)
                              ttCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTc                                                                                                                >  1:611809/1‑148 (MQ=255)
                               tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCCACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCa                                                                                                               >  1:1603727/1‑148 (MQ=255)
                               tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCa                                                                                                               <  1:2834600/148‑1 (MQ=255)
                               tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCa                                                                                                               <  1:2770724/148‑1 (MQ=255)
                               tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCa                                                                                                               <  2:2101325/148‑1 (MQ=255)
                                 aGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACt                                                                                                             <  1:1454737/148‑1 (MQ=255)
                                   gCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTgg                                                                                                           >  2:2182597/1‑148 (MQ=255)
                                           aaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCAt                                                                                                   <  1:1842108/148‑1 (MQ=255)
                                           aaaCTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCAt                                                                                                   <  2:2733631/148‑1 (MQ=255)
                                                  ttGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGcaca                                                                                            >  2:1743498/1‑148 (MQ=255)
                                                     aGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGt                                                                                         <  2:3302177/148‑1 (MQ=255)
                                                      gTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGtt                                                                                        <  2:2624670/148‑1 (MQ=255)
                                                          ccGTACTTAGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGc                                                                                    >  2:3227581/1‑148 (MQ=255)
                                                                  cGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgtg                                                                            <  1:483695/148‑1 (MQ=255)
                                                                  cGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgtg                                                                            <  2:766956/148‑1 (MQ=255)
                                                                   gTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGTCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgtgt                                                                           <  2:3079133/148‑1 (MQ=255)
                                                                   gTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgtgt                                                                           <  1:444034/148‑1 (MQ=255)
                                                                   gTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCgtgt                                                                           <  1:2334865/148‑1 (MQ=255)
                                                                      cccGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGGTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGt                                                                        >  1:3265323/1‑148 (MQ=255)
                                                                       ccGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGtt                                                                       <  2:1057702/148‑1 (MQ=255)
                                                                       ccGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGtt                                                                       <  1:3034754/148‑1 (MQ=255)
                                                                       ccGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGtt                                                                       >  2:912134/1‑148 (MQ=255)
                                                                        cGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGttt                                                                      <  2:2736093/148‑1 (MQ=255)
                                                                           aTGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGc                                                                   >  1:1546471/1‑148 (MQ=255)
                                                                            tGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCt                                                                  <  2:652348/148‑1 (MQ=255)
                                                                            tGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCt                                                                  <  2:1775572/148‑1 (MQ=255)
                                                                              cTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCttt                                                                <  2:1241398/148‑1 (MQ=255)
                                                                                   agCGAGTGGGTATTTTTATCACTTAAGGGTAATCCGCTTTCTCGGCAACACTTTTACCATTTTATATCGACTTGCGGTGGTAATTCCGGGTTGTCACTGGAGATTAATCCGCACATGTTACGCCATTCGTGTGTTTTTGCTTTGgaga                                                           >  2:941015/1‑148 (MQ=255)
                                                                                     cGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTTCTTTGGCGAat                                                         >  1:2694616/1‑148 (MQ=255)
                                                                                       aGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAatat                                                       >  1:3220055/1‑148 (MQ=255)
                                                                                       aGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAatat                                                       >  2:1128693/1‑148 (MQ=255)
                                                                                        gTGGGTATTTTTATCACGTCAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATg                                                      <  2:3159547/148‑1 (MQ=255)
                                                                                        gTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATg                                                      >  2:631513/1‑148 (MQ=255)
                                                                                         tGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATgg                                                     <  2:506265/148‑1 (MQ=255)
                                                                                         tGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATgg                                                     >  2:1927641/1‑148 (MQ=255)
                                                                                          gggTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATggg                                                    >  1:122431/1‑148 (MQ=255)
                                                                                             tATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATCTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAAt                                                 >  2:1896853/1‑148 (MQ=255)
                                                                                              aTTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATa                                                >  2:2132182/1‑148 (MQ=255)
                                                                                               tttttATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAg                                               <  2:3371518/148‑1 (MQ=255)
                                                                                               tttttATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAg                                               <  1:914245/148‑1 (MQ=255)
                                                                                                   tATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATAc                                           <  2:2415370/148‑1 (MQ=255)
                                                                                                    aTCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATAcg                                          <  2:381538/148‑1 (MQ=255)
                                                                                                        cGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGAc                                      >  1:705948/1‑148 (MQ=255)
                                                                                                             ggggAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATc                                 <  1:1005296/148‑1 (MQ=255)
                                                                                                             ggggAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATc                                 >  1:308737/1‑148 (MQ=255)
                                                                                                             ggggAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATc                                 <  1:1004551/148‑1 (MQ=255)
                                                                                                               ggAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCa                               <  2:2547156/148‑1 (MQ=255)
                                                                                                               ggAATCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCa                               <  2:488022/148‑1 (MQ=255)
                                                                                                                 aaTCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAgg                             <  2:3391559/148‑1 (MQ=255)
                                                                                                                  aTCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGa                            >  1:3085866/1‑148 (MQ=255)
                                                                                                                   tCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGAt                           <  1:1061026/148‑1 (MQ=255)
                                                                                                                   tCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGAt                           <  1:3394344/148‑1 (MQ=255)
                                                                                                                   tCCGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGAt                           <  2:3146563/148‑1 (MQ=255)
                                                                                                                     cGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTa                         <  2:1568177/148‑1 (MQ=255)
                                                                                                                     cGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTa                         <  2:59615/148‑1 (MQ=255)
                                                                                                                     cGCTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTa                         <  2:1960483/148‑1 (MQ=255)
                                                                                                                       cTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATc                       <  2:1223105/148‑1 (MQ=255)
                                                                                                                       cTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATc                       >  1:604329/1‑148 (MQ=255)
                                                                                                                       cTTTCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATc                       >  2:1202015/1‑148 (MQ=255)
                                                                                                                        tttCTCGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATACGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCt                      <  1:3434682/148‑1 (MQ=255)
                                                                                                                          tctcGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACGGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTg                    <  2:115654/148‑1 (MQ=255)
                                                                                                                            tcGGCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTggg                  <  2:3466625/148‑1 (MQ=255)
                                                                                                                               gCAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCAt               <  2:617855/148‑1 (MQ=255)
                                                                                                                                cAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATc              >  2:1032633/1‑148 (MQ=255)
                                                                                                                                cAACAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATc              >  1:3479533/1‑148 (MQ=255)
                                                                                                                                  aCAGTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGc            >  1:2004837/1‑148 (MQ=255)
                                                                                                                                     gTTTTACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAat         <  2:201068/148‑1 (MQ=255)
                                                                                                                                      ttttACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAata        >  1:3199125/1‑148 (MQ=255)
                                                                                                                                       tttACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAatat       <  2:291626/148‑1 (MQ=255)
                                                                                                                                       tttACCATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAatat       >  2:148541/1‑148 (MQ=255)
                                                                                                                                        ttACCATTTTCTCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATAtt      <  2:2787504/148‑1 (MQ=255)
                                                                                                                                           ccATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGt   >  2:1216288/1‑148 (MQ=255)
                                                                                                                                            cATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTc  <  2:3276778/148‑1 (MQ=255)
                                                                                                                                            cATTTTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTc  >  1:1887178/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                
CAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTC  >  CP009273/4531000‑4531287

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: