Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,077,911 A→C 92.6% 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,9110AC92.6% 713.3 / 23.1 207K71N (AAA→AACefeOinactive ferrous ion transporter EfeUOB
Reads supporting (aligned to +/- strand):  ref base A (6/11);  new base C (114/76);  total (120/87)
Fisher's exact test for biased strand distribution p-value = 7.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

ATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAA  >  CP009273/1077765‑1078056
                                                                                                                                                  |                                                                                                                                                 
aTGGCTAACGCCGCTGATGTGCCGCAGTTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACTCCGGGAAAACACATTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGc                                                                                                                                                <  1:1844653/150‑1 (MQ=255)
aTGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGc                                                                                                                                                <  1:3249033/150‑1 (MQ=255)
 tGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCg                                                                                                                                               >  1:410687/1‑150 (MQ=255)
 tGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCg                                                                                                                                               >  1:1607161/1‑150 (MQ=255)
 tGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCg                                                                                                                                               >  2:1832804/1‑150 (MQ=255)
     tAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGAt                                                                                                                                           >  1:3412668/1‑150 (MQ=255)
         gccgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggtg                                                                                                                                       >  1:1664099/1‑150 (MQ=255)
         gccgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggtg                                                                                                                                       >  1:1664256/1‑150 (MQ=255)
         gccgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGAtggtg                                                                                                                                       >  2:4639286/1‑150 (MQ=255)
          ccgcTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggtgg                                                                                                                                      >  2:3529284/1‑150 (MQ=255)
            acTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAtggtggtg                                                                                                                                    >  1:1223238/2‑150 (MQ=255)
              tGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGa                                                                                                                                  <  1:4613700/150‑1 (MQ=255)
              tGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGa                                                                                                                                  <  1:4614531/150‑1 (MQ=255)
                 tgtgCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGAGGGTGGTGGAaga                                                                                                                               >  2:3210188/1‑150 (MQ=255)
                  gtgCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAagag                                                                                                                              >  2:4996203/1‑150 (MQ=255)
                   tgCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGc                                                                                                                             <  1:118877/150‑1 (MQ=255)
                    gCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCg                                                                                                                            <  1:3790560/150‑1 (MQ=255)
                    gCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCg                                                                                                                            <  2:2047892/150‑1 (MQ=255)
                       gCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGa                                                                                                                         <  2:4423560/150‑1 (MQ=255)
                        cAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGaa                                                                                                                        <  2:3617406/150‑1 (MQ=255)
                          ggTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGaaaa                                                                                                                      >  1:1945572/1‑150 (MQ=255)
                              aaaGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATc                                                                                                                  <  1:414486/150‑1 (MQ=255)
                              aaaGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATc                                                                                                                  >  2:1376939/1‑150 (MQ=255)
                                   gACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGcccc                                                                                                             <  2:1365180/150‑1 (MQ=255)
                                   gACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGcccc                                                                                                             >  1:2244339/1‑150 (MQ=255)
                                    aCCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCt                                                                                                            <  2:2560995/150‑1 (MQ=255)
                                    aCCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCt                                                                                                            <  1:407334/150‑1 (MQ=255)
                                    aCCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCt                                                                                                            <  2:3097999/150‑1 (MQ=255)
                                     cnGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTg                                                                                                           >  1:2766620/1‑150 (MQ=255)
                                      cGTGACGGATAAGCAGTGCGAACCGATGACAATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTgg                                                                                                          <  1:2301460/150‑1 (MQ=255)
                                        tGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCt                                                                                                        <  2:3240267/150‑1 (MQ=255)
                                          aCGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCttt                                                                                                      <  1:52899/150‑1 (MQ=255)
                                          aCGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCttt                                                                                                      <  2:2273324/150‑1 (MQ=255)
                                           cGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTa                                                                                                     <  2:2333332/150‑1 (MQ=255)
                                            ggATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAg                                                                                                    >  2:1406078/1‑150 (MQ=255)
                                              aTAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGcc                                                                                                  <  1:125837/150‑1 (MQ=255)
                                                aaGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAg                                                                                                >  2:2175800/1‑150 (MQ=255)
                                                 aGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGa                                                                                               <  2:4143646/150‑1 (MQ=255)
                                                   cAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATCTCGCCCCTGGCTTTAGCCAGaaa                                                                                             >  2:105167/1‑150 (MQ=255)
                                                     gTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAAt                                                                                           >  1:2646902/1‑150 (MQ=255)
                                                        cGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGAc                                                                                        >  1:4678751/1‑150 (MQ=255)
                                                        cGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGAc                                                                                        >  1:3162339/1‑150 (MQ=255)
                                                         gAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACg                                                                                       >  2:3185685/1‑150 (MQ=255)
                                                           aCCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGc                                                                                     <  1:2447823/150‑1 (MQ=255)
                                                             ggATGACCATTCCGGTAAACGCCGGGAAAACACATTTCCTTTTTCATAACCACCGCCTGAAGGCGCTGGAGTGGTAGATCCTAAACGGCGTGATGGTGGTGGAAGTGGGGGAAAATATCGCACCTGTCTTTAGCCAGAAAATGACGGcta                                                                                   <  2:1398223/149‑3 (MQ=255)
                                                               aTGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAAt                                                                                 <  1:4504548/150‑1 (MQ=255)
                                                               aTGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAAt                                                                                 >  1:772563/1‑150 (MQ=255)
                                                                tGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAAtt                                                                                >  2:63359/1‑150 (MQ=255)
                                                                 gACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAAttt                                                                               <  2:4860817/150‑1 (MQ=255)
                                                                      ttACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGc                                                                          >  1:750546/1‑150 (MQ=255)
                                                                       tACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGcc                                                                         >  2:1679968/1‑150 (MQ=255)
                                                                        aCGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCt                                                                        >  1:2197291/1‑150 (MQ=255)
                                                                        aCGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCt                                                                        >  1:1761566/1‑150 (MQ=255)
                                                                        aCGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCt                                                                        >  2:1582044/1‑150 (MQ=255)
                                                                          ggTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTgg                                                                      <  1:3534273/150‑1 (MQ=255)
                                                                             tAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGa                                                                   >  2:340207/1‑150 (MQ=255)
                                                                               aCGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAAt                                                                 <  2:5030649/150‑1 (MQ=255)
                                                                               aCGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAAt                                                                 <  2:2747240/150‑1 (MQ=255)
                                                                                cGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATa                                                                <  1:600023/150‑1 (MQ=255)
                                                                                  ccGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACg                                                              >  2:3303462/1‑150 (MQ=255)
                                                                                  ccGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACg                                                              >  2:4916066/1‑150 (MQ=255)
                                                                                    gggAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGag                                                            >  1:3291221/1‑149 (MQ=255)
                                                                                      gAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGatat                                                          >  2:4480885/1‑150 (MQ=255)
                                                                                       aaaaCACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGACTACGATATg                                                         >  1:2843551/1‑150 (MQ=255)
                                                                                        aaaCACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGa                                                        <  2:674548/150‑1 (MQ=255)
                                                                                          acacAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGAcc                                                      <  1:869474/150‑1 (MQ=255)
                                                                                             caGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGc                                                   <  2:3890184/150‑1 (MQ=255)
                                                                                                ttCATTATTCAGAACCACAGCCAGAAGGCTCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGt                                                <  2:3028394/150‑1 (MQ=255)
                                                                                                  cATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTct                                              <  2:2431786/150‑1 (MQ=255)
                                                                                                  cATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTct                                              >  2:435919/1‑150 (MQ=255)
                                                                                                  cATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTct                                              <  1:2143779/150‑1 (MQ=255)
                                                                                                   attattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctg                                             <  2:3827745/150‑1 (MQ=255)
                                                                                                     tattCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTctgct                                           <  1:2160329/150‑1 (MQ=255)
                                                                                                        tCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGAc                                        <  2:5130082/150‑1 (MQ=255)
                                                                                                         cAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACt                                       >  1:3238382/1‑150 (MQ=255)
                                                                                                         cAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACt                                       >  2:2833102/1‑150 (MQ=255)
                                                                                                         cAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACt                                       >  1:4965192/1‑150 (MQ=255)
                                                                                                           gAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTaa                                     >  2:2702003/1‑150 (MQ=255)
                                                                                                            aaCCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAAc                                    <  2:1693813/150‑1 (MQ=255)
                                                                                                             aCCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAAcc                                   >  2:2377587/1‑150 (MQ=255)
                                                                                                             aCCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAAcc                                   <  1:1999314/150‑1 (MQ=255)
                                                                                                              caacaGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGTCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCAAGAAAATGACGGCGAATTTACAGACTGGCGAATACGATATGACCTGCTGTCTGCTGACTAAccc                                  >  1:3959630/3‑150 (MQ=255)
                                                                                                               cacaGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCg                                 >  2:419917/1‑150 (MQ=255)
                                                                                                                acaGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGa                                >  2:4108866/1‑150 (MQ=255)
                                                                                                                  aGCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGaaa                              >  1:3983005/1‑150 (MQ=255)
                                                                                                                   gCCAGAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAg                             <  2:1797546/150‑1 (MQ=255)
                                                                                                                      aGAAGGCGCTGGAGTGGGAGGTCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGa                          >  1:1499743/1‑150 (MQ=255)
                                                                                                                       gAAGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGaa                         >  1:1628426/1‑150 (MQ=255)
                                                                                                                        aaGGCGCTGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACATGCGGTCTGCTGACTAACCCGAAAGGGAAg                        >  2:3133870/1‑150 (MQ=255)
                                                                                                                               tGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCg                 <  1:1448042/150‑1 (MQ=255)
                                                                                                                               tGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCg                 >  1:1021141/1‑150 (MQ=255)
                                                                                                                               tGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCg                 >  2:3713096/1‑150 (MQ=255)
                                                                                                                               tGGAGTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCg                 >  1:4239318/1‑150 (MQ=255)
                                                                                                                                   gTGGGAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCaa             >  2:1220303/1‑150 (MQ=255)
                                                                                                                                     gggAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAg           >  2:3967024/1‑150 (MQ=255)
                                                                                                                                      ggAGATCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAgg          <  1:1063868/150‑1 (MQ=255)
                                                                                                                                       gagaTCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGt         <  2:2090228/150‑1 (MQ=255)
                                                                                                                                         gaTCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGa       >  2:2437645/1‑150 (MQ=255)
                                                                                                                                         gaTCCTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGa       <  1:2056245/150‑1 (MQ=255)
                                                                                                                                            ccTCAACGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGc    >  2:1253211/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATGGCTAACGCCGCTGATGTGCCGCAGGTCAAAGTGACCGTGACGGATAAGCAGTGCGAACCGATGACCATTACGGTTAACGCCGGGAAAACACAGTTCATTATTCAGAACCACAGCCAGAAGGCGCTGGAGTGGGAGATCCTCAAAGGCGTGATGGTGGTGGAAGAGCGGGAAAATATCGCCCCTGGCTTTAGCCAGAAAATGACGGCGAATTTACAGCCTGGCGAATACGATATGACCTGCGGTCTGCTGACTAACCCGAAAGGGAAGTTGATCGTCAAAGGTGAGGCAA  >  CP009273/1077765‑1078056

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: