Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 788,902 G→T 94.8% 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,9020GT94.8% 1110.0 / 41.9 342R115S (CGT→AGT) modFfused molybdate transporter subunits of ABC superfamily: ATP‑binding components
Reads supporting (aligned to +/- strand):  ref base G (6/12);  new base T (171/153);  total (177/165)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 343 total aligned reads displayed.

CGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGT  >  CP009273/788755‑789048
                                                                                                                                                   |                                                                                                                                                  
cGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAActcg                                                                                                                                                  >  1:253508/1‑150 (MQ=255)
  aaCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGt                                                                                                                                                >  1:616293/1‑150 (MQ=255)
   aCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTg                                                                                                                                               >  1:3338588/1‑150 (MQ=255)
    cGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGGAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGc                                                                                                                                              >  1:448678/1‑150 (MQ=255)
     ggCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCa                                                                                                                                             >  2:3168477/1‑150 (MQ=255)
     ggCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCa                                                                                                                                             <  2:3577355/150‑1 (MQ=255)
      gCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCAt                                                                                                                                            <  2:2072138/150‑1 (MQ=255)
       cTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATc                                                                                                                                           >  1:3279309/1‑150 (MQ=255)
       cTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATc                                                                                                                                           >  1:3749789/1‑150 (MQ=255)
         catcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCt                                                                                                                                         <  1:4352686/150‑1 (MQ=255)
           tcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCttt                                                                                                                                       <  1:809230/150‑1 (MQ=255)
           tcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCttt                                                                                                                                       >  2:1528529/1‑150 (MQ=255)
            caAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTa                                                                                                                                      >  1:2600725/1‑150 (MQ=255)
                aaTCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTGGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTc                                                                                                                                  >  1:3926127/1‑150 (MQ=255)
                   caacaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATc                                                                                                                               <  2:437267/150‑1 (MQ=255)
                      ctaGTCAGGCTCCGACATCTGCGCCTGACACAGCTGTGTTTTTCGCGTCTCGCGAGTGGAAAGTTATTTTAATCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGAATCCTTTACTTCATCCTg                                                                                                                            <  1:4047255/148‑1 (MQ=255)
                           cAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGa                                                                                                                       <  2:686965/150‑1 (MQ=255)
                             ggCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtc                                                                                                                     >  1:2398260/1‑150 (MQ=255)
                             ggCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtc                                                                                                                     <  2:2122982/150‑1 (MQ=255)
                              gCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtct                                                                                                                    >  2:1572241/1‑150 (MQ=255)
                              gCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtct                                                                                                                    <  2:5089417/150‑1 (MQ=255)
                              gCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtct                                                                                                                    <  1:3438088/150‑1 (MQ=255)
                                tCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCa                                                                                                                  >  2:1910747/1‑150 (MQ=255)
                                tCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCa                                                                                                                  <  1:4691008/150‑1 (MQ=255)
                                tCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCa                                                                                                                  >  2:3915880/1‑150 (MQ=255)
                                   gACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGcc                                                                                                               <  2:3876574/150‑1 (MQ=255)
                                   gACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGcc                                                                                                               <  2:2638494/150‑1 (MQ=255)
                                    aCATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCg                                                                                                              >  1:4122978/1‑150 (MQ=255)
                                     cATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGt                                                                                                             >  2:1168524/1‑150 (MQ=255)
                                      aTCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTa                                                                                                            <  1:1000150/150‑1 (MQ=255)
                                      aTCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTa                                                                                                            <  1:4841189/150‑1 (MQ=255)
                                       tCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAg                                                                                                           <  2:2870947/150‑1 (MQ=255)
                                       tCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTAGAGGAGGGCGGTAATACAGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAg                                                                                                           <  2:3695904/150‑1 (MQ=255)
                                        cAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGt                                                                                                          >  1:4236713/1‑150 (MQ=255)
                                        cAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGt                                                                                                          >  1:1488943/1‑150 (MQ=255)
                                         aGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg                                                                                                         >  2:1483975/1‑150 (MQ=255)
                                         aGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg                                                                                                         >  1:1531625/1‑150 (MQ=255)
                                         aGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg                                                                                                         >  2:4043642/1‑150 (MQ=255)
                                         aGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg                                                                                                         >  2:1888484/1‑150 (MQ=255)
                                         aGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGGAGTg                                                                                                         >  2:3242201/1‑150 (MQ=255)
                                                 acacAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGgt                                                                                                 >  1:2646113/1‑150 (MQ=255)
                                                  cacaGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGgtg                                                                                                <  2:2537795/150‑1 (MQ=255)
                                                   acaGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGgtgt                                                                                               >  2:1028361/1‑150 (MQ=255)
                                                    cagcagGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTc                                                                                              <  2:376846/150‑1 (MQ=255)
                                                       cagGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATc                                                                                           >  1:270591/1‑150 (MQ=255)
                                                         gggTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCtt                                                                                         <  2:3113087/150‑1 (MQ=255)
                                                         gggTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCtt                                                                                         >  1:229010/1‑150 (MQ=255)
                                                          ggTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTc                                                                                        <  1:1150734/150‑1 (MQ=255)
                                                           gTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCg                                                                                       <  1:939260/150‑1 (MQ=255)
                                                                tCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAgg                                                                                  >  2:4468561/1‑150 (MQ=255)
                                                                   cgTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGcc                                                                               >  1:2091654/1‑150 (MQ=255)
                                                                     tctcGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGa                                                                             <  2:1781307/150‑1 (MQ=255)
                                                                        cGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCa                                                                          >  2:1056142/1‑150 (MQ=255)
                                                                           cAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTCGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCatat                                                                       >  1:4171427/1‑150 (MQ=255)
                                                                           cAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCatat                                                                       <  1:1988571/150‑1 (MQ=255)
                                                                              tGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCgg                                                                    >  2:4076632/1‑150 (MQ=255)
                                                                               ggAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgt                                                                   >  1:1413297/1‑150 (MQ=255)
                                                                               ggAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgt                                                                   >  2:235807/1‑150 (MQ=255)
                                                                                 aaaGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgtgt                                                                 >  2:1299885/1‑150 (MQ=255)
                                                                                 aaaGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgtgt                                                                 >  1:4788525/1‑150 (MQ=255)
                                                                                   aGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTa                                                               >  2:3482423/1‑150 (MQ=255)
                                                                                      tATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTc                                                            <  1:539243/150‑1 (MQ=255)
                                                                                       aTTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTcc                                                           <  1:319872/150‑1 (MQ=255)
                                                                                           aaaGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTg                                                       <  1:2232189/150‑1 (MQ=255)
                                                                                            aaGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGc                                                      <  2:4633854/150‑1 (MQ=255)
                                                                                            aaGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGc                                                      >  1:1237956/1‑150 (MQ=255)
                                                                                              gCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCa                                                    <  2:1622031/150‑1 (MQ=255)
                                                                                               cGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCAt                                                   <  2:4279023/150‑1 (MQ=255)
                                                                                                gtcggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCAtt                                                  >  2:1709284/1‑150 (MQ=255)
                                                                                                 tcggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtc                                                 <  1:3089535/150‑1 (MQ=255)
                                                                                                  cggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcg                                                >  2:2423206/1‑150 (MQ=255)
                                                                                                  cggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcg                                                >  1:1250920/1‑150 (MQ=255)
                                                                                                    gtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgtc                                              >  2:357447/1‑150 (MQ=255)
                                                                                                     tcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgtcg                                             >  1:1695718/1‑150 (MQ=255)
                                                                                                       gaggagGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCt                                           >  1:109099/1‑150 (MQ=255)
                                                                                                       gaggagGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCt                                           >  2:5142260/1‑150 (MQ=255)
                                                                                                          gaggGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGAc                                        <  2:74481/150‑1 (MQ=255)
                                                                                                               cGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGtt                                   <  2:220711/150‑1 (MQ=255)
                                                                                                                ggTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGttt                                  >  2:2771685/1‑150 (MQ=255)
                                                                                                                  tAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGttttt                                <  1:3169790/150‑1 (MQ=255)
                                                                                                                      aCCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCaa                            >  1:380770/1‑150 (MQ=255)
                                                                                                                       ccGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTTAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAAt                           <  2:1829769/150‑1 (MQ=255)
                                                                                                                        cGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAAtt                          >  1:5071250/1‑150 (MQ=255)
                                                                                                                         gAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGATGAGTTTTTGCAATTg                         <  2:4824317/150‑1 (MQ=255)
                                                                                                                         gAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTg                         >  2:4286556/1‑150 (MQ=255)
                                                                                                                            ctgctgCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTc                      <  1:752555/150‑1 (MQ=255)
                                                                                                                              gctgcGCCAGGTGCATGCAACTCGGTGCATCCTGTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGa                    <  2:4248408/150‑1 (MQ=255)
                                                                                                                                tgcGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAg                  <  2:1111085/150‑1 (MQ=255)
                                                                                                                                    ccAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgaga              <  2:299406/150‑1 (MQ=255)
                                                                                                                                    ccAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgaga              >  1:2954400/1‑150 (MQ=255)
                                                                                                                                     cAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgagag             >  1:5148633/1‑150 (MQ=255)
                                                                                                                                     cAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgagag             >  2:4906219/1‑150 (MQ=255)
                                                                                                                                      aGTTTCATGCAACTCGTTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATATTCGACTTGGCAGAGCTTATCGGTGTTATTCCGATGCAATTAGTCTCTGACGAGTTTTTTCAATTGCTCTAAGgagaga            <  2:4964914/150‑1 (MQ=255)
                                                                                                                                      aGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCACTTGCTCGAAGgagaga            >  2:904188/1‑150 (MQ=255)
                                                                                                                                      aGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgagaga            >  1:3054947/1‑150 (MQ=255)
                                                                                                                                       gTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGAc           <  1:1659012/150‑1 (MQ=255)
                                                                                                                                        ttGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACg          <  1:3568345/150‑1 (MQ=255)
                                                                                                                                             tgcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGa     <  1:4419828/150‑1 (MQ=255)
                                                                                                                                                aaCTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGAtgt  >  1:713801/1‑150 (MQ=255)
                                                                                                                                                aaCTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGAtgt  >  2:2817844/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGT  >  CP009273/788755‑789048

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: