Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 788,902 G→T 89.5% R115S (CGT→AGT)  modF ← fused molybdate transporter subunits of ABC superfamily: ATP‑binding components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273788,9020GT89.5% 895.8 / 94.8 295R115S (CGT→AGT) modFfused molybdate transporter subunits of ABC superfamily: ATP‑binding components
Reads supporting (aligned to +/- strand):  ref base G (16/15);  new base T (125/139);  total (141/154)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 295 total aligned reads displayed.

AACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGAT  >  CP009273/788757‑789046
                                                                                                                                                 |                                                                                                                                                
aaCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCgg                                                                                                                                               <  1:1244043/149‑1 (MQ=255)
 aCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGt                                                                                                                                              >  1:1875051/1‑149 (MQ=255)
  cGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTg                                                                                                                                             >  2:284259/1‑149 (MQ=255)
    gCTCATCAAGAATCTACAACTCAGGCTCCGACATCAGAGCCTGACACTGCAGGGATTTTATCGTCTCTCCTGTTGAAAGGTATTTACTGCGTCGGTCTAGGAGGGCGGTAATACCTAACTGCTGCGCCAGTTGCATGCAACTCGGTGCa                                                                                                                                           <  2:3913149/149‑1 (MQ=255)
    gCTCATCAAGAATCCACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCa                                                                                                                                           <  2:1295595/149‑1 (MQ=255)
    gCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGACAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCa                                                                                                                                           >  1:437919/1‑149 (MQ=255)
     cTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCAt                                                                                                                                          <  1:932240/149‑1 (MQ=255)
       catcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATcc                                                                                                                                        <  1:334692/149‑1 (MQ=255)
         tcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCtt                                                                                                                                      <  1:2183493/149‑1 (MQ=255)
          caAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCttt                                                                                                                                     <  2:4915232/149‑1 (MQ=255)
            aGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTAc                                                                                                                                   <  1:2566/149‑1 (MQ=255)
                tCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCa                                                                                                                               <  2:2761559/149‑1 (MQ=255)
                   acaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATcc                                                                                                                            <  2:956138/149‑1 (MQ=255)
                   acaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATcc                                                                                                                            <  1:1776427/149‑1 (MQ=255)
                     aaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTg                                                                                                                          >  2:2257033/1‑149 (MQ=255)
                      aGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGa                                                                                                                         >  1:1686216/1‑149 (MQ=255)
                         cAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATg                                                                                                                      <  1:3114480/149‑1 (MQ=255)
                            gCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTCATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtc                                                                                                                   >  2:1359673/1‑149 (MQ=255)
                            gCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtc                                                                                                                   <  2:1568149/149‑1 (MQ=255)
                             cTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtct                                                                                                                  <  1:2483253/149‑1 (MQ=255)
                               ccGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCa                                                                                                                <  2:3042254/149‑1 (MQ=255)
                               ccGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCa                                                                                                                >  1:4857575/1‑149 (MQ=255)
                               ccGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCa                                                                                                                <  2:1479357/149‑1 (MQ=255)
                                 gACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGc                                                                                                              >  2:3649916/1‑149 (MQ=255)
                                 gACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGc                                                                                                              <  2:1176755/149‑1 (MQ=255)
                                      cAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAg                                                                                                         <  2:1495480/149‑1 (MQ=255)
                                       aGCGCCTGACACAGCCGGGTTTTTCGCGTCTCGCAAGTGGATTGTTATTTAAATCGTCGGTCGAGGAGGTCGCTATTACGTTAATGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGt                                                                                                        <  2:869067/149‑1 (MQ=255)
                                        gcgcCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg                                                                                                       >  2:1888846/1‑149 (MQ=255)
                                           ccTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGt                                                                                                    <  1:4187225/149‑1 (MQ=255)
                                            cTGACACATCAGGGTTTTTCGCGTCTCGCCAGTTGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGAGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTc                                                                                                   <  2:1083891/149‑1 (MQ=255)
                                            cTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTc                                                                                                   >  1:4024985/1‑149 (MQ=255)
                                            cTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTc                                                                                                   <  2:4270804/149‑1 (MQ=255)
                                              gACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCg                                                                                                 >  1:3838873/1‑149 (MQ=255)
                                                cacaGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGgt                                                                                               >  2:3456907/1‑149 (MQ=255)
                                                 acaacaGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGgtg                                                                                              >  2:2689685/1‑149 (MQ=255)
                                                   agcagGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTc                                                                                            >  1:2793108/1‑149 (MQ=255)
                                                      agGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATc                                                                                         <  2:3698047/149‑1 (MQ=255)
                                                        ggTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCtt                                                                                       <  2:180860/149‑1 (MQ=255)
                                                        ggTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCtt                                                                                       <  2:753211/149‑1 (MQ=255)
                                                           ttttCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGc                                                                                    <  2:843002/149‑1 (MQ=255)
                                                           ttttCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGc                                                                                    <  2:1257569/149‑1 (MQ=255)
                                                            tttCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGcc                                                                                   >  2:2124501/1‑149 (MQ=255)
                                                            tttCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGcc                                                                                   <  2:3275152/149‑1 (MQ=255)
                                                              tCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAg                                                                                 <  1:5174582/149‑1 (MQ=255)
                                                                     tcGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAg                                                                          >  2:3347320/1‑149 (MQ=255)
                                                                      cGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGc                                                                         >  2:3675990/1‑149 (MQ=255)
                                                                       gCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCa                                                                        <  2:2700722/149‑1 (MQ=255)
                                                                        ccAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCat                                                                       >  2:3749989/1‑149 (MQ=255)
                                                                         cAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCata                                                                      >  2:418916/1‑149 (MQ=255)
                                                                          aGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCatat                                                                     >  1:2604333/1‑149 (MQ=255)
                                                                           gTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATc                                                                    >  1:4408841/1‑149 (MQ=255)
                                                                             ggAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCgg                                                                  >  1:3274089/1‑149 (MQ=255)
                                                                               aaaGGTATGTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgtg                                                                <  2:5092536/149‑1 (MQ=255)
                                                                                aaGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgtgt                                                               >  1:224702/1‑149 (MQ=255)
                                                                                   gTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTAt                                                            <  2:2517534/149‑1 (MQ=255)
                                                                                    tATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTAtt                                                           <  2:376376/149‑1 (MQ=255)
                                                                                      tttAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTcc                                                         <  1:2774502/149‑1 (MQ=255)
                                                                                      tttAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTcc                                                         <  2:3838873/149‑1 (MQ=255)
                                                                                         aaaGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCt                                                      <  1:464954/149‑1 (MQ=255)
                                                                                          aaGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTg                                                     >  2:553794/1‑149 (MQ=255)
                                                                                          aaGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTg                                                     >  2:688544/1‑149 (MQ=255)
                                                                                           aGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGc                                                    >  2:417975/1‑149 (MQ=255)
                                                                                           aGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGc                                                    >  2:4989614/1‑149 (MQ=255)
                                                                                             cGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCa                                                  <  2:4615895/149‑1 (MQ=255)
                                                                                               tcggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCAtt                                                >  1:1068504/1‑149 (MQ=255)
                                                                                               tcggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCCTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATc                                                                    <  1:644704/129‑1 (MQ=255)
                                                                                                cggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtc                                               <  1:5103279/149‑1 (MQ=255)
                                                                                                  gtcgAGGAGGGCGGTAATACGGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgt                                             >  1:213121/1‑149 (MQ=255)
                                                                                                  gtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgt                                             <  2:261657/149‑1 (MQ=255)
                                                                                                  gtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgt                                             <  1:2466263/149‑1 (MQ=255)
                                                                                                  gtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgt                                             >  2:1615664/1‑149 (MQ=255)
                                                                                                     gaggagGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTACATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGc                                          <  2:4349502/149‑1 (MQ=255)
                                                                                                      aggaggGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCt                                         >  1:1208660/1‑149 (MQ=255)
                                                                                                      aggaggGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCt                                         >  1:2938230/1‑149 (MQ=255)
                                                                                                       ggaggGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTg                                        >  2:2199525/1‑149 (MQ=255)
                                                                                                         aggGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGAc                                      <  1:38076/149‑1 (MQ=255)
                                                                                                          gggCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACg                                     <  1:2683611/149‑1 (MQ=255)
                                                                                                          gggCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACg                                     >  2:4733246/1‑149 (MQ=255)
                                                                                                                tAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGtttt                               <  1:1062194/149‑1 (MQ=255)
                                                                                                                      cGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAAt                         >  1:2886165/1‑149 (MQ=255)
                                                                                                                      cGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAAt                         <  1:1509670/149‑1 (MQ=255)
                                                                                                                      cGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAAt                         <  1:773770/149‑1 (MQ=255)
                                                                                                                       gAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAAtt                        <  1:4554725/149‑1 (MQ=255)
                                                                                                                         aCTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGc                      <  2:1223033/149‑1 (MQ=255)
                                                                                                                          ctgctgCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCt                     <  2:2404833/149‑1 (MQ=255)
                                                                                                                          ctgctgCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCt                     <  2:1695501/149‑1 (MQ=255)
                                                                                                                          ctgctgCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCt                     <  1:4986253/149‑1 (MQ=255)
                                                                                                                               gcgcCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAg                <  2:4816783/149‑1 (MQ=255)
                                                                                                                                    aGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgagag           <  1:3549471/149‑1 (MQ=255)
                                                                                                                                    aGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgagag           >  2:636399/1‑149 (MQ=255)
                                                                                                                                    aGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgagag           <  1:3318501/149‑1 (MQ=255)
                                                                                                                                     gTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgagaga          <  1:3413473/149‑1 (MQ=255)
                                                                                                                                      ttGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGAc         <  1:4281597/149‑1 (MQ=255)
                                                                                                                                      ttGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGAc         <  2:3435235/149‑1 (MQ=255)
                                                                                                                                      ttGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGAc         <  1:2384050/149‑1 (MQ=255)
                                                                                                                                       tgcatgcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACg        >  2:4845113/1‑149 (MQ=255)
                                                                                                                                            gcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGa   <  1:4040513/149‑1 (MQ=255)
                                                                                                                                            gcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGa   <  2:4024985/149‑1 (MQ=255)
                                                                                                                                            gcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGa   <  1:1421311/149‑1 (MQ=255)
                                                                                                                                             caACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGAt  <  1:3908538/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGAT  >  CP009273/788757‑789046

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: