Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 319,381 G→T 100% D445Y (GAC→TAC)  ykgF → ferridoxin‑like LutB family protein; putative electron transport chain YkgEFG component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273319,3810GT99.5% 661.2 / ‑4.2 199D445Y (GAC→TAC) ykgFferridoxin‑like LutB family protein; putative electron transport chain YkgEFG component
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (72/126);  total (72/127)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
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 199 total aligned reads displayed.

AAAAGCAGAGCAACGGGCGATAAAAATGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCGACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGTTGCTCAGGCACTGGGTCGCCCGC  >  CP009273/319236‑319523
                                                                                                                                                 |                                                                                                                                              
aaaaGCAGAGCAACGGGCGATAAAAATGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTAc                                                                                                                                              >  2:682661/1‑148 (MQ=255)
aaaaGCAGAGCAACGGGCGATAAAAATGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTAc                                                                                                                                              <  1:3690045/148‑1 (MQ=255)
   aGCAGAGCAACGGGCGATAAAAATGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACtgg                                                                                                                                           >  1:743922/1‑148 (MQ=255)
    gCAGAGCAACGGGCGATAAAAATGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACtgga                                                                                                                                          <  1:301388/148‑1 (MQ=255)
           aaCGGGCGATAAAAATGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGc                                                                                                                                   <  2:839066/148‑1 (MQ=255)
            aCGGGCGATAAAAATGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCa                                                                                                                                  >  2:149564/1‑148 (MQ=255)
            aCGGGCGATAAAAATGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCa                                                                                                                                  <  2:2242908/148‑1 (MQ=255)
              gggCGATAAAAATGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCAcg                                                                                                                                <  1:1207929/148‑1 (MQ=255)
                      aaaaTGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTc                                                                                                                        <  1:1100638/148‑1 (MQ=255)
                      aaaaTGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTc                                                                                                                        >  2:2187552/1‑148 (MQ=255)
                       aaaTGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTcc                                                                                                                       >  2:2545927/1‑148 (MQ=255)
                        aaTGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCt                                                                                                                      >  1:934651/1‑148 (MQ=255)
                          tGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGa                                                                                                                    <  1:3322562/148‑1 (MQ=255)
                          tGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGa                                                                                                                    >  2:1930736/1‑148 (MQ=255)
                          tGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGa                                                                                                                    <  1:3396126/148‑1 (MQ=255)
                          tGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGa                                                                                                                    <  2:3337564/148‑1 (MQ=255)
                                cTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGa                                                                                                              <  2:2128469/148‑1 (MQ=255)
                                 ttATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTCCTGGATGGAAGCACGCGATCTTCCTGAAGCTGAc                                                                                                             >  1:1337900/1‑148 (MQ=255)
                                 ttATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGAc                                                                                                             >  2:425774/1‑148 (MQ=255)
                                   aTGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACgg                                                                                                           >  1:3103118/1‑148 (MQ=255)
                                      ccAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGACCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGgaga                                                                                                        <  2:3394730/148‑1 (MQ=255)
                                       cAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGgagag                                                                                                       >  1:1461805/1‑148 (MQ=255)
                                        aaTAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGgagaga                                                                                                      >  1:3659531/1‑148 (MQ=255)
                                         aTAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGagagag                                                                                                     >  2:880367/1‑148 (MQ=255)
                                         aTAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGagagag                                                                                                     <  1:484324/148‑1 (MQ=255)
                                         aTAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGagagag                                                                                                     <  2:2733853/148‑1 (MQ=255)
                                          tAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGt                                                                                                    <  1:1589937/148‑1 (MQ=255)
                                              cATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTcc                                                                                                <  2:209740/148‑1 (MQ=255)
                                               aTCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCg                                                                                               >  2:3089987/1‑148 (MQ=255)
                                                 ccAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTa                                                                                             >  2:3456336/1‑148 (MQ=255)
                                                     gATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTg                                                                                         <  2:2445931/148‑1 (MQ=255)
                                                     gATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTg                                                                                         <  2:525208/148‑1 (MQ=255)
                                                        tgtgGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGtt                                                                                      <  2:1131949/148‑1 (MQ=255)
                                                          tgGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTa                                                                                    >  2:3171047/1‑148 (MQ=255)
                                                          tgGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTa                                                                                    >  2:2526331/1‑148 (MQ=255)
                                                             aaaGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGa                                                                                 <  1:2856478/148‑1 (MQ=255)
                                                             aaaGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGa                                                                                 >  2:3659169/1‑148 (MQ=255)
                                                               aGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGaaa                                                                               >  1:2174362/1‑148 (MQ=255)
                                                               aGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGaaa                                                                               <  2:2710325/148‑1 (MQ=255)
                                                                  cGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACAt                                                                            <  1:252034/148‑1 (MQ=255)
                                                                   gggATGATGGCCGGTGCTAATGCGGCAAGCTGGTTTATCAATGTCGGCAAAACAACACTCAAATTTGGCGCGATTAGCTACTGGATTGAAGCACGAGATCTTCCTGAAGCTTACGGAGAGAGTTTCCGTAGTTGGTTTAATAAACATc                                                                           <  2:571978/148‑1 (MQ=255)
                                                                    ggATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCa                                                                          <  1:533056/148‑1 (MQ=255)
                                                                    ggATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCa                                                                          <  1:2301340/148‑1 (MQ=255)
                                                                      atgatgGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAgg                                                                        <  1:475855/148‑1 (MQ=255)
                                                                      atgatgGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAgg                                                                        <  1:360339/148‑1 (MQ=255)
                                                                         atgGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGgcgc                                                                     <  2:1086369/148‑1 (MQ=255)
                                                                         atgGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGgcgc                                                                     <  2:1086402/148‑1 (MQ=255)
                                                                          tgGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCa                                                                    <  1:2145051/148‑1 (MQ=255)
                                                                           ggCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAg                                                                   <  1:2345420/148‑1 (MQ=255)
                                                                            gCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAgg                                                                  <  2:2056931/148‑1 (MQ=255)
                                                                             ccGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGga                                                                 <  1:2541796/148‑1 (MQ=255)
                                                                                  gCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGaaaa                                                            <  1:875282/148‑1 (MQ=255)
                                                                                  gCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGaaaa                                                            <  2:853548/148‑1 (MQ=255)
                                                                                    tCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAg                                                          <  1:2652038/148‑1 (MQ=255)
                                                                                     cATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGa                                                         <  2:743922/148‑1 (MQ=255)
                                                                                     cATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGa                                                         <  2:934681/148‑1 (MQ=255)
                                                                                     cATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGa                                                         <  1:1930736/148‑1 (MQ=255)
                                                                                     cATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGa                                                         <  1:79469/148‑1 (MQ=255)
                                                                                       tGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAAt                                                       <  2:1349299/148‑1 (MQ=255)
                                                                                       tGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAAt                                                       <  2:1023930/148‑1 (MQ=255)
                                                                                        gCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATg                                                      <  2:2038098/148‑1 (MQ=255)
                                                                                         cGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATgg                                                     <  1:834637/148‑1 (MQ=255)
                                                                                          ggCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGa                                                    <  1:3174922/148‑1 (MQ=255)
                                                                                             aaGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATaa                                                 <  1:254069/148‑1 (MQ=255)
                                                                                             aaGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATaa                                                 <  2:1462176/148‑1 (MQ=255)
                                                                                                cTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCg                                              <  1:363812/148‑1 (MQ=255)
                                                                                                cTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCg                                              <  1:1003523/148‑1 (MQ=255)
                                                                                                cTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCg                                              <  1:2503665/148‑1 (MQ=255)
                                                                                                  ggTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAg                                            >  2:854347/1‑148 (MQ=255)
                                                                                                  ggTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAg                                            <  2:3586367/148‑1 (MQ=255)
                                                                                                      tATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGa                                        <  1:2711629/148‑1 (MQ=255)
                                                                                                      tATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGa                                        >  1:3675884/1‑148 (MQ=255)
                                                                                                       aTCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGaa                                       <  2:1784123/148‑1 (MQ=255)
                                                                                                       aTCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGaa                                       <  2:3674715/148‑1 (MQ=255)
                                                                                                       aTCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGaa                                       <  2:3659531/148‑1 (MQ=255)
                                                                                                       aTCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGaa                                       >  2:1030631/1‑148 (MQ=255)
                                                                                                         cAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAAtt                                     >  1:2811742/1‑148 (MQ=255)
                                                                                                          aaTGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAAttt                                    <  1:677672/148‑1 (MQ=255)
                                                                                                           aTGGCGGCAAAACACCACTCAAATTTGGTGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAAtttt                                   >  1:1724930/1‑148 (MQ=255)
                                                                                                           aTGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAAtttt                                   >  1:1794309/1‑148 (MQ=255)
                                                                                                           aTGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAAtttt                                   <  1:2882171/148‑1 (MQ=255)
                                                                                                               cggcAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGaa                               <  2:1842973/148‑1 (MQ=255)
                                                                                                               cggcAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGaa                               <  2:3732453/148‑1 (MQ=255)
                                                                                                               cggcAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGaa                               <  2:2601466/148‑1 (MQ=255)
                                                                                                                 gcAAAACACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATa                             <  2:3143532/148‑1 (MQ=255)
                                                                                                                   aaaaCACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAAc                           <  2:1356488/148‑1 (MQ=255)
                                                                                                                    aaaCACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACg                          <  2:2357577/148‑1 (MQ=255)
                                                                                                                     cgcaccacTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGt                         <  1:2016850/146‑1 (MQ=255)
                                                                                                                     aaCACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGt                         <  2:3015851/148‑1 (MQ=255)
                                                                                                                     aaCACCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGt                         >  1:564676/1‑148 (MQ=255)
                                                                                                                      acacCACTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGtt                        >  1:1054693/1‑148 (MQ=255)
                                                                                                                        accacTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGTTGc                      >  2:1688280/1‑148 (MQ=255)
                                                                                                                         ccacTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGTTGCt                     <  2:589632/148‑1 (MQ=255)
                                                                                                                         ccacTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGTTGCt                     <  1:1527996/148‑1 (MQ=255)
                                                                                                                          cacTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGTTGCTc                    <  2:2582689/148‑1 (MQ=255)
                                                                                                                            cTCAAATTTGGCGCGATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGTTGCTCAg                  >  2:1654474/1‑148 (MQ=255)
                                                                                                                                       cgcgATTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGTTGCTCAGGCACTGGGTCg       >  2:279030/1‑148 (MQ=255)
                                                                                                                                           aTTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGTTGCTCAGGCACTGGGTCGCCCg   >  1:3503936/1‑148 (MQ=255)
                                                                                                                                           aTTAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGTTGCTCAGGCACTGGGTCGCCCg   <  2:2519540/148‑1 (MQ=255)
                                                                                                                                            ttAGCTACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGTTGCTCAGGCACTGGGTCGCCCGc  >  2:1922269/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
AAAAGCAGAGCAACGGGCGATAAAAATGTTCGCTTATGCCAATAGTCATCCAGGATTGTGGAAAGTCGGGATGATGGCCGGTGCTCATGCGGCAAGCTGGTTTATCAATGGCGGCAAAACACCACTCAAATTTGGCGCGATTAGCGACTGGATGGAAGCACGCGATCTTCCTGAAGCTGACGGAGAGAGTTTCCGTAGTTGGTTTAAGAAACATCAGGCGCAGGAGAAAAAGAATGGATAATCGAGGCGAATTTTTGAATAACGTTGCTCAGGCACTGGGTCGCCCGC  >  CP009273/319236‑319523

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: