Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,275,339 G→T 14.2% R7L (CGG→CTG)  narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,275,3390GT14.2% 987.4 / 112.9 307R7L (CGG→CTG) narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (156/107);  new base T (32/11);  total (189/118)
Fisher's exact test for biased strand distribution p-value = 6.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

Only 100 of 312 total aligned reads displayed.

ACTTTATTTTTCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCT  >  CP009273/1275190‑1275488
                                                                                                                                                     |                                                                                                                                                     
   ttATTTTTCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGtt                                                                                                                                                    >  1:607871/1‑150 (MQ=255)
    tATTTTTCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGttt                                                                                                                                                   >  1:3071104/1‑150 (MQ=255)
    tATTTTTCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGttt                                                                                                                                                   >  1:2677284/1‑150 (MQ=255)
      tttttCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCg                                                                                                                                                 >  1:2874655/1‑150 (MQ=255)
      tttttCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCg                                                                                                                                                 >  1:2327315/1‑150 (MQ=255)
          tCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTAc                                                                                                                                             >  2:4095640/1‑150 (MQ=255)
           cATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACt                                                                                                                                            >  2:3670160/1‑150 (MQ=255)
           cATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACt                                                                                                                                            <  2:3747909/150‑1 (MQ=255)
             tCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTc                                                                                                                                          <  2:1435389/150‑1 (MQ=255)
              ccccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCa                                                                                                                                         <  1:2450544/150‑1 (MQ=255)
              ccccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCa                                                                                                                                         <  2:509924/150‑1 (MQ=255)
                ccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCaaa                                                                                                                                       >  1:4628365/1‑150 (MQ=255)
                ccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCaaa                                                                                                                                       >  1:4007570/1‑150 (MQ=255)
                  aTCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACa                                                                                                                                     >  2:2366121/1‑150 (MQ=255)
                     aCTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGaa                                                                                                                                  >  2:701078/1‑150 (MQ=255)
                     aCTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAAGTCCTGGACCGGTTTCGCTACTTCAAACAGaa                                                                                                                                  >  2:1276320/1‑150 (MQ=255)
                      ctctTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAg                                                                                                                                 >  1:1815587/1‑150 (MQ=255)
                       tctTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAgg                                                                                                                                >  2:385890/1‑150 (MQ=255)
                          tGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTg                                                                                                                             <  1:3679211/150‑1 (MQ=255)
                          tGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTg                                                                                                                             >  1:726353/1‑150 (MQ=255)
                               gTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAAcc                                                                                                                        >  1:2713930/1‑150 (MQ=255)
                                ttATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCt                                                                                                                       <  1:3506858/150‑1 (MQ=255)
                                 tATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCtt                                                                                                                      >  2:2755810/1‑150 (MQ=255)
                                  aTCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCttt                                                                                                                     >  2:4868403/1‑150 (MQ=255)
                                  aTCAAGTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCttt                                                                                                                     <  2:3576253/150‑1 (MQ=255)
                                   tCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTg                                                                                                                    <  2:1509557/150‑1 (MQ=255)
                                     aaTTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGcc                                                                                                                  >  2:3247563/1‑150 (MQ=255)
                                     aaTTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGcc                                                                                                                  <  2:3393430/150‑1 (MQ=255)
                                     aaTTCCCACGCAGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTACATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGcc                                                                                                                  >  2:4793746/1‑150 (MQ=255)
                                       ttCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGa                                                                                                                <  2:3541983/150‑1 (MQ=255)
                                        tCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGAt                                                                                                               >  2:4503919/1‑150 (MQ=255)
                                        tCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGAt                                                                                                               >  1:4080775/1‑150 (MQ=255)
                                        tCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGAt                                                                                                               >  1:4081815/1‑150 (MQ=255)
                                         cccACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATg                                                                                                              <  2:4776719/150‑1 (MQ=255)
                                                tGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATgg                                                                                                       <  2:3089130/150‑1 (MQ=255)
                                                 gTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGc                                                                                                      >  2:1464849/1‑150 (MQ=255)
                                                   ttCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCa                                                                                                    <  1:3770074/150‑1 (MQ=255)
                                                   ttCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCa                                                                                                    >  1:1806015/1‑150 (MQ=255)
                                                       gagCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCtt                                                                                                >  2:37030/1‑150 (MQ=255)
                                                        agCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTtc                                                                                               >  2:4356711/1‑150 (MQ=255)
                                                          cGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTtctc                                                                                             >  1:2505022/1‑150 (MQ=255)
                                                           gTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCa                                                                                            >  2:1529985/1‑150 (MQ=255)
                                                           gTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCa                                                                                            >  2:381234/1‑150 (MQ=255)
                                                           gTTACCTTGCCATTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTAAAACAGAAGGGTTAAACCTTTGCCGAAGGCCATGGCAAGCTACTCa                                                                                            >  2:4696498/1‑150 (MQ=255)
                                                            ttACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCaa                                                                                           >  1:2386981/1‑150 (MQ=255)
                                                              aCCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATa                                                                                         >  1:4216932/1‑150 (MQ=255)
                                                               ccTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATAc                                                                                        >  2:2164917/1‑150 (MQ=255)
                                                               ccTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATAc                                                                                        >  2:593883/1‑150 (MQ=255)
                                                               ccTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATAc                                                                                        <  1:1926205/150‑1 (MQ=255)
                                                                cTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATAcc                                                                                       <  2:3017477/150‑1 (MQ=255)
                                                                 ttGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCa                                                                                      <  2:2663565/150‑1 (MQ=255)
                                                                  tGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCaa                                                                                     >  1:1102357/1‑150 (MQ=255)
                                                                   gCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAAc                                                                                    >  1:4517727/1‑150 (MQ=255)
                                                                     ccTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCg                                                                                  >  1:3625516/1‑150 (MQ=255)
                                                                        tAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGa                                                                               <  1:1017529/150‑1 (MQ=255)
                                                                         aaaCATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGAc                                                                              <  1:1255308/150‑1 (MQ=255)
                                                                         aaaCATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGAc                                                                              <  2:3183477/150‑1 (MQ=255)
                                                                         aaaCATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGAc                                                                              >  1:2216654/1‑150 (MQ=255)
                                                                           aCATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTg                                                                            <  2:1930223/150‑1 (MQ=255)
                                                                            cATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTgg                                                                           >  2:4476440/1‑150 (MQ=255)
                                                                              ttAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGgga                                                                         <  1:303651/150‑1 (MQ=255)
                                                                              ttAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGgga                                                                         >  1:2289500/1‑150 (MQ=255)
                                                                                aGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGggagg                                                                       <  2:1145567/150‑1 (MQ=255)
                                                                                aGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGggagg                                                                       >  2:1987829/1‑150 (MQ=255)
                                                                                 gCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGAAGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGggagga                                                                      >  2:2936275/1‑150 (MQ=255)
                                                                                    aTGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAggatgg                                                                   >  2:3002631/1‑150 (MQ=255)
                                                                                     tGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAggatgga                                                                  >  2:1071096/1‑150 (MQ=255)
                                                                                     tGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAggatgga                                                                  >  1:2025027/1‑150 (MQ=255)
                                                                                       tCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGata                                                                <  1:1388829/150‑1 (MQ=255)
                                                                                          aTTTCTCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCg                                                             <  1:920374/150‑1 (MQ=255)
                                                                                            ttATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGcc                                                           >  2:1370289/1‑150 (MQ=255)
                                                                                            ttATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGcc                                                           <  1:4356726/150‑1 (MQ=255)
                                                                                            ttATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGcc                                                           >  2:852887/1‑150 (MQ=255)
                                                                                             tATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCa                                                          <  1:4704744/150‑1 (MQ=255)
                                                                                                cAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCg                                                       <  1:4018228/150‑1 (MQ=255)
                                                                                                  gagGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGtt                                                     >  1:518680/1‑150 (MQ=255)
                                                                                                   agGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTg                                                    <  2:2426426/150‑1 (MQ=255)
                                                                                                   agGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTg                                                    <  1:2816488/150‑1 (MQ=255)
                                                                                                    gggCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTgg                                                   >  2:4240054/1‑150 (MQ=255)
                                                                                                       ccGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGgcag                                                <  2:2514959/150‑1 (MQ=255)
                                                                                                        cGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGgcagc                                               >  2:4410163/1‑150 (MQ=255)
                                                                                                        cGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGgcagc                                               <  1:3842250/150‑1 (MQ=255)
                                                                                                         gACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGgcagca                                              >  2:170667/1‑150 (MQ=255)
                                                                                                           cAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATg                                            >  1:1129637/1‑150 (MQ=255)
                                                                                                                tCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACaaa                                       >  1:4237625/1‑150 (MQ=255)
                                                                                                                tCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACaaa                                       >  1:955886/1‑150 (MQ=255)
                                                                                                                              aCCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTAc                         >  1:3454761/1‑150 (MQ=255)
                                                                                                                                cGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTAccc                       >  1:1576447/1‑150 (MQ=255)
                                                                                                                                      gTAAATTCCTGGACCGGTTTCGCTACTTAAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGAAAAAATCGTCCGCTCTACCCACgggg                 <  1:828364/150‑1 (MQ=255)
                                                                                                                                        aaaTTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTa               <  1:3651082/150‑1 (MQ=255)
                                                                                                                                        aaaTTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTa               <  2:4627343/150‑1 (MQ=255)
                                                                                                                                           ttCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAAc            >  2:2804719/1‑150 (MQ=255)
                                                                                                                                             ccTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTg          >  2:1196320/1‑150 (MQ=255)
                                                                                                                                               tGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCa        >  1:4571940/1‑150 (MQ=255)
                                                                                                                                                ggACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCATTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCAc       <  1:443450/150‑1 (MQ=255)
                                                                                                                                                 gACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCAcc      <  1:96610/150‑1 (MQ=255)
                                                                                                                                                  aCCTGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCg     >  1:1749196/1‑150 (MQ=255)
                                                                                                                                                   ccGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCgg    >  2:4110681/1‑150 (MQ=255)
                                                                                                                                                   ccGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCgg    <  2:2583078/150‑1 (MQ=255)
                                                                                                                                                     ggTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCt  >  1:689591/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTTTATTTTTCATCCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATTCCTGGACCGGTTTCGCTACTTCAAACAGAAGGGTGAAACCTTTGCCGATGGGCATGGCCAGCTTCTCAATACCAACCGTGACTGGGAGGATGGATATCGCCAGCGTTGGCAGCATGACAAAATCGTCCGCTCTACCCACGGGGTAAACTGCACCGGCT  >  CP009273/1275190‑1275488

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: