Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 607,660 G→T 8.2% R97R (CGC→CGA)  fepA ← iron‑enterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273607,6600GT 8.2% 1136.2 / 58.8 318R97R (CGC→CGA) fepAiron‑enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (159/133);  new base T (10/16);  total (169/149)
Fisher's exact test for biased strand distribution p-value = 1.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

Only 100 of 324 total aligned reads displayed.

AGAACTTCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAACGCCAGG  >  CP009273/607511‑607809
                                                                                                                                                     |                                                                                                                                                     
 gAACTTCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCAGGACCCATACcgc                                                                                                                                                      >  2:4856571/1‑150 (MQ=255)
  aaCTTCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACcgcg                                                                                                                                                     <  2:934194/150‑1 (MQ=255)
   aCTTCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGa                                                                                                                                                    >  2:2283509/1‑150 (MQ=255)
     ttCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAat                                                                                                                                                  <  2:4983003/150‑1 (MQ=255)
     ttCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAat                                                                                                                                                  <  1:3712875/150‑1 (MQ=255)
     ttCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAat                                                                                                                                                  <  1:3658972/150‑1 (MQ=255)
       cAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCCATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAatat                                                                                                                                                >  2:2042046/1‑150 (MQ=255)
       cAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATAAGCGTGTTTTCAGGACCCATACCGCGAatat                                                                                                                                                >  2:477432/1‑150 (MQ=255)
        aaTACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATc                                                                                                                                               >  1:419078/1‑150 (MQ=255)
         aTACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCa                                                                                                                                              >  2:4129507/1‑150 (MQ=255)
         aTACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCa                                                                                                                                              <  2:2004521/150‑1 (MQ=255)
           aagTTCAATCATTTCAGGTGGCACCCTGGAAGTATCACCTCGGGTATCGCGCTCGCCAAGCCAGCCCTGAAGCACCGAGTTAAGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAAt                                                                                                                                            <  2:3482647/148‑1 (MQ=255)
            cGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCTCGAATATCAATc                                                                                                                                           <  1:5017517/150‑1 (MQ=255)
                    cATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAAAATCAATCTGTCGGtt                                                                                                                                   >  2:4954003/1‑150 (MQ=255)
                          aGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTccc                                                                                                                             <  1:4973849/150‑1 (MQ=255)
                          aGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTccc                                                                                                                             >  2:1021190/1‑150 (MQ=255)
                           ggTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCa                                                                                                                            <  2:1708247/150‑1 (MQ=255)
                              ggCACACAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGc                                                                                                                         >  1:705123/1‑150 (MQ=255)
                               gCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCt                                                                                                                        >  2:2514796/1‑150 (MQ=255)
                                 aCCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGa                                                                                                                      >  1:2024308/1‑150 (MQ=255)
                                  cccAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGAc                                                                                                                     >  2:4221520/1‑150 (MQ=255)
                                    cAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCa                                                                                                                   >  1:4752187/1‑150 (MQ=255)
                                    cAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCa                                                                                                                   <  1:4906199/150‑1 (MQ=255)
                                    cAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTGATTCCCACGCTGAccc                                                                                                                   >  2:1262021/1‑149 (MQ=255)
                                        aaGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACtgg                                                                                                               <  1:2845603/150‑1 (MQ=255)
                                            aTCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGa                                                                                                           >  2:2562372/1‑150 (MQ=255)
                                            aTCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGa                                                                                                           >  2:854648/1‑150 (MQ=255)
                                             tCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGATATTCCCACGCTGACCACTGGTGGAg                                                                                                          >  1:3049793/1‑150 (MQ=255)
                                                  acGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTAcc                                                                                                     <  1:751972/150‑1 (MQ=255)
                                                  acGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTAcc                                                                                                     <  1:1497501/150‑1 (MQ=255)
                                                    gggTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCgg                                                                                                   >  2:2660836/1‑150 (MQ=255)
                                                    gggTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCgg                                                                                                   >  1:1267144/1‑150 (MQ=255)
                                                    gggTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCgg                                                                                                   >  1:5178261/1‑150 (MQ=255)
                                                      gTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTc                                                                                                 <  2:717911/150‑1 (MQ=255)
                                                       tATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCa                                                                                                >  1:3235721/1‑150 (MQ=255)
                                                       tATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCa                                                                                                >  2:472649/1‑150 (MQ=255)
                                                       tATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCa                                                                                                >  2:2441197/1‑150 (MQ=255)
                                                       tATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCa                                                                                                <  2:4661722/150‑1 (MQ=255)
                                                        aTCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAg                                                                                               <  1:1145365/150‑1 (MQ=255)
                                                         tCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAgg                                                                                              >  2:1377842/1‑150 (MQ=255)
                                                          cgcgCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGt                                                                                             <  1:1713186/150‑1 (MQ=255)
                                                          cgcgCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGt                                                                                             <  1:1896213/150‑1 (MQ=255)
                                                          cgcgCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGt                                                                                             <  1:3017702/150‑1 (MQ=255)
                                                          cgcgCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGt                                                                                             <  1:3873766/150‑1 (MQ=255)
                                                            cgcTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGACTATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGGTACCGGTCAGGTTa                                                                                           >  2:3103082/1‑150 (MQ=255)
                                                            cgcTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTa                                                                                           >  1:513284/1‑150 (MQ=255)
                                                             gcTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTaa                                                                                          >  2:3012968/1‑150 (MQ=255)
                                                             gcTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTaa                                                                                          <  2:510308/150‑1 (MQ=255)
                                                                 gccacgccaGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGcc                                                                                      >  2:4265613/1‑150 (MQ=255)
                                                                 gccacgccaGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGcc                                                                                      >  2:1128613/1‑150 (MQ=255)
                                                                     cgccaGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGc                                                                                  >  2:158831/1‑150 (MQ=255)
                                                                       ccaGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCAt                                                                                >  2:3912013/1‑150 (MQ=255)
                                                                        caGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATg                                                                               <  2:357014/150‑1 (MQ=255)
                                                                             cTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACg                                                                          <  1:973748/150‑1 (MQ=255)
                                                                               gACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGga                                                                        <  2:3466431/150‑1 (MQ=255)
                                                                                aCGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGgat                                                                       <  2:3164902/150‑1 (MQ=255)
                                                                                aCGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGgat                                                                       <  1:4368853/150‑1 (MQ=255)
                                                                                   cACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGTTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGgatgat                                                                    >  1:4885454/1‑150 (MQ=255)
                                                                                   cACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGgatgat                                                                    >  1:2080463/1‑150 (MQ=255)
                                                                                   cACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGgatgat                                                                    >  1:1378392/1‑150 (MQ=255)
                                                                                   cACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGgatgat                                                                    >  1:4714000/1‑150 (MQ=255)
                                                                                   cACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGgatgat                                                                    <  2:2183963/150‑1 (MQ=255)
                                                                                     ccGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCt                                                                  <  1:3848829/150‑1 (MQ=255)
                                                                                          ttACGGCTGCTTACCGGCTTGCCGTCAATCAACATCAGCGTGTTTTCCGGACCCCTACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGac                                                             <  1:1029142/150‑1 (MQ=255)
                                                                                           tACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTAGTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGaca                                                            <  2:4000852/150‑1 (MQ=255)
                                                                                            aCGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGacac                                                           >  2:2417233/1‑150 (MQ=255)
                                                                                              ggCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACAt                                                         <  2:4541973/150‑1 (MQ=255)
                                                                                              ggCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACAt                                                         <  1:1610530/150‑1 (MQ=255)
                                                                                              ggCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACAt                                                         <  1:3621534/150‑1 (MQ=255)
                                                                                                   ctTACCGGCTTGCCGTCAAACAAAATCAGCGTGTTTTCCGGACCCATACCGCGAAAATCAATCTGTCGGTTATTCCAAAGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATATTCGAAACATcgcgt                                                    >  2:4399540/1‑149 (MQ=255)
                                                                                                    ttACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCggg                                                   >  2:1826085/1‑150 (MQ=255)
                                                                                                    ttACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCggg                                                   >  1:692223/1‑150 (MQ=255)
                                                                                                     tACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGc                                                  >  1:3679108/1‑150 (MQ=255)
                                                                                                      aCCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCAGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCa                                                 >  1:2352541/1‑150 (MQ=255)
                                                                                                        cGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAAc                                               >  2:2852302/1‑150 (MQ=255)
                                                                                                          gCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCg                                             >  1:2638780/1‑150 (MQ=255)
                                                                                                           cTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCgg                                            <  2:659805/150‑1 (MQ=255)
                                                                                                             tGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGct                                                 <  2:3720588/143‑2 (MQ=255)
                                                                                                             tGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGt                                          >  1:1425984/1‑150 (MQ=255)
                                                                                                             tGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGt                                          >  1:4036747/1‑150 (MQ=255)
                                                                                                              gCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGtt                                         <  1:4129507/150‑1 (MQ=255)
                                                                                                               ccGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGttt                                        >  2:1619860/1‑150 (MQ=255)
                                                                                                                cGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGtttt                                       >  1:3206139/1‑150 (MQ=255)
                                                                                                                 gTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGttttt                                      >  2:5182551/1‑150 (MQ=255)
                                                                                                                     atcaaAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCgg                                  <  1:1982114/150‑1 (MQ=255)
                                                                                                                     atcaaAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTAGTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCgg                                  <  2:5063073/150‑1 (MQ=255)
                                                                                                                         aaaTCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGTATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGAttt                              <  1:933080/150‑1 (MQ=255)
                                                                                                                             cAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATc                          >  1:3457805/1‑150 (MQ=255)
                                                                                                                              aGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCt                         >  1:2353083/1‑150 (MQ=255)
                                                                                                                               gCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTg                        <  2:279119/150‑1 (MQ=255)
                                                                                                                                  tgtTTTCCGGACCCATACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCgg                     >  1:4290230/1‑150 (MQ=255)
                                                                                                                                   gtTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGt                    >  1:2403196/1‑150 (MQ=255)
                                                                                                                                      ttCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGAt                 <  1:2748617/150‑1 (MQ=255)
                                                                                                                                       tCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATCGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATg                <  1:1177831/150‑1 (MQ=255)
                                                                                                                                         cGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGt              <  2:268958/150‑1 (MQ=255)
                                                                                                                                           gACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCg            >  1:594556/1‑150 (MQ=255)
                                                                                                                                               cATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAAc        <  1:4629903/150‑1 (MQ=255)
                                                                                                                                                  aCCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAACGcc     >  1:3565204/1‑150 (MQ=255)
                                                                                                                                                  aCCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAACGcc     >  1:2859677/1‑150 (MQ=255)
                                                                                                                                                   ccGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAACGCCa    >  2:3248911/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGAACTTCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAACGCCAGG  >  CP009273/607511‑607809

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: ‑