Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 619,048 C→A 44.9% A307S (GCC→TCC)  fepB ← iron‑enterobactin transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273619,0480CA44.9% 222.9 / inf 296A307S (GCC→TCC) fepBiron‑enterobactin transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (97/66);  new base A (56/77);  total (153/143)
Fisher's exact test for biased strand distribution p-value = 3.45e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

CGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATTTA  >  CP009273/618900‑619197
                                                                                                                                                    |                                                                                                                                                     
cgAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCAATCCAGCACTTGCATGgcg                                                                                                                                                      >  2:3012770/1‑150 (MQ=255)
  aGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCt                                                                                                                                                    >  1:1549269/1‑150 (MQ=255)
         ttGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTaa                                                                                                                                             >  2:1027065/1‑150 (MQ=255)
          tGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAAt                                                                                                                                            >  2:59095/1‑150 (MQ=255)
           ggTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATc                                                                                                                                           <  2:3434511/150‑1 (MQ=255)
           ggTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATc                                                                                                                                           >  2:315628/1‑150 (MQ=255)
            gTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATcc                                                                                                                                          >  2:1369710/1‑150 (MQ=255)
             tttGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCa                                                                                                                                         >  1:3185712/1‑150 (MQ=255)
                gttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGAc                                                                                                                                      <  1:527358/150‑1 (MQ=255)
                gttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGAc                                                                                                                                      <  1:2386555/150‑1 (MQ=255)
                      tATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACg                                                                                                                                <  2:940346/150‑1 (MQ=255)
                        tCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGtc                                                                                                                              >  1:2324121/1‑150 (MQ=255)
                               gttgttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTc                                                                                                                       >  2:3510319/1‑150 (MQ=255)
                               gttgttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTc                                                                                                                       >  2:3468590/1‑150 (MQ=255)
                                 tgttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTccc                                                                                                                     <  1:954753/150‑1 (MQ=255)
                                    tATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAgc                                                                                                                  >  1:2511691/1‑150 (MQ=255)
                                     aTTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAgcg                                                                                                                 >  2:3950658/1‑150 (MQ=255)
                                     aTTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAgcg                                                                                                                 >  1:1435068/1‑150 (MQ=255)
                                         ctggtgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATa                                                                                                             >  2:1758468/1‑150 (MQ=255)
                                           ggtgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATaaa                                                                                                           >  2:2507488/1‑150 (MQ=255)
                                             tgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAAcc                                                                                                         <  1:2988757/150‑1 (MQ=255)
                                              gctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCt                                                                                                        >  2:2329237/1‑150 (MQ=255)
                                                tggtggAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGc                                                                                                      >  2:3708043/1‑150 (MQ=255)
                                                  gtggAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCtt                                                                                                    <  1:3946587/150‑1 (MQ=255)
                                                   tggAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTg                                                                                                   <  1:1231056/150‑1 (MQ=255)
                                                   tggAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTg                                                                                                   >  1:1690090/1‑150 (MQ=255)
                                                       gTTGCGACATTTTCGCCAGACGCAGCCGCATGTGTCAGCGTCCGACAGTTAATTCTTAAAACAGCGCCTTAAGCCTATACAGCACTTGCATGGCGCTGTAGTAATCCTGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGtttt                                                                                               <  1:3385982/150‑1 (MQ=255)
                                                         tGCGACATTTTCGCCAGACGCCGCCGCTGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGt                                                                                             >  2:2269405/1‑150 (MQ=255)
                                                           cGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTAc                                                                                           <  2:2029022/150‑1 (MQ=255)
                                                            gACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACt                                                                                          >  2:179141/1‑150 (MQ=255)
                                                            gACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACt                                                                                          >  1:3589361/1‑150 (MQ=255)
                                                            gACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACt                                                                                          >  1:3834993/1‑150 (MQ=255)
                                                            gACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACt                                                                                          <  2:3757200/150‑1 (MQ=255)
                                                              cATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTgc                                                                                        >  2:4614483/1‑150 (MQ=255)
                                                              cATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTgc                                                                                        >  2:3359791/1‑150 (MQ=255)
                                                                   tCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTgcaggca                                                                                   <  1:280157/150‑1 (MQ=255)
                                                                     gCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGcaggcagg                                                                                 >  1:2945497/1‑150 (MQ=255)
                                                                      ccAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGt                                                                                <  1:1795896/150‑1 (MQ=255)
                                                                        aGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTgc                                                                              >  2:3851119/1‑150 (MQ=255)
                                                                         gACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTgcg                                                                             <  1:163984/150‑1 (MQ=255)
                                                                            gccgccgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGa                                                                          <  1:1133266/150‑1 (MQ=255)
                                                                              cgccgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGagc                                                                        <  2:1726758/150‑1 (MQ=255)
                                                                                    aGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGa                                                                  <  2:1431727/150‑1 (MQ=255)
                                                                                     ggTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGAt                                                                 <  1:2434377/150‑1 (MQ=255)
                                                                                       tGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTa                                                               >  1:3630360/1‑150 (MQ=255)
                                                                                        gACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAg                                                              >  2:494722/1‑150 (MQ=255)
                                                                                        gACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAg                                                              <  1:1093580/150‑1 (MQ=255)
                                                                                        gACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAg                                                              <  2:718998/150‑1 (MQ=255)
                                                                                            gCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATa                                                          <  1:4672306/150‑1 (MQ=255)
                                                                                                  gACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGc                                                    <  1:835221/150‑1 (MQ=255)
                                                                                                  gACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGc                                                    >  1:4564559/1‑150 (MQ=255)
                                                                                                     aGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAta                                                 <  2:1309825/150‑2 (MQ=255)
                                                                                                      gTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCg                                                <  1:473657/150‑1 (MQ=255)
                                                                                                      gTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCg                                                <  2:1673285/150‑1 (MQ=255)
                                                                                                       ttAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCgg                                               <  2:2824691/150‑1 (MQ=255)
                                                                                                         aaTGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCa                                             >  1:2057702/1‑150 (MQ=255)
                                                                                                         aaTGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCa                                             >  2:4771423/1‑150 (MQ=255)
                                                                                                         aaTGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCa                                             <  1:1940578/150‑1 (MQ=255)
                                                                                                          aTGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCAt                                            >  2:1817045/1‑150 (MQ=255)
                                                                                                            gCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCt                                          <  2:2730141/150‑1 (MQ=255)
                                                                                                            gCTTAAAACAGCGCCTTAAGCCTATCCAGCACGTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCt                                          <  2:4283/150‑1 (MQ=255)
                                                                                                              ttAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCttt                                        <  1:3310982/150‑1 (MQ=255)
                                                                                                              ttAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCttt                                        >  2:3301215/1‑150 (MQ=255)
                                                                                                              ttAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCttt                                        >  2:1149737/1‑150 (MQ=255)
                                                                                                              ttAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCttt                                        <  2:1091351/150‑1 (MQ=255)
                                                                                                                aaaaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCt                                      <  1:484087/150‑1 (MQ=255)
                                                                                                                aaaaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCt                                      >  2:446862/1‑150 (MQ=255)
                                                                                                                  aaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCACAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGa                                    <  1:3798657/150‑1 (MQ=255)
                                                                                                                     aGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCa                                 <  2:1930194/150‑1 (MQ=255)
                                                                                                                       cgcCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCAcc                               >  1:1366753/1‑150 (MQ=255)
                                                                                                                        gcCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCg                              >  1:865081/1‑150 (MQ=255)
                                                                                                                          cTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGc                            <  1:3967633/150‑1 (MQ=255)
                                                                                                                           ttAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCg                           <  2:4182771/150‑1 (MQ=255)
                                                                                                                            tAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGa                          >  2:1386618/1‑150 (MQ=255)
                                                                                                                             aaGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGaa                         >  2:287461/1‑150 (MQ=255)
                                                                                                                              aGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAAc                        >  1:2575827/1‑150 (MQ=255)
                                                                                                                              aGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAAc                        >  1:4777444/1‑150 (MQ=255)
                                                                                                                               gCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACa                       >  1:1007306/1‑150 (MQ=255)
                                                                                                                                ccTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAg                      <  2:419895/150‑1 (MQ=255)
                                                                                                                                 cTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGACCAgg                     >  1:674221/1‑150 (MQ=255)
                                                                                                                                 cTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAgg                     >  1:3063857/1‑150 (MQ=255)
                                                                                                                                  tATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGa                    >  2:100027/1‑150 (MQ=255)
                                                                                                                                    tCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAAt                  >  2:976075/1‑150 (MQ=255)
                                                                                                                                     ccAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATa                 <  2:2278411/150‑1 (MQ=255)
                                                                                                                                      cAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAg                >  1:1200136/1‑150 (MQ=255)
                                                                                                                                      cAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAg                >  1:350995/1‑150 (MQ=255)
                                                                                                                                        gCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTg              <  1:1582377/150‑1 (MQ=255)
                                                                                                                                         cACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGa             >  2:3974517/1‑150 (MQ=255)
                                                                                                                                         cACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGa             <  2:4337620/150‑1 (MQ=255)
                                                                                                                                          aCTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGAc            <  1:315628/150‑1 (MQ=255)
                                                                                                                                            ttGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTc          >  2:2067418/1‑150 (MQ=255)
                                                                                                                                               cATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGcc       >  1:4139903/1‑150 (MQ=255)
                                                                                                                                               cATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGcc       <  1:4378543/150‑1 (MQ=255)
                                                                                                                                                aTGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCa      <  1:1010638/150‑1 (MQ=255)
                                                                                                                                                aTGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCa      >  1:4280292/1‑150 (MQ=255)
                                                                                                                                                 tGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCAt     <  1:2269405/150‑1 (MQ=255)
                                                                                                                                                 tGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCAt     <  2:1507143/150‑1 (MQ=255)
                                                                                                                                                 tGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCAt     <  2:4088535/150‑1 (MQ=255)
                                                                                                                                                 tGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCAt     <  2:3916139/150‑1 (MQ=255)
                                                                                                                                                    cgcTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATTTa  <  1:1758468/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCCATTTA  >  CP009273/618900‑619197

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: