Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 619,048 C→A 100% 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,0480CA97.2% 440.7 / 4.1 142A307S (GCC→TCC) fepBiron‑enterobactin transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base A (66/71);  total (68/74)
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.90e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 142 total aligned reads displayed.

CGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCC  >  CP009273/618904‑619192
                                                                                                                                                |                                                                                                                                                
cGGTTTTGGTTTGTTGATGATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCt                                                                                                                                               <  2:3425149/148‑1 (MQ=255)
   ttttGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgta                                                                                                                                            >  1:1655099/1‑148 (MQ=255)
    tttGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgtag                                                                                                                                           >  2:1126565/1‑148 (MQ=255)
    tttGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTgtag                                                                                                                                           >  2:1126746/1‑148 (MQ=255)
        gTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAAt                                                                                                                                       <  2:1785873/148‑1 (MQ=255)
         tttGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATc                                                                                                                                      >  1:2621844/1‑148 (MQ=255)
         tttGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCATATCCAGCACAAGCATGGAGCTGTAGTAATc                                                                                                                                      >  2:1774502/1‑148 (MQ=255)
         tttGTTGATTATCGGCGTGGTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATc                                                                                                                                      >  1:2513502/1‑148 (MQ=255)
          ttgttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATTCTTAAAAAAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATcc                                                                                                                                     >  1:321767/1‑148 (MQ=255)
          ttgttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATcc                                                                                                                                     >  1:192188/1‑148 (MQ=255)
          ttgttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATcc                                                                                                                                     >  2:1238399/1‑148 (MQ=255)
            gttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAg                                                                                                                                   <  1:3444384/148‑1 (MQ=255)
            gttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAg                                                                                                                                   <  2:1062695/148‑1 (MQ=255)
             ttgATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGa                                                                                                                                  >  1:532662/1‑148 (MQ=255)
                 ttATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGa                                                                                                                              >  2:2527989/1‑148 (MQ=255)
                    tCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACg                                                                                                                           <  1:1822786/148‑1 (MQ=255)
                     cGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGt                                                                                                                          <  1:1811783/148‑1 (MQ=255)
                     cGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGt                                                                                                                          <  1:1757699/148‑1 (MQ=255)
                     cGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACACCGCCTTAATCCTATCCACCACTTGCAAGGAGCTGTAGTAATCCAGACGGAACGt                                                                                                                          >  1:1084011/1‑148 (MQ=255)
                     cGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGAGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGt                                                                                                                          <  1:2903862/148‑1 (MQ=255)
                     cGGCGTGGTGTTATTTCTGGTGCTGGTGGAGTTGCGACAATTTCGCCAGACGCCGCCGCACGTGACAGCGTCCGACACTTAATGCTTAAAAAAGCGCCATAAGCCAATACACCACATGCATGGAGCTGTAGTAATCCAGACGGAACGt                                                                                                                          >  1:1083927/1‑148 (MQ=255)
                       gCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGtct                                                                                                                        <  2:797059/148‑1 (MQ=255)
                         gtgtTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCg                                                                                                                      <  1:1313270/148‑1 (MQ=255)
                         gtgtTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCg                                                                                                                      <  1:3058781/148‑1 (MQ=255)
                          tgttgtTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCgg                                                                                                                     <  1:2169233/148‑1 (MQ=255)
                           gttgttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGt                                                                                                                    >  1:3258900/1‑148 (MQ=255)
                            ttgttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGtt                                                                                                                   >  2:2639795/1‑148 (MQ=255)
                              gttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTcc                                                                                                                 <  1:1237156/148‑1 (MQ=255)
                              gttATTGCTGGTGCTGGTGGAGGTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTcc                                                                                                                 <  2:851896/148‑1 (MQ=255)
                               ttATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTccc                                                                                                                <  1:3093684/148‑1 (MQ=255)
                                tATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCa                                                                                                               >  1:26530/1‑148 (MQ=255)
                                tATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCa                                                                                                               >  1:1205475/1‑148 (MQ=255)
                                       ggtgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATa                                                                                                        >  1:1363067/1‑148 (MQ=255)
                                        gtgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATaa                                                                                                       >  1:46589/1‑148 (MQ=255)
                                        gtgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATaa                                                                                                       <  2:263330/148‑1 (MQ=255)
                                         tgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATaaa                                                                                                      >  1:768017/1‑148 (MQ=255)
                                         tgctggtgGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATaaa                                                                                                      >  1:8679/1‑148 (MQ=255)
                                                ggAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCtt                                                                                               >  1:1150935/1‑148 (MQ=255)
                                                      gCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTg                                                                                         <  1:2639795/148‑1 (MQ=255)
                                                        gACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTa                                                                                       <  2:3116887/148‑1 (MQ=255)
                                                             tttCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTgca                                                                                  <  2:2343290/148‑1 (MQ=255)
                                                             tttCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTgca                                                                                  >  2:1299293/1‑148 (MQ=255)
                                                             tttCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTgca                                                                                  >  2:1299166/1‑148 (MQ=255)
                                                                    aGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGt                                                                           <  2:2101522/148‑1 (MQ=255)
                                                                    aGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGt                                                                           <  2:2101498/148‑1 (MQ=255)
                                                                     gACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTg                                                                          <  1:992144/148‑1 (MQ=255)
                                                                      aCGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTgc                                                                         >  2:2449774/1‑148 (MQ=255)
                                                                      aCGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTgc                                                                         >  2:2449961/1‑148 (MQ=255)
                                                                      aCGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTgc                                                                         >  1:3478807/1‑148 (MQ=255)
                                                                           gccgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGag                                                                    <  1:1890664/148‑1 (MQ=255)
                                                                            ccgcAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGagc                                                                   >  1:1375978/1‑148 (MQ=255)
                                                                                 ggTGACAGCTTAAGAAATTTAATGCTTAAAACAGCGCCTTAAGCCTTTCCAGCACTTGCTTGGTGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATATACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCgg                                                              <  2:2633592/148‑1 (MQ=255)
                                                                                  gTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGa                                                             <  2:1682159/148‑1 (MQ=255)
                                                                                  gTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGa                                                             <  1:1774502/148‑1 (MQ=255)
                                                                                    gACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGAtt                                                           <  2:3271400/148‑1 (MQ=255)
                                                                                     aCAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTa                                                          <  1:1689399/148‑1 (MQ=255)
                                                                                     aCAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTa                                                          <  2:3096979/148‑1 (MQ=255)
                                                                                     aCAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTa                                                          <  1:46949/148‑1 (MQ=255)
                                                                                       aGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGc                                                        <  2:3306024/148‑1 (MQ=255)
                                                                                        gCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCa                                                       <  2:1069928/148‑1 (MQ=255)
                                                                                        gCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCa                                                       <  1:2805652/148‑1 (MQ=255)
                                                                                           tCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGANCGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATaa                                                    <  1:1783842/148‑1 (MQ=255)
                                                                                           tCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATaa                                                    <  1:1156878/148‑1 (MQ=255)
                                                                                           tCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATaa                                                    <  1:2996305/148‑1 (MQ=255)
                                                                                            ccGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATaaa                                                   <  2:532662/148‑1 (MQ=255)
                                                                                              gACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATa                                                 >  2:1521046/1‑148 (MQ=255)
                                                                                                cAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGc                                               <  2:509099/148‑1 (MQ=255)
                                                                                                cAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGc                                               >  1:1809862/1‑148 (MQ=255)
                                                                                                cAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGc                                               <  1:1917728/148‑1 (MQ=255)
                                                                                                 aGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCa                                              <  2:1493025/148‑1 (MQ=255)
                                                                                                  gTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAt                                             >  2:1885242/1‑148 (MQ=255)
                                                                                                  gTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAt                                             <  2:209117/148‑1 (MQ=255)
                                                                                                  gTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAt                                             <  1:561126/148‑1 (MQ=255)
                                                                                                  gTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCAt                                             >  2:765770/1‑148 (MQ=255)
                                                                                                   ttAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATc                                            <  2:239954/148‑1 (MQ=255)
                                                                                                       tGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCa                                        >  2:1872618/1‑148 (MQ=255)
                                                                                                       tGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCa                                        >  2:1808191/1‑148 (MQ=255)
                                                                                                        gCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCAt                                       >  2:1833636/1‑148 (MQ=255)
                                                                                                        gCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCAt                                       >  2:2829240/1‑148 (MQ=255)
                                                                                                         cTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATc                                      <  2:657427/148‑1 (MQ=255)
                                                                                                          ttAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCt                                     >  1:3042597/1‑148 (MQ=255)
                                                                                                          ttAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCt                                     >  2:551547/1‑148 (MQ=255)
                                                                                                          ttAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCt                                     >  1:1295924/1‑148 (MQ=255)
                                                                                                            aaaaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCttt                                   <  2:2825157/148‑1 (MQ=255)
                                                                                                            aaaaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCttt                                   <  2:1072145/148‑1 (MQ=255)
                                                                                                              aaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCt                                 >  1:2486655/1‑148 (MQ=255)
                                                                                                              aaCAGCGCCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCt                                 >  1:2012325/1‑148 (MQ=255)
                                                                                                                  gcgcCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATc                             >  1:349908/1‑148 (MQ=255)
                                                                                                                   cgcCTTAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCa                            >  2:2960039/1‑148 (MQ=255)
                                                                                                                       ttAAGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCgg                        <  1:1238399/148‑1 (MQ=255)
                                                                                                                         aaGCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCg                      <  2:2360199/148‑1 (MQ=255)
                                                                                                                           gCCTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGaa                    <  1:1126565/148‑1 (MQ=255)
                                                                                                                            ccTATCCAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAAc                   <  2:2493643/148‑1 (MQ=255)
                                                                                                                                 ccAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGaa              <  1:2896702/148‑1 (MQ=255)
                                                                                                                                  cAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAAt             <  1:2154912/148‑1 (MQ=255)
                                                                                                                                  cAGCACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAAt             >  1:753216/1‑148 (MQ=255)
                                                                                                                                     cACTTGCATGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGt          >  1:2858483/1‑148 (MQ=255)
                                                                                                                                            aTGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGc   <  2:1655099/148‑1 (MQ=255)
                                                                                                                                             tGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGcc  >  2:3424876/1‑148 (MQ=255)
                                                                                                                                             tGGAGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGcc  >  2:1655752/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCATGGCGCTGTAGTAATCCAGACGGAACGTCTCGGTTCCCAGCGCATAAACCTGCTTGTTTTGTACTGCAGGCAGGTGCGCGAGCAGCGGATTAGCATAAATAGCATCGGCATCTTTCTGATCACCGGCGAACAGGAATAGTGACTCGCC  >  CP009273/618904‑619192

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: