Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 788,902 G→T 100% 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,9020GT97.4% 823.5 / 8.8 248R115S (CGT→AGT) modFfused molybdate transporter subunits of ABC superfamily: ATP‑binding components
Reads supporting (aligned to +/- strand):  ref base G (5/2);  new base T (124/117);  total (129/119)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 255 total aligned reads displayed.

CGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGG  >  CP009273/788755‑789050
                                                                                                                                                   |                                                                                                                                                    
  aaCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGt                                                                                                                                                  >  2:3467241/1‑150 (MQ=255)
  aaCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGt                                                                                                                                                  <  1:4161902/150‑1 (MQ=255)
  aaCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGt                                                                                                                                                  >  1:3624157/1‑150 (MQ=255)
  aaCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGt                                                                                                                                                  <  2:758513/150‑1 (MQ=255)
  aaCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGt                                                                                                                                                  <  2:759010/150‑1 (MQ=255)
   aCGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTg                                                                                                                                                 <  1:3218258/150‑1 (MQ=255)
        tcatcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATcc                                                                                                                                            >  2:3189487/1‑150 (MQ=255)
        tcatcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATcc                                                                                                                                            >  1:3706977/1‑150 (MQ=255)
        tcatcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATcc                                                                                                                                            <  1:1812079/150‑1 (MQ=255)
         catcaAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCt                                                                                                                                           <  1:486637/150‑1 (MQ=255)
          atcaATAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCtt                                                                                                                                          >  1:4461168/1‑150 (MQ=255)
           tctaGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCttt                                                                                                                                         <  1:2053054/147‑1 (MQ=255)
             aaGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTAc                                                                                                                                       >  2:884732/1‑150 (MQ=255)
             aaGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTAc                                                                                                                                       >  1:1912327/1‑150 (MQ=255)
             aaGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTAc                                                                                                                                       >  1:2642175/1‑150 (MQ=255)
                aaTCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACtta                                                                                                                                    <  2:3378641/150‑2 (MQ=255)
                aaTCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTc                                                                                                                                    >  2:1321245/1‑150 (MQ=255)
                    aacaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATcc                                                                                                                                >  2:2001830/1‑150 (MQ=255)
                    aacaaGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATcc                                                                                                                                <  1:4024476/150‑1 (MQ=255)
                        aGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGaa                                                                                                                            >  1:1208952/1‑150 (MQ=255)
                         gTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAAt                                                                                                                           <  2:3574492/150‑1 (MQ=255)
                           cAGGCTCCGACATCAGCTCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACAGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCAACCTGAATGa                                                                                                                         >  1:4363321/1‑150 (MQ=255)
                           cAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGa                                                                                                                         <  1:1735946/150‑1 (MQ=255)
                            aGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAt                                                                                                                        <  2:243982/150‑1 (MQ=255)
                             ggCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtc                                                                                                                       >  2:3189315/1‑150 (MQ=255)
                              gCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtct                                                                                                                      <  2:3029725/150‑1 (MQ=255)
                               cTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtctc                                                                                                                     <  1:4412788/150‑1 (MQ=255)
                               cTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtctc                                                                                                                     >  1:1313721/1‑150 (MQ=255)
                               cTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtctc                                                                                                                     >  2:2092227/1‑150 (MQ=255)
                               cTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGAtctc                                                                                                                     >  2:4436311/1‑150 (MQ=255)
                                  cGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGc                                                                                                                  >  1:2837440/1‑150 (MQ=255)
                                      aTCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTa                                                                                                              <  2:4181890/150‑1 (MQ=255)
                                         aGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTATAGCGTCGGTCGAGGAGTGCGTTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg                                                                                                           <  2:866613/150‑1 (MQ=255)
                                         aGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTg                                                                                                           <  2:3294245/150‑1 (MQ=255)
                                            gcCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGt                                                                                                        >  2:2622690/1‑150 (MQ=255)
                                            gcCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGt                                                                                                        >  1:3061771/1‑150 (MQ=255)
                                             ccTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTc                                                                                                       <  1:1491546/150‑1 (MQ=255)
                                             ccTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTc                                                                                                       >  1:3784771/1‑150 (MQ=255)
                                               tGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCg                                                                                                     >  1:2520523/1‑150 (MQ=255)
                                                gACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCgg                                                                                                    >  1:93428/1‑150 (MQ=255)
                                                gACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCgg                                                                                                    >  1:642332/1‑150 (MQ=255)
                                                   acaGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGgtgt                                                                                                 >  1:4617404/1‑150 (MQ=255)
                                                        agGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCt                                                                                            >  1:3465063/1‑150 (MQ=255)
                                                        agGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCt                                                                                            >  2:1394840/1‑150 (MQ=255)
                                                         gggTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCtt                                                                                           <  2:2720454/150‑1 (MQ=255)
                                                         gggTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCtt                                                                                           <  2:401082/150‑1 (MQ=255)
                                                          ggTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTc                                                                                          >  1:107667/1‑150 (MQ=255)
                                                            tttttCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGc                                                                                        >  1:4027554/1‑150 (MQ=255)
                                                              tttCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCa                                                                                      <  1:1341498/150‑1 (MQ=255)
                                                               ttCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAg                                                                                     >  1:889241/1‑150 (MQ=255)
                                                                tCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAgg                                                                                    >  2:551155/1‑150 (MQ=255)
                                                                tCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAgg                                                                                    <  1:3222683/150‑1 (MQ=255)
                                                                tCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAgg                                                                                    <  1:4067334/150‑1 (MQ=255)
                                                                 cgcgTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAggg                                                                                   >  1:1827555/1‑150 (MQ=255)
                                                                   cgTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGcc                                                                                 <  1:2629187/150‑1 (MQ=255)
                                                                    gTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCg                                                                                >  2:3424611/1‑150 (MQ=255)
                                                                       tcGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGAAGTGCGTCCGGTGACATCTTCGCCAGGGCCGAGc                                                                             >  1:2265498/1‑150 (MQ=255)
                                                                       tcGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGc                                                                             >  2:3545112/1‑150 (MQ=255)
                                                                            aGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATc                                                                        <  1:212233/150‑1 (MQ=255)
                                                                             gTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCg                                                                       >  2:3962706/1‑150 (MQ=255)
                                                                             gTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCg                                                                       >  2:3246178/1‑150 (MQ=255)
                                                                               ggAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgt                                                                     >  1:2139403/1‑150 (MQ=255)
                                                                                gAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgtg                                                                    <  1:3610497/150‑1 (MQ=255)
                                                                                 aaaGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGgtgt                                                                   <  1:2364977/150‑1 (MQ=255)
                                                                                  aaGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGtt                                                                  >  1:4068939/1‑150 (MQ=255)
                                                                                    ggTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTAt                                                                <  1:310136/150‑1 (MQ=255)
                                                                                       aTTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTcc                                                             >  2:3775580/1‑150 (MQ=255)
                                                                                           aaaGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTg                                                         >  2:2125031/1‑150 (MQ=255)
                                                                                           aaaGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTg                                                         <  2:1042552/150‑1 (MQ=255)
                                                                                               cGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCAt                                                     >  2:2239136/1‑150 (MQ=255)
                                                                                                 tcggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtc                                                   <  2:2524101/150‑1 (MQ=255)
                                                                                                   ggtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgt                                                 >  2:3719922/1‑150 (MQ=255)
                                                                                                    gtcgAGGAGGGCGTTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgtc                                                <  2:3082033/150‑1 (MQ=255)
                                                                                                    gtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgtc                                                <  2:2825170/150‑1 (MQ=255)
                                                                                                    gtcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgtc                                                <  1:682976/150‑1 (MQ=255)
                                                                                                     tcgAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATtcgtcg                                               <  1:1001959/150‑1 (MQ=255)
                                                                                                       gaggagGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCt                                             <  2:176805/150‑1 (MQ=255)
                                                                                                        aggaggGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTg                                            >  2:2769905/1‑150 (MQ=255)
                                                                                                        aggaggGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTg                                            >  2:1956056/1‑150 (MQ=255)
                                                                                                        aggaggGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTg                                            <  1:167583/150‑1 (MQ=255)
                                                                                                          gaggGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGAc                                          <  2:4555468/150‑1 (MQ=255)
                                                                                                             ggCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAg                                       <  2:93428/150‑1 (MQ=255)
                                                                                                             ggCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAg                                       <  2:642264/150‑1 (MQ=255)
                                                                                                                  tAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGttttt                                  >  1:788350/1‑150 (MQ=255)
                                                                                                                  tAATACCGAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGttttt                                  <  1:1945628/150‑1 (MQ=255)
                                                                                                                         gAACTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTNCCATTCGTCGCTGACGAGTTTTTGCAATTg                           <  1:4279782/150‑1 (MQ=255)
                                                                                                                          aaCTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGc                          <  1:1624597/150‑1 (MQ=255)
                                                                                                                           aCTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCt                         <  1:1053359/150‑1 (MQ=255)
                                                                                                                           aCTGCTGCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCt                         <  1:3732843/150‑1 (MQ=255)
                                                                                                                            ctgctgCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTc                        >  2:1300144/1‑150 (MQ=255)
                                                                                                                            ctgctgCGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTc                        <  1:3933852/150‑1 (MQ=255)
                                                                                                                             tgctgcGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCCTCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCg                       <  2:1122433/150‑1 (MQ=255)
                                                                                                                             tgctgcGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATATGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCg                       >  2:2053082/1‑150 (MQ=255)
                                                                                                                               ctgcGCCAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGaa                     >  1:145869/1‑150 (MQ=255)
                                                                                                                                    ccAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgaga                <  2:2916968/150‑1 (MQ=255)
                                                                                                                                     cAGTTGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGgagag               <  2:1619614/150‑1 (MQ=255)
                                                                                                                                        ttGCATGCAACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACg            <  1:4023252/150‑1 (MQ=255)
                                                                                                                                          gcatgcaACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAg          >  1:2966725/1‑150 (MQ=255)
                                                                                                                                               caACTCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGAtg     >  1:532428/1‑150 (MQ=255)
                                                                                                                                                  cgcgGTGAATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCTGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTgg  >  1:2226885/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGAACGGCTCATCAAGAATCAACAAGTCAGGCTCCGACATCAGCGCCTGACACAGCAGGGTTTTTCGCGTCTCGCCAGTGGAAAGGTATTTAAAGCGTCGGTCGAGGAGGGCGGTAATACCGAACTGCTGCGCCAGTTGCATGCAACGCGGTGCATCCTTTACTTCATCCTGAATGATCTCAGCCGTAGTGCGTCCGGTGTCATCTTCGCCAGGGCCGAGCATATCGGTGTTATTCCGCTGCCATTCGTCGCTGACGAGTTTTTGCAATTGCTCGAAGGAGAGACGAGTGATGTGG  >  CP009273/788755‑789050

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: