Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,240,361 C→A 5.9% intergenic (‑113/+181) cirA ← / ← lysP colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,240,3610CA 5.9% 916.2 / 23.9 237intergenic (‑113/+181)cirA/lysPcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter
Reads supporting (aligned to +/- strand):  ref base C (106/117);  new base A (7/7);  total (113/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 243 total aligned reads displayed.

CAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCG  >  CP009273/2240213‑2240509
                                                                                                                                                    |                                                                                                                                                    
caGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatca                                                                                                                                                     >  2:2850072/1‑150 (MQ=255)
caGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatca                                                                                                                                                     >  2:679539/1‑150 (MQ=255)
    ccGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATccc                                                                                                                                                 >  1:3253160/1‑150 (MQ=255)
    ccGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATccc                                                                                                                                                 <  2:3792116/150‑1 (MQ=255)
      gACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAg                                                                                                                                               >  2:39748/1‑150 (MQ=255)
      gACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAg                                                                                                                                               >  1:4149843/1‑150 (MQ=255)
      gACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAg                                                                                                                                               >  1:4149909/1‑150 (MQ=255)
          cGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATAATCCCAGAAAg                                                                                                                                           >  2:3034949/1‑150 (MQ=255)
            tACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGcc                                                                                                                                         >  1:1394693/1‑150 (MQ=255)
             aCGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCg                                                                                                                                        >  2:1570132/1‑150 (MQ=255)
              cGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCgg                                                                                                                                       <  1:2445057/150‑1 (MQ=255)
              cGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCgg                                                                                                                                       >  2:3235781/1‑150 (MQ=255)
                aaaGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGaa                                                                                                                                     >  1:1274798/1‑150 (MQ=255)
                     gTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGcc                                                                                                                                >  2:20109/1‑150 (MQ=255)
                      ttCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCct                                                                                                                               >  1:4494715/1‑150 (MQ=255)
                       tCAACCTAAACATTCCATATCACCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCctc                                                                                                                              <  2:783359/150‑1 (MQ=255)
                           ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCGGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGAt                                                                                                                          <  1:1634829/150‑1 (MQ=255)
                             tAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGAttt                                                                                                                        >  2:4234686/1‑150 (MQ=255)
                               aaCATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGt                                                                                                                      >  2:4169149/1‑150 (MQ=255)
                                aCATTCCATATCTCCATGAGGTAACTACGAAAATAAACTGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTc                                                                                                                     <  2:3492607/150‑1 (MQ=255)
                                     ccATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCt                                                                                                                >  2:1009092/1‑150 (MQ=255)
                                     ccATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCt                                                                                                                <  1:2597587/150‑1 (MQ=255)
                                       atatCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATAATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCttt                                                                                                              <  1:1963762/150‑1 (MQ=255)
                                        tatCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCtttt                                                                                                             >  1:4107978/1‑150 (MQ=255)
                                         atCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCttttt                                                                                                            <  1:1973270/150‑1 (MQ=255)
                                           ctcCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTg                                                                                                          >  1:2760497/1‑150 (MQ=255)
                                           ctcCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATAATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTg                                                                                                          <  1:3918945/150‑1 (MQ=255)
                                            tcCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGa                                                                                                         >  2:239944/1‑150 (MQ=255)
                                              cATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGAtt                                                                                                       <  1:4203239/150‑1 (MQ=255)
                                               aTGAGGTAACTACGAAAATAAAAGGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTg                                                                                                      <  2:3398903/150‑1 (MQ=255)
                                                tGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTgc                                                                                                     <  1:12085/150‑1 (MQ=255)
                                                 gAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTgcg                                                                                                    >  1:2048743/1‑150 (MQ=255)
                                                   ggTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTgcggc                                                                                                  >  1:327967/1‑150 (MQ=255)
                                                     tAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGt                                                                                                >  1:4262428/1‑150 (MQ=255)
                                                      aaCTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTa                                                                                               >  2:4706674/1‑150 (MQ=255)
                                                       aCTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTaa                                                                                              <  1:3034295/150‑1 (MQ=255)
                                                        cTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAAc                                                                                             <  2:917257/150‑1 (MQ=255)
                                                         tACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACg                                                                                            >  2:4314317/1‑150 (MQ=255)
                                                         tACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGAGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACg                                                                                            >  2:2728791/1‑150 (MQ=255)
                                                           cGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGAt                                                                                          <  2:2637593/150‑1 (MQ=255)
                                                           cGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGAt                                                                                          <  1:1570132/150‑1 (MQ=255)
                                                            gAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATa                                                                                         <  2:230599/150‑1 (MQ=255)
                                                            gAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATa                                                                                         >  1:3767785/1‑150 (MQ=255)
                                                            gAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATa                                                                                         <  1:4540163/150‑1 (MQ=255)
                                                            gAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATa                                                                                         >  1:858323/1‑150 (MQ=255)
                                                              aaataaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTGATTGCGGCGTAACGATAATg                                                                                      <  1:840437/150‑1 (MQ=255)
                                                                ataaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTGATTGCGGCGTAACGATAATGCa                                                                                    >  1:3076687/1‑150 (MQ=255)
                                                                 taaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCa                                                                                    <  1:3342038/150‑1 (MQ=255)
                                                                  aaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCaa                                                                                   >  1:4759856/1‑150 (MQ=255)
                                                                    aatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAAc                                                                                 >  1:1111653/1‑150 (MQ=255)
                                                                    aatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAAc                                                                                 >  1:3359918/1‑150 (MQ=255)
                                                                       gggTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATAATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGAt                                                                              >  2:1235467/1‑150 (MQ=255)
                                                                         gTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATaa                                                                            <  1:1470190/150‑1 (MQ=255)
                                                                         gTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATaa                                                                            <  1:4003752/150‑1 (MQ=255)
                                                                          ttATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAAc                                                                           <  2:3880689/150‑1 (MQ=255)
                                                                           tATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACa                                                                          >  2:2321395/1‑150 (MQ=255)
                                                                             tCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAAt                                                                        >  1:288189/1‑150 (MQ=255)
                                                                               gCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTa                                                                      >  1:657173/1‑150 (MQ=255)
                                                                               gCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTa                                                                      >  2:3449981/1‑150 (MQ=255)
                                                                               gCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTa                                                                      >  2:1812788/1‑150 (MQ=255)
                                                                                cTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATAATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTAt                                                                     >  1:359624/1‑150 (MQ=255)
                                                                                cTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATAATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTGTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTAt                                                                     >  1:2188391/1‑150 (MQ=255)
                                                                                 tCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATc                                                                    >  2:1602214/1‑150 (MQ=255)
                                                                                   acaTCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCaa                                                                  <  1:4486059/150‑1 (MQ=255)
                                                                                   acaTCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCaa                                                                  >  1:685080/1‑150 (MQ=255)
                                                                                   acaTCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCaa                                                                  >  2:1452210/1‑150 (MQ=255)
                                                                                    caTCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAAt                                                                 >  2:1942518/1‑150 (MQ=255)
                                                                                    caTCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAAt                                                                 >  2:2787644/1‑150 (MQ=255)
                                                                                        ttcttcACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCaa                                                             >  2:3714141/1‑150 (MQ=255)
                                                                                          cttcACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAAt                                                           <  2:886225/150‑1 (MQ=255)
                                                                                            tcACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGt                                                         <  1:2666423/150‑1 (MQ=255)
                                                                                            tcACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGt                                                         <  2:1523156/150‑1 (MQ=255)
                                                                                            tcACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGt                                                         <  1:4314317/150‑1 (MQ=255)
                                                                                              aCGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTa                                                       <  2:3527916/150‑1 (MQ=255)
                                                                                                   cccTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAAttt                                                  <  2:3324239/150‑1 (MQ=255)
                                                                                                    ccTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAAtttt                                                 <  2:4149843/150‑1 (MQ=255)
                                                                                                    ccTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAAtttt                                                 <  2:4149909/150‑1 (MQ=255)
                                                                                                     cTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTta                                                <  2:2584056/150‑1 (MQ=255)
                                                                                                       tGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTtata                                              <  1:3182013/150‑1 (MQ=255)
                                                                                                          gTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATcc                                           <  2:178564/150‑1 (MQ=255)
                                                                                                          gTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATcc                                           <  2:1120919/150‑1 (MQ=255)
                                                                                                             cGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATcgcat                                        >  1:1035643/1‑146 (MQ=255)
                                                                                                             cGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCa                                        >  1:1035713/1‑150 (MQ=255)
                                                                                                              ggCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAGCGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCAt                                       <  1:793454/150‑1 (MQ=255)
                                                                                                              ggCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCAt                                       <  2:2576309/150‑1 (MQ=255)
                                                                                                                 agcaTCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGAt                                    <  2:1833003/150‑1 (MQ=255)
                                                                                                                   caTCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGAttt                                  >  2:1666107/1‑150 (MQ=255)
                                                                                                                      cgcgGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTc                               <  1:4680043/150‑1 (MQ=255)
                                                                                                                         ggTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAgg                            <  1:2835810/150‑1 (MQ=255)
                                                                                                                           tGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGa                          >  2:822058/1‑150 (MQ=255)
                                                                                                                              aCTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAAt                       <  1:1429002/150‑1 (MQ=255)
                                                                                                                              aCTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAAt                       <  1:3166049/150‑1 (MQ=255)
                                                                                                                                ttATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATgg                     >  2:1557274/1‑150 (MQ=255)
                                                                                                                                ttATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATgg                     <  1:436448/150‑1 (MQ=255)
                                                                                                                                ttATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATgg                     <  1:3908560/150‑1 (MQ=255)
                                                                                                                                 tATTTTTTATGCAGGTGATAATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGa                    <  2:4302855/150‑1 (MQ=255)
                                                                                                                                   ttttttATGCAGGTGATAATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAATCGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGact                  >  1:2566203/1‑148 (MQ=255)
                                                                                                                                        tATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGaaa             >  2:2520912/1‑150 (MQ=255)
                                                                                                                                             aGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGcc        >  1:2056545/1‑150 (MQ=255)
                                                                                                                                                tGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCctc     <  2:4494715/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCCGGAACAGCCTCTGATTTGTCAGCCTTTTTTGATTGCGGCGTAACGATAATGCAAACGATAACAATTATCAATCCAAATGTTAAATTTTATATCCGCATGATTTGTCAGGGAAATGGAATTGAAAAAGCCCTCTCG  >  CP009273/2240213‑2240509

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: