Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,281,042 C→A 9.0% R150S (CGT→AGT)  narJ → molybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,281,0420CA 9.0% 768.7 / 34.1 221R150S (CGT→AGT) narJmolybdenum‑cofactor‑assembly chaperone subunit (delta subunit) of nitrate reductase 1
Reads supporting (aligned to +/- strand):  ref base C (100/101);  new base A (13/7);  total (113/108)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 224 total aligned reads displayed.

GCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGT  >  CP009273/1280894‑1281190
                                                                                                                                                    |                                                                                                                                                    
gcaCGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGc                                                                                                                                                      <  1:3887846/149‑1 (MQ=255)
  aCGGCTTGCAGTTACACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGt                                                                                                                                                    <  2:3941271/149‑1 (MQ=255)
  aCGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGt                                                                                                                                                    <  2:1529522/149‑1 (MQ=255)
  aCGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGt                                                                                                                                                    <  1:3439495/149‑1 (MQ=255)
   cGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTg                                                                                                                                                   <  2:2685951/149‑1 (MQ=255)
    gtcTTGCAGTTACACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGa                                                                                                                                                  <  1:1873467/147‑1 (MQ=255)
      cTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGaaa                                                                                                                                                <  1:2183829/149‑1 (MQ=255)
        tGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGc                                                                                                                                              <  2:4067408/149‑1 (MQ=255)
         gCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGcc                                                                                                                                             <  1:3849190/149‑1 (MQ=255)
            gTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGtt                                                                                                                                          <  1:1605246/149‑1 (MQ=255)
            gTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGtt                                                                                                                                          <  1:842616/149‑1 (MQ=255)
             ttAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGTAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTa                                                                                                                                         <  2:3827535/149‑1 (MQ=255)
             ttAAACAGCCGCGAAGTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTa                                                                                                                                         <  2:3825528/149‑1 (MQ=255)
               aaaCAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATg                                                                                                                                       >  2:2085675/1‑149 (MQ=255)
                aaCAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGc                                                                                                                                      <  2:2821357/149‑1 (MQ=255)
                 aCAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGcc                                                                                                                                     <  1:2443212/149‑1 (MQ=255)
                   aGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGt                                                                                                                                   <  2:365077/149‑1 (MQ=255)
                    gCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTg                                                                                                                                  >  1:350750/1‑149 (MQ=255)
                     ccGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGAGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGc                                                                                                                                 <  1:746692/149‑1 (MQ=255)
                      cgcgAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCt                                                                                                                                <  1:213125/149‑1 (MQ=255)
                         gAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGtt                                                                                                                             <  2:2268615/149‑1 (MQ=255)
                          aaTTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGttt                                                                                                                            <  2:624688/149‑1 (MQ=255)
                          aaTTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGttt                                                                                                                            >  1:3933247/1‑149 (MQ=255)
                           aTTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTg                                                                                                                           >  1:2004175/1‑149 (MQ=255)
                            ttGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGa                                                                                                                          >  1:1805964/1‑149 (MQ=255)
                             tGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGAt                                                                                                                         <  2:772254/149‑1 (MQ=255)
                              gCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATc                                                                                                                        >  1:1385467/1‑149 (MQ=255)
                              gCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATc                                                                                                                        <  2:331992/149‑1 (MQ=255)
                               ccGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATct                                                                                                                       <  2:2930046/149‑1 (MQ=255)
                                cGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATctg                                                                                                                      <  1:2224331/149‑1 (MQ=255)
                                      aTCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATctgctgctg                                                                                                                >  2:494926/1‑149 (MQ=255)
                                      aTCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATctgctgctg                                                                                                                >  1:3019828/1‑149 (MQ=255)
                                          gCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAAc                                                                                                            >  1:3531631/1‑149 (MQ=255)
                                           ccGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACt                                                                                                           <  2:3692626/149‑1 (MQ=255)
                                             gCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTgg                                                                                                         >  2:3161506/1‑149 (MQ=255)
                                              cTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGc                                                                                                        <  1:1915376/149‑1 (MQ=255)
                                                 tATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGaa                                                                                                     >  1:1200332/1‑149 (MQ=255)
                                                 tATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGaa                                                                                                     <  1:3219439/149‑1 (MQ=255)
                                                 tATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGaa                                                                                                     >  2:292095/1‑149 (MQ=255)
                                                  aTCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAAt                                                                                                    <  2:2283205/149‑1 (MQ=255)
                                                  aTCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAAt                                                                                                    <  2:3089946/149‑1 (MQ=255)
                                                   tCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATa                                                                                                   <  2:849060/149‑1 (MQ=255)
                                                               ggCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACa                                                                                       >  1:779602/1‑149 (MQ=255)
                                                                 cgcAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGc                                                                                     <  2:1768789/149‑1 (MQ=255)
                                                                  gcAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCg                                                                                    >  2:1503661/1‑149 (MQ=255)
                                                                   cAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGa                                                                                   <  2:779602/149‑1 (MQ=255)
                                                                    aGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGAc                                                                                  >  1:720745/1‑149 (MQ=255)
                                                                     gCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACa                                                                                 >  1:2839654/1‑149 (MQ=255)
                                                                     gCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACa                                                                                 >  2:1753163/1‑149 (MQ=255)
                                                                         ccGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGt                                                                             <  2:2050480/149‑1 (MQ=255)
                                                                         ccGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGt                                                                             >  2:3749187/1‑149 (MQ=255)
                                                                             aaaGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCg                                                                         >  2:3363943/1‑149 (MQ=255)
                                                                             aaaGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCg                                                                         >  2:3704827/1‑149 (MQ=255)
                                                                              aaGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCgg                                                                        <  2:2658150/149‑1 (MQ=255)
                                                                              aaGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCgg                                                                        >  2:1030665/1‑149 (MQ=255)
                                                                               aGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGa                                                                       >  2:4110147/1‑149 (MQ=255)
                                                                                   aaGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGaaaaa                                                                   >  1:3650559/1‑149 (MQ=255)
                                                                                    aGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGaaaaaa                                                                  >  2:1885229/1‑149 (MQ=255)
                                                                                    aGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGaaaaaa                                                                  <  1:506541/149‑1 (MQ=255)
                                                                                    aGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGaaaaaa                                                                  >  1:456196/1‑149 (MQ=255)
                                                                                    aGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGaaaaaa                                                                  <  1:1086191/149‑1 (MQ=255)
                                                                                     gCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTAAGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAAt                                                                 >  1:134732/1‑149 (MQ=255)
                                                                                        gTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGc                                                              <  1:373965/149‑1 (MQ=255)
                                                                                         tGGAAGTTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGATGTTTGATCTGCTGCTGAAACTGGCGAATACCGATATCGACAGCGACAAAGTGGCGGAAAAAATTGcc                                                             >  2:1137511/1‑149 (MQ=255)
                                                                                         tGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGcc                                                             >  2:874931/1‑149 (MQ=255)
                                                                                         tGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGcc                                                             >  1:2415862/1‑149 (MQ=255)
                                                                                             aGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCcgacg                                                         >  2:1983467/1‑149 (MQ=255)
                                                                                                      aGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGcgcgcg                                                >  1:1007822/1‑149 (MQ=255)
                                                                                                      aGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGcgcgcg                                                >  2:415205/1‑149 (MQ=255)
                                                                                                          atCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGACGCGCGCgatga                                            >  2:1536791/1‑149 (MQ=255)
                                                                                                          atCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCgatga                                            >  2:469950/1‑149 (MQ=255)
                                                                                                           tCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCgatgat                                           >  2:1672734/1‑149 (MQ=255)
                                                                                                               gcCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATAcgc                                       >  2:3443205/1‑149 (MQ=255)
                                                                                                                 cGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATAcgccg                                     >  1:2534536/1‑149 (MQ=255)
                                                                                                                   aTTCTGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCa                                   >  1:959562/1‑149 (MQ=255)
                                                                                                                       tGGCATTGCTGAGCGCGCGTCTGCAACAGAGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGgcg                               >  1:1237347/1‑149 (MQ=255)
                                                                                                                        ggCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGgcgc                              <  1:1302805/149‑1 (MQ=255)
                                                                                                                           aTTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTgg                           >  2:3339979/1‑149 (MQ=255)
                                                                                                                            ttGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGa                          >  1:2275273/1‑149 (MQ=255)
                                                                                                                            ttGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGa                          >  1:2379612/1‑149 (MQ=255)
                                                                                                                                tGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCt                      <  2:2509148/149‑1 (MQ=255)
                                                                                                                                  aGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGt                    <  1:2820765/149‑1 (MQ=255)
                                                                                                                                  aGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGt                    <  1:2923750/149‑1 (MQ=255)
                                                                                                                                   gcgcgcGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTc                   >  1:4018001/1‑149 (MQ=255)
                                                                                                                                   gcgcgcGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTc                   >  1:3103383/1‑149 (MQ=255)
                                                                                                                                    cgcgcgTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCt                  <  2:2739025/149‑1 (MQ=255)
                                                                                                                                    cgcgcgTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCt                  <  2:2361773/149‑1 (MQ=255)
                                                                                                                                     gcgcgTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGACGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTg                 >  2:3767756/1‑149 (MQ=255)
                                                                                                                                       gcgTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTggg               >  2:3027482/1‑149 (MQ=255)
                                                                                                                                       gcgTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTggg               <  2:2873851/149‑1 (MQ=255)
                                                                                                                                         gTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGgaa             <  1:494926/149‑1 (MQ=255)
                                                                                                                                         gTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGgaa             >  1:1442802/1‑149 (MQ=255)
                                                                                                                                          tCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGACGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGgaag            >  2:1636904/1‑149 (MQ=255)
                                                                                                                                              cAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAagaaga        <  2:1805964/149‑1 (MQ=255)
                                                                                                                                              cAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGAGGCTGTCTGGGAagaaga        <  2:2395488/149‑1 (MQ=255)
                                                                                                                                                aCAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGc      <  1:4110147/149‑1 (MQ=255)
                                                                                                                                                 cAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCa     >  2:3177126/1‑149 (MQ=255)
                                                                                                                                                   gCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAgg   >  2:4105005/1‑149 (MQ=255)
                                                                                                                                                   gCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAgg   <  2:3244455/149‑1 (MQ=255)
                                                                                                                                                    cGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGt  <  2:43920/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCACGGCTTGCAGTTAAACAGCCGCGAATTGCCGGACCATCTGCCGCTGTATCTGGAGTACCTGGCGCAGCTGCCGCAAAGCGAAGCCGTGGAAGGTTTGAAAGATATCGCGCCGATTCTGGCATTGCTGAGCGCGCGTCTGCAACAGCGTGAAAGCCGTTATGCCGTGCTGTTTGATCTGCTGCTGAAACTGGCGAATACCGCTATCGACAGCGACAAAGTGGCGGAAAAAATTGCCGACGAAGCGCGCGATGATACGCCGCAGGCGCTGGATGCTGTCTGGGAAGAAGAGCAGGT  >  CP009273/1280894‑1281190

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: