Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,506,430 A→T 8.3% L22Q (CTG→CAG)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,506,4300AT 8.3% 1019.7 / 63.8 290L22Q (CTG→CAG) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (93/173);  new base T (12/12);  total (105/185)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 311 total aligned reads displayed.

CGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTC  >  CP009273/4506281‑4506579
                                                                                                                                                     |                                                                                                                                                     
cGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGgca                                                                                                                                                       <  2:4672665/150‑1 (MQ=255)
cGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGgca                                                                                                                                                       >  2:3044349/1‑150 (MQ=255)
 gTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGgcag                                                                                                                                                      >  1:356424/1‑150 (MQ=255)
   tGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGgcagca                                                                                                                                                    <  1:2265214/150‑1 (MQ=255)
        tGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTc                                                                                                                                               >  2:751280/1‑150 (MQ=255)
         gCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCa                                                                                                                                              >  1:1719924/1‑150 (MQ=255)
          cTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAACTCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAg                                                                                                                                             <  2:2814877/150‑1 (MQ=255)
          cTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAg                                                                                                                                             <  2:63440/150‑1 (MQ=255)
          cTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAg                                                                                                                                             <  2:2939575/150‑1 (MQ=255)
          cTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAg                                                                                                                                             <  2:2228348/150‑1 (MQ=255)
           ttGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAgg                                                                                                                                            <  2:772742/150‑1 (MQ=255)
           ttGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAgg                                                                                                                                            <  2:246770/150‑1 (MQ=255)
            tGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGc                                                                                                                                           <  1:1054769/150‑1 (MQ=255)
             gCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGCGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCg                                                                                                                                          <  1:2002039/150‑1 (MQ=255)
              ccGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGa                                                                                                                                         <  1:3057994/150‑1 (MQ=255)
              ccGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGa                                                                                                                                         <  1:3057633/150‑1 (MQ=255)
               cgcgcgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGaa                                                                                                                                        <  2:1322195/150‑1 (MQ=255)
               cgcgcgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGaa                                                                                                                                        <  2:4100789/150‑1 (MQ=255)
               cgcgcgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGaa                                                                                                                                        <  2:417300/150‑1 (MQ=255)
               cgcgcgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGaa                                                                                                                                        >  2:558321/1‑150 (MQ=255)
                 cgcgTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATg                                                                                                                                      <  2:3972379/150‑1 (MQ=255)
                      cAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGtt                                                                                                                                 <  1:3780189/150‑1 (MQ=255)
                      cAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGtt                                                                                                                                 <  1:796118/150‑1 (MQ=255)
                      cAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGtt                                                                                                                                 >  1:613444/1‑150 (MQ=255)
                         gCTGGCGTCAACGGTTAGGTTAAATCCGCTGTGTGCGGCATACTGATGGTGCGCTTTATCGCGCGATCCCGGTGCGATTTTAACGTGTGCAGCAAAAGCGTAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAAc                                                                                                                              <  1:1186487/150‑1 (MQ=255)
                         gCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAAc                                                                                                                              <  1:345119/150‑1 (MQ=255)
                          cTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAAcc                                                                                                                             >  2:3420196/1‑150 (MQ=255)
                           tGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCa                                                                                                                            <  1:2625996/150‑1 (MQ=255)
                                     cGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGtgttgt                                                                                                                  >  2:608405/1‑150 (MQ=255)
                                       agagGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGttt                                                                                                                <  2:3850827/150‑1 (MQ=255)
                                       agagGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGttt                                                                                                                >  2:4136240/1‑150 (MQ=255)
                                        gagGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGtttt                                                                                                               <  1:2585202/150‑1 (MQ=255)
                                            gTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGa                                                                                                           <  2:1652436/150‑1 (MQ=255)
                                                aTCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGaaaaa                                                                                                       <  1:4193410/150‑1 (MQ=255)
                                                 tCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAAc                                                                                                      <  2:1245094/150‑1 (MQ=255)
                                                  ccGCTGTGTGCGGCATACTGATTGAGAGATTAATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGAAAAAGCGGAAAAAGAGAGACCGGCCAGAGGCAGAAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAAct                                                                                                     >  2:532220/1‑149 (MQ=255)
                                                      tgtgtgCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTa                                                                                                 <  2:3596961/150‑1 (MQ=255)
                                                      tgtgtgCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTa                                                                                                 <  1:499938/150‑1 (MQ=255)
                                                       gtgtgCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTaa                                                                                                >  1:4028117/1‑150 (MQ=255)
                                                          tgCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACgg                                                                                             <  2:4185812/150‑1 (MQ=255)
                                                          tgCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACgg                                                                                             >  1:3519440/1‑150 (MQ=255)
                                                            cGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCg                                                                                           >  1:3323651/1‑150 (MQ=255)
                                                             ggCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGt                                                                                          >  2:2432410/1‑150 (MQ=255)
                                                               cATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCa                                                                                        <  1:4687665/150‑1 (MQ=255)
                                                                aTACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCAt                                                                                       <  1:2181030/150‑1 (MQ=255)
                                                                aTACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCAt                                                                                       >  1:3694540/1‑150 (MQ=255)
                                                                 tACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATa                                                                                      <  1:2117410/150‑1 (MQ=255)
                                                                 tACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATa                                                                                      <  2:1123861/150‑1 (MQ=255)
                                                                  aCTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATAc                                                                                     >  1:2759127/1‑150 (MQ=255)
                                                                                  aTCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCAtt                                                                     <  2:4752129/150‑1 (MQ=255)
                                                                                    cGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCAtttt                                                                   <  1:751280/150‑1 (MQ=255)
                                                                                          aTCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTtgttgt                                                             <  2:1378824/150‑1 (MQ=255)
                                                                                          aTCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTtgttgt                                                             >  2:1867198/1‑150 (MQ=255)
                                                                                           tCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTtgtg                                                            >  1:1317225/1‑150 (MQ=255)
                                                                                           tCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTtgtg                                                            >  1:2493587/1‑150 (MQ=255)
                                                                                           tCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTtgtg                                                            >  1:1684354/1‑150 (MQ=255)
                                                                                            cccGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTgtgt                                                           <  2:3205967/150‑1 (MQ=255)
                                                                                               ggTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCa                                                        >  1:4250904/1‑150 (MQ=255)
                                                                                               ggTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCa                                                        <  1:4568175/150‑1 (MQ=255)
                                                                                                   cGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTa                                                    <  2:3156689/150‑1 (MQ=255)
                                                                                                         tAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTc                                              >  2:139050/1‑150 (MQ=255)
                                                                                                          aaCCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCg                                             <  2:1580216/150‑1 (MQ=255)
                                                                                                          aaCCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCg                                             <  2:3079878/150‑1 (MQ=255)
                                                                                                          aaCCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCg                                             <  1:2441632/150‑1 (MQ=255)
                                                                                                             cTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAAc                                          <  1:3164262/150‑1 (MQ=255)
                                                                                                              tgtgCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACg                                         >  1:4341304/1‑150 (MQ=255)
                                                                                                               gtgCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACga                                        >  2:2695567/1‑150 (MQ=255)
                                                                                                                tgCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACgag                                       >  1:1543768/1‑150 (MQ=255)
                                                                                                                 gcagcaAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACgaga                                      <  2:3344929/150‑1 (MQ=255)
                                                                                                                  cagcaAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGaa                                     >  2:1485103/1‑150 (MQ=255)
                                                                                                                   agcaAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAg                                    <  1:614355/150‑1 (MQ=255)
                                                                                                                   agcaAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAg                                    >  2:1447368/1‑150 (MQ=255)
                                                                                                                    gcaAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAgg                                   <  2:2668474/150‑1 (MQ=255)
                                                                                                                     caAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAggg                                  <  1:3920488/150‑1 (MQ=255)
                                                                                                                       aaaGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTa                                <  2:4595059/150‑1 (MQ=255)
                                                                                                                        aaGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTaa                               >  1:4492866/1‑150 (MQ=255)
                                                                                                                          gCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTaaaa                             >  2:1013212/1‑150 (MQ=255)
                                                                                                                           cGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTaaaaa                            <  1:1865306/150‑1 (MQ=255)
                                                                                                                            ggAAAACGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAg                           >  2:4440777/1‑150 (MQ=255)
                                                                                                                            ggAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAg                           >  2:1523880/1‑150 (MQ=255)
                                                                                                                             gAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAgg                          <  1:927029/150‑1 (MQ=255)
                                                                                                                              aaaaCGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGa                         <  2:2521096/150‑1 (MQ=255)
                                                                                                                              aaaaCGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGa                         >  2:1078317/1‑150 (MQ=255)
                                                                                                                              aaaaCGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGa                         <  2:4260482/150‑1 (MQ=255)
                                                                                                                               aaaCGAGAGACCGGCCAGCGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGAc                        >  2:2193938/1‑150 (MQ=255)
                                                                                                                                     gagaCCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTGTTGTTGTGTTCAGCTATGAGTCGAACGAGACGGGTAAAAAGGACAATCGa                  >  2:3015412/1‑150 (MQ=255)
                                                                                                                                        aCCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAAt               <  1:2936551/150‑1 (MQ=255)
                                                                                                                                          cGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATaa             <  1:1592795/150‑1 (MQ=255)
                                                                                                                                             ccAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGaa          <  1:558321/150‑1 (MQ=255)
                                                                                                                                               aGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAatt        >  1:1754947/1‑150 (MQ=255)
                                                                                                                                                gCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAatta       >  1:4240295/1‑150 (MQ=255)
                                                                                                                                                gCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAatta       <  1:2878271/150‑1 (MQ=255)
                                                                                                                                                 cGGCTGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAattat      <  1:2334880/150‑1 (MQ=255)
                                                                                                                                                  ggCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAattatt     <  1:1139515/150‑1 (MQ=255)
                                                                                                                                                   gcagcaGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTAttt    <  2:4388386/150‑1 (MQ=255)
                                                                                                                                                   gcagcaGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTAttt    <  2:1951669/150‑1 (MQ=255)
                                                                                                                                                   gcagcaGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTAttt    <  2:1952138/150‑1 (MQ=255)
                                                                                                                                                   gcagcaGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTAttt    <  2:2054363/150‑1 (MQ=255)
                                                                                                                                                     agcagGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTc  <  2:1759741/150‑1 (MQ=255)
                                                                                                                                                     agcagGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTc  >  2:2742663/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTTGCTCTGCTTGCCGCGCGTCAGGCTGGCGTCAACCGAGAGGGTAAATCCGCTGTGTGCGGCATACTGATTGAGCGCTTTATCGAGCGATCCCGGTGCGATATTAACCTGTGCAGCAAAAGCGGAAAACGAGAGACCGGCCAGCGGCAGCAGGCTCAGGCGAATGGTGTTAACCAAAGGTGTTGTTTTACGAAAAACGCGTAACGGCGTCATACCTTCCCCATCATCATTTTTGTTGTGTTCAGCTATGAGTCGAACGAGAAGGGTAAAAAGGACAATCGAAATAAGAATTATTTTC  >  CP009273/4506281‑4506579

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: