Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,077,911 A→C 100% K71N (AAA→AAC efeO → inactive ferrous ion transporter EfeUOB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,077,9110AC99.6% 848.0 / ‑5.4 210K71N (AAA→AACefeOinactive ferrous ion transporter EfeUOB
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (103/106);  total (104/106)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 210 total aligned reads displayed.

TGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCA  >  CP009273/1077766‑1078055
                                                                                                                                                 |                                                                                                                                                
tGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACgg                                                                                                                                                >  2:818683/1‑148 (MQ=255)
tGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACgg                                                                                                                                                <  2:2446280/148‑1 (MQ=255)
  gCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCg                                                                                                                                              <  2:1813017/148‑1 (MQ=255)
    tAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTg                                                                                                                                            <  2:2510603/148‑1 (MQ=255)
    tAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTg                                                                                                                                            <  1:853924/148‑1 (MQ=255)
      aCGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAt                                                                                                                                          >  1:2914977/1‑148 (MQ=255)
       cgccgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtg                                                                                                                                         <  1:1420246/148‑1 (MQ=255)
       cgccgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtg                                                                                                                                         <  1:2415059/148‑1 (MQ=255)
         ccgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAgggt                                                                                                                                       >  1:1073951/1‑148 (MQ=255)
            cTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggtggt                                                                                                                                    >  2:989027/1‑148 (MQ=255)
                  tgCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAaga                                                                                                                              >  2:1782874/1‑148 (MQ=255)
                  tgCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAaga                                                                                                                              <  2:1228428/148‑1 (MQ=255)
                  tgCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAaga                                                                                                                              >  2:3543553/1‑148 (MQ=255)
                   gCCGCAGGTCAACGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAagag                                                                                                                             <  1:1301390/148‑1 (MQ=255)
                   gCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGTGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAagag                                                                                                                             <  1:2220672/148‑1 (MQ=255)
                       cAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCggg                                                                                                                         >  2:783265/1‑148 (MQ=255)
                         ggTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGaa                                                                                                                       >  1:2444457/1‑148 (MQ=255)
                         ggTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGaa                                                                                                                       >  1:3551563/1‑148 (MQ=255)
                          gTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGaaa                                                                                                                      <  2:47288/148‑1 (MQ=255)
                            cAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAat                                                                                                                    <  1:41049/148‑1 (MQ=255)
                            cAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAat                                                                                                                    <  2:613757/148‑1 (MQ=255)
                             aaaGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAata                                                                                                                   <  1:1505616/148‑1 (MQ=255)
                             aaaGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAata                                                                                                                   <  1:3319348/148‑1 (MQ=255)
                                gTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCg                                                                                                                <  1:1818607/148‑1 (MQ=255)
                                    ccGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGcccc                                                                                                            <  2:1946376/148‑1 (MQ=255)
                                      gTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTg                                                                                                          <  2:2954482/148‑1 (MQ=255)
                                       tGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTgg                                                                                                         >  1:2705064/1‑148 (MQ=255)
                                       tGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTgg                                                                                                         <  1:806013/148‑1 (MQ=255)
                                          cGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCtt                                                                                                      <  2:175608/148‑1 (MQ=255)
                                             aTAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAg                                                                                                   >  2:318039/1‑148 (MQ=255)
                                              tAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGc                                                                                                  <  1:3514323/148‑1 (MQ=255)
                                              tAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGc                                                                                                  >  2:910951/1‑148 (MQ=255)
                                               aaGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGcc                                                                                                 >  1:1172140/1‑148 (MQ=255)
                                                aGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCa                                                                                                >  1:484674/1‑148 (MQ=255)
                                                aGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCa                                                                                                >  1:914664/1‑148 (MQ=255)
                                                 gCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAg                                                                                               >  2:3246846/1‑148 (MQ=255)
                                                  cAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGa                                                                                              <  1:172048/148‑1 (MQ=255)
                                                    gTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGaaa                                                                                            <  2:1846051/148‑1 (MQ=255)
                                                    gTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGaaa                                                                                            <  2:2620294/148‑1 (MQ=255)
                                                    gTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGaaa                                                                                            >  1:203391/1‑148 (MQ=255)
                                                     tGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGaaaa                                                                                           >  1:497887/1‑148 (MQ=255)
                                                         aaCCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGAc                                                                                       >  1:1067516/1‑148 (MQ=255)
                                                          cccGATTACATTTACGTTTAACGCCGTTAAAACACATTTAATTATTCACACCCACCGCCAGAAGGGGGTGGAGTGGTATATCCTAAACGGGGTTTTGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACg                                                                                      <  1:1230055/147‑1 (MQ=255)
                                                             gATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCg                                                                                   <  2:1828294/148‑1 (MQ=255)
                                                              aTGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGa                                                                                  >  2:1223574/1‑148 (MQ=255)
                                                               tGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGaa                                                                                 <  2:2126564/148‑1 (MQ=255)
                                                                  ccATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAAttt                                                                              <  2:2415399/148‑1 (MQ=255)
                                                                     ttACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACa                                                                           >  1:60985/1‑148 (MQ=255)
                                                                      tACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAg                                                                          <  2:52232/148‑1 (MQ=255)
                                                                      tACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAg                                                                          >  2:1131159/1‑148 (MQ=255)
                                                                        cGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGcc                                                                        >  2:1109647/1‑148 (MQ=255)
                                                                          gTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTg                                                                      >  2:689974/1‑148 (MQ=255)
                                                                               cGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGaa                                                                 <  1:891598/148‑1 (MQ=255)
                                                                                gCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAAt                                                                <  2:101550/148‑1 (MQ=255)
                                                                                 ccGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATa                                                               <  2:2897974/148‑1 (MQ=255)
                                                                                   gggAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACg                                                             >  1:1980444/1‑148 (MQ=255)
                                                                                   gggAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACg                                                             <  2:284527/148‑1 (MQ=255)
                                                                                    ggAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGa                                                            <  2:2222622/148‑1 (MQ=255)
                                                                                    ggAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGa                                                            <  2:1273749/148‑1 (MQ=255)
                                                                                     gAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGat                                                           >  1:3017600/1‑148 (MQ=255)
                                                                                      aaaaCACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGata                                                          >  2:2148435/1‑148 (MQ=255)
                                                                                        aaCACAGTTCATTATTCAGAACCACAGCCATAAGGCGCTGGAGTGTGAGATCCTCAACGGCGTGATGGTGGTGTAAGAGCGGGAAAATATCTCCCCTTGCTTTAGCCAGAAAATGACGGCGAATTTACATCCTGGCTAATACGATATg                                                        >  1:2051672/1‑148 (MQ=255)
                                                                                        aaCACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATg                                                        >  2:1500854/1‑148 (MQ=255)
                                                                                        aaCACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATg                                                        >  2:2204927/1‑148 (MQ=255)
                                                                                             aGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTg                                                   <  1:2811388/148‑1 (MQ=255)
                                                                                               ttCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCg                                                 >  1:2448228/1‑148 (MQ=255)
                                                                                               ttCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCg                                                 >  1:3252996/1‑148 (MQ=255)
                                                                                                 cATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGt                                               <  1:2117507/148‑1 (MQ=255)
                                                                                                 cATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGt                                               <  1:687811/148‑1 (MQ=255)
                                                                                                   ttattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGTCGAATATACAGCCTGGCGAATACGATATGACCTGCGGTct                                             >  1:3042845/1‑148 (MQ=255)
                                                                                                   ttattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTct                                             >  1:2599791/1‑148 (MQ=255)
                                                                                                    tattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCACCGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctg                                            <  1:759945/148‑1 (MQ=255)
                                                                                                    tattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctg                                            >  1:2160492/1‑148 (MQ=255)
                                                                                                    tattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctg                                            <  2:3357972/148‑1 (MQ=255)
                                                                                                      ttCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctgct                                          <  1:1782874/148‑1 (MQ=255)
                                                                                                      ttCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctgct                                          <  2:1985203/148‑1 (MQ=255)
                                                                                                        cAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGa                                        <  1:3402541/148‑1 (MQ=255)
                                                                                                        cAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGAGAATACGATATGACCTGCGGTCTGCTGa                                        >  2:2083943/1‑148 (MQ=255)
                                                                                                         aGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGAc                                       <  2:485695/148‑1 (MQ=255)
                                                                                                          gAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACt                                      >  2:3027562/1‑148 (MQ=255)
                                                                                                            aCCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTaa                                    <  2:302520/148‑1 (MQ=255)
                                                                                                             ccACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAAc                                   <  1:2496338/148‑1 (MQ=255)
                                                                                                               acaGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAAccc                                 >  1:1651356/1‑148 (MQ=255)
                                                                                                                  gCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGaa                              <  1:1893918/148‑1 (MQ=255)
                                                                                                                    cAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAg                            <  2:3551563/148‑1 (MQ=255)
                                                                                                                      gAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATTACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAggg                          >  1:1600009/1‑148 (MQ=255)
                                                                                                                      gAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAggg                          >  1:2210471/1‑148 (MQ=255)
                                                                                                                             cTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGACGTTGa                   >  1:2274770/1‑148 (MQ=255)
                                                                                                                             cTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGa                   <  2:2590028/148‑1 (MQ=255)
                                                                                                                             cTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGa                   <  1:61719/148‑1 (MQ=255)
                                                                                                                                 aGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGt               <  1:386496/148‑1 (MQ=255)
                                                                                                                                     ggAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCaaa           <  2:2548472/148‑1 (MQ=255)
                                                                                                                                     ggAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCaaa           >  2:941091/1‑148 (MQ=255)
                                                                                                                                     ggAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCaaa           >  2:2723542/1‑148 (MQ=255)
                                                                                                                                       agaTCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAgg         >  2:3463079/1‑148 (MQ=255)
                                                                                                                                       agaTCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAgg         <  2:252579/148‑1 (MQ=255)
                                                                                                                                         aTCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTg       >  1:973363/1‑148 (MQ=255)
                                                                                                                                          tCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGa      >  1:1718821/1‑148 (MQ=255)
                                                                                                                                              cAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCa  >  1:156174/1‑148 (MQ=255)
                                                                                                                                              cAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCa  >  1:1681242/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCA  >  CP009273/1077766‑1078055

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: