Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,534,644 C→A 6.2% S175Y (TCT→TAT)  fimC → periplasmic chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,534,6440CA 6.2% 801.7 / 27.1 227S175Y (TCT→TAT) fimCperiplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base C (109/104);  new base A (3/11);  total (112/115)
Fisher's exact test for biased strand distribution p-value = 5.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 246 total aligned reads displayed.

GATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAA  >  CP009273/4534495‑4534793
                                                                                                                                                     |                                                                                                                                                     
gATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATtc                                                                                                                                                       <  1:4096833/150‑1 (MQ=255)
gATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATtc                                                                                                                                                       <  1:2197397/150‑1 (MQ=255)
gATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATtc                                                                                                                                                       <  1:4217956/150‑1 (MQ=255)
 aTTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATtct                                                                                                                                                      <  2:1960191/150‑1 (MQ=255)
   tCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTctct                                                                                                                                                    >  1:4544244/1‑150 (MQ=255)
   tCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTctct                                                                                                                                                    <  2:4519566/150‑1 (MQ=255)
     cGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGa                                                                                                                                                  <  1:1867586/150‑1 (MQ=255)
     cGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGa                                                                                                                                                  <  2:1997068/150‑1 (MQ=255)
       tCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCCGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACg                                                                                                                                                >  2:18843/1‑150 (MQ=255)
        cAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGc                                                                                                                                               >  1:3292100/1‑150 (MQ=255)
        cAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGc                                                                                                                                               >  2:2655157/1‑150 (MQ=255)
         aaTGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCt                                                                                                                                              <  2:2257474/150‑1 (MQ=255)
           tGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATGTCGTCGTAGCGCGAATTCTCTGACGCTGa                                                                                                                                            >  2:2686780/1‑150 (MQ=255)
            ggATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAACTGAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGAt                                                                                                                                           <  1:957867/150‑1 (MQ=255)
            ggATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGAt                                                                                                                                           <  1:1027973/150‑1 (MQ=255)
             gATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGAtt                                                                                                                                          >  2:1127618/1‑150 (MQ=255)
             gATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGAtt                                                                                                                                          <  1:1195503/150‑1 (MQ=255)
                aaaTCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAAc                                                                                                                                       >  1:3394335/1‑150 (MQ=255)
                aaaTCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAAc                                                                                                                                       >  1:1138497/1‑150 (MQ=255)
                  aTCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAAccc                                                                                                                                     <  2:690887/150‑1 (MQ=255)
                  aTCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAAccc                                                                                                                                     <  1:5059025/150‑1 (MQ=255)
                   tCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCg                                                                                                                                    >  1:2306177/1‑150 (MQ=255)
                    cAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGa                                                                                                                                   <  1:4945556/150‑1 (MQ=255)
                    cAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGa                                                                                                                                   >  1:944517/1‑150 (MQ=255)
                       aaTTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGacac                                                                                                                                <  2:3019876/150‑1 (MQ=255)
                        aTTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACAcc                                                                                                                               >  1:700664/1‑150 (MQ=255)
                          tGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCt                                                                                                                             >  1:3751535/1‑150 (MQ=255)
                              tGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTa                                                                                                                         >  1:4481838/1‑150 (MQ=255)
                                 gaATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCt                                                                                                                      <  1:3159943/150‑1 (MQ=255)
                                 gaATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCt                                                                                                                      <  1:647230/150‑1 (MQ=255)
                                     aCGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACg                                                                                                                  >  1:2233517/1‑150 (MQ=255)
                                       gCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGt                                                                                                                >  1:865447/1‑150 (MQ=255)
                                       gCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGt                                                                                                                <  1:4229881/150‑1 (MQ=255)
                                       gCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGt                                                                                                                <  1:1982433/150‑1 (MQ=255)
                                          aCAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAAc                                                                                                             >  2:2951759/1‑150 (MQ=255)
                                           cAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACa                                                                                                            >  1:2874389/1‑150 (MQ=255)
                                                    aTTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAAt                                                                                                   <  2:27236/150‑1 (MQ=255)
                                                    aTTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAAt                                                                                                   <  1:322880/150‑1 (MQ=255)
                                                    aTTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAAt                                                                                                   >  1:294162/1‑150 (MQ=255)
                                                        tCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCg                                                                                               >  2:68178/1‑150 (MQ=255)
                                                         cAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCgg                                                                                              <  2:869421/150‑1 (MQ=255)
                                                            ccGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAAc                                                                                           >  2:4957670/1‑150 (MQ=255)
                                                               cATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCg                                                                                        <  2:2893021/150‑1 (MQ=255)
                                                                aTTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGACCCCgg                                                                                       >  2:1258346/1‑150 (MQ=255)
                                                                aTTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCgg                                                                                       >  1:1168631/1‑150 (MQ=255)
                                                                 ttAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTATCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGACCCCggg                                                                                      >  2:4355734/1‑150 (MQ=255)
                                                                  tAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGt                                                                                     >  1:1372234/1‑150 (MQ=255)
                                                                    aaCTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTc                                                                                   <  1:3692339/150‑1 (MQ=255)
                                                                      cTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGACCCCGGGTTCtt                                                                                 >  2:4786565/1‑150 (MQ=255)
                                                                      cTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCtt                                                                                 >  1:4318351/1‑150 (MQ=255)
                                                                         tACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGaa                                                                              <  2:4891878/150‑1 (MQ=255)
                                                                         tACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGaa                                                                              >  1:98605/1‑150 (MQ=255)
                                                                            tATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAAt                                                                           <  2:17555/150‑1 (MQ=255)
                                                                             aTCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATg                                                                          >  2:3688348/1‑150 (MQ=255)
                                                                              tCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGc                                                                         >  1:2534290/1‑150 (MQ=255)
                                                                              tCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGc                                                                         >  2:1452389/1‑150 (MQ=255)
                                                                              tCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGc                                                                         >  1:2800502/1‑150 (MQ=255)
                                                                                gCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCAt                                                                       <  1:4444458/150‑1 (MQ=255)
                                                                                 cccGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCAtt                                                                      <  2:1067396/150‑1 (MQ=255)
                                                                                   cGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGACAATGCATTgg                                                                    >  1:4646923/1‑150 (MQ=255)
                                                                                   cGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTgg                                                                    >  1:4646974/1‑150 (MQ=255)
                                                                                       tAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGcc                                                                >  1:5087591/1‑150 (MQ=255)
                                                                                       tAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGcc                                                                <  2:5016893/150‑1 (MQ=255)
                                                                                         aaTTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTc                                                              <  1:2126848/150‑1 (MQ=255)
                                                                                         aaTTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTATCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTc                                                              >  1:2986728/1‑150 (MQ=255)
                                                                                           ttAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCa                                                            <  2:1598862/150‑1 (MQ=255)
                                                                                                gTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATggg                                                       >  1:386031/1‑150 (MQ=255)
                                                                                                      aCCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAg                                                 >  1:537436/1‑150 (MQ=255)
                                                                                                          gATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACg                                             >  1:2367946/1‑150 (MQ=255)
                                                                                                          gATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACg                                             <  1:2555241/150‑1 (MQ=255)
                                                                                                          gATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACg                                             >  1:3677756/1‑150 (MQ=255)
                                                                                                          gATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACg                                             >  1:1466262/1‑150 (MQ=255)
                                                                                                           aTCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACgg                                            <  2:797943/150‑1 (MQ=255)
                                                                                                           aTCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACgg                                            <  2:4391675/150‑1 (MQ=255)
                                                                                                              aGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTa                                         >  1:1877031/1‑150 (MQ=255)
                                                                                                                   gCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTg                                    <  1:5089246/150‑1 (MQ=255)
                                                                                                                       aaaaaTTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCtt                                >  1:4542015/1‑150 (MQ=255)
                                                                                                                        aaaaTTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTc                               <  1:1718241/150‑1 (MQ=255)
                                                                                                                        aaaaTTAAGATTTCGTCGTAGCGCGAATTCTCAGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGATGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGACAGCACGGTTAAATTGACTTc                               >  2:3073057/1‑150 (MQ=255)
                                                                                                                         aaaTTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCt                              >  1:1001548/1‑150 (MQ=255)
                                                                                                                         aaaTTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCt                              <  2:2013723/150‑1 (MQ=255)
                                                                                                                          aaTTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTg                             >  1:3000704/1‑150 (MQ=255)
                                                                                                                          aaTTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTg                             >  1:2076540/1‑150 (MQ=255)
                                                                                                                           aTTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGa                            <  2:4115352/150‑1 (MQ=255)
                                                                                                                           aTTAAGATTTCGTCGTAGCGCGAATTATCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGa                            <  1:1013798/150‑1 (MQ=255)
                                                                                                                              aaGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGGAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGc                         >  2:2227907/1‑150 (MQ=255)
                                                                                                                              aaGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGc                         <  2:3724668/150‑1 (MQ=255)
                                                                                                                              aaGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGc                         <  2:2367435/150‑1 (MQ=255)
                                                                                                                                   ttCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGaa                    >  2:2386947/1‑150 (MQ=255)
                                                                                                                                   ttCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTAGTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGaa                    <  2:3259968/150‑1 (MQ=255)
                                                                                                                                           aGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATAtt            <  2:339113/150‑1 (MQ=255)
                                                                                                                                           aGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATAtt            >  1:1563679/1‑150 (MQ=255)
                                                                                                                                              gcgAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATAttact         <  1:4271081/150‑1 (MQ=255)
                                                                                                                                               cgAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATAttactt        >  2:4192154/1‑150 (MQ=255)
                                                                                                                                                gAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATAttactta       <  2:551299/150‑1 (MQ=255)
                                                                                                                                                  aTTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTAATTAcc     >  2:254855/1‑150 (MQ=255)
                                                                                                                                                  aTTATCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTAcc     <  1:1513113/150‑1 (MQ=255)
                                                                                                                                                   ttCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCg    >  2:2441702/1‑150 (MQ=255)
                                                                                                                                                   ttCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCg    >  1:5186290/1‑150 (MQ=255)
                                                                                                                                                     ctctGACGCTGATTAACCCGACACGCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCATGAAAATGCATTGGTGCCTCCAATGGTCGAAATCACGGTTAAATTGCCTTCTGATGCTGGAATCAATATTTCTTAACGaa  <  2:4791610/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATTCCGTCAATGGATAAATCAAAATTGACTGAGAATACGCTACAGCTCGCAATTATCAGCCGCATTAAACTGTACTATCGCCCGGCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAA  >  CP009273/4534495‑4534793

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: