Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,033,430 G→T 14.8% D73Y (GAC→TAC)  fldB → flavodoxin 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,033,4300GT14.8% 924.6 / 111.7 303D73Y (GAC→TAC) fldBflavodoxin 2
Reads supporting (aligned to +/- strand):  ref base G (125/133);  new base T (15/30);  total (140/163)
Fisher's exact test for biased strand distribution p-value = 7.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 311 total aligned reads displayed.

ATATTATCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCAACGGA  >  CP009273/3033281‑3033578
                                                                                                                                                     |                                                                                                                                                    
atatTATCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTcg                                                                                                                                                      <  2:3773651/150‑1 (MQ=255)
   ttatCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTcgacg                                                                                                                                                   <  2:4812861/150‑1 (MQ=255)
   ttatCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTcgacg                                                                                                                                                   <  2:2905052/150‑1 (MQ=255)
       cGGCCCAGCACTGGTGACCTTACATAACCTCAAGGACGTCTCCCCGAAATTAATGGAGCAGTACGATGTGGTGATTCTGGGTATCCCGACCTGGGATTTGGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGTCGACCt                                                                                                                                               <  2:2004498/150‑1 (MQ=255)
          cccAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATACTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGaa                                                                                                                                            >  2:260218/1‑150 (MQ=255)
              gAACTTGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCtt                                                                                                                                        <  2:2724347/150‑1 (MQ=255)
              gAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCtt                                                                                                                                        <  2:2921978/150‑1 (MQ=255)
              gAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCtt                                                                                                                                        <  2:1129264/150‑1 (MQ=255)
               aaCTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTg                                                                                                                                       >  1:1474474/1‑150 (MQ=255)
                aCTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGa                                                                                                                                      >  2:777189/1‑150 (MQ=255)
                 cTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGaa                                                                                                                                     <  1:2376081/150‑1 (MQ=255)
                  tGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAg                                                                                                                                    <  2:3968443/150‑1 (MQ=255)
                    gTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGt                                                                                                                                  <  1:1240184/150‑1 (MQ=255)
                     tGACCTTACATAACCTCACGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTa                                                                                                                                 <  2:840976/150‑1 (MQ=255)
                        ccTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTaaaa                                                                                                                              <  1:343953/150‑1 (MQ=255)
                         cTTACATACCCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAAt                                                                                                                             <  2:1016207/150‑1 (MQ=255)
                          ttACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAAtt                                                                                                                            <  1:4435810/150‑1 (MQ=255)
                             cATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAAttgtt                                                                                                                         <  2:2113585/150‑1 (MQ=255)
                             cATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAAttgtt                                                                                                                         <  1:4064795/150‑1 (MQ=255)
                                  ccTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCt                                                                                                                    >  1:1040656/1‑150 (MQ=255)
                                   cTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTg                                                                                                                   <  1:3987776/150‑1 (MQ=255)
                                    tCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGt                                                                                                                  >  1:3117170/1‑150 (MQ=255)
                                        ggACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATgg                                                                                                              >  1:3375375/1‑150 (MQ=255)
                                        ggACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATgg                                                                                                              <  1:4044328/150‑1 (MQ=255)
                                        ggACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATgg                                                                                                              <  2:4594859/150‑1 (MQ=255)
                                        ggACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATgg                                                                                                              <  1:2896766/150‑1 (MQ=255)
                                         gacgacTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATggg                                                                                                             <  2:383409/150‑1 (MQ=255)
                                                  ccGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGAt                                                                                                    <  1:2257533/150‑1 (MQ=255)
                                                    gAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCa                                                                                                  <  1:3191645/150‑1 (MQ=255)
                                                    gAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCa                                                                                                  >  1:437883/1‑150 (MQ=255)
                                                     aaaTTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCaa                                                                                                 <  2:3125715/150‑1 (MQ=255)
                                                     aaaTTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCaa                                                                                                 <  2:2674024/150‑1 (MQ=255)
                                                       aTTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCaaat                                                                                               >  1:4166388/1‑148 (MQ=255)
                                                         tAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTgg                                                                                             >  2:654876/1‑150 (MQ=255)
                                                         tAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTgg                                                                                             >  2:3050844/1‑150 (MQ=255)
                                                          aaTGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTggg                                                                                            <  2:1736706/150‑1 (MQ=255)
                                                          aaTGGAGCAGTACGATGTGCTCATTCTGGGGATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGCCCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTggg                                                                                            <  2:2471412/150‑1 (MQ=255)
                                                             ggAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATa                                                                                         <  1:620143/150‑1 (MQ=255)
                                                               aGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGAATTTGGTGAAATACAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGAAGACCTGAACATTGAAGGTAAAATTGTTGCGATGTATGGGCTTGGCGATCAACTGGGATACg                                                                                       >  2:3604553/1‑150 (MQ=255)
                                                                 cAGTACGATGTGCTCATTCTGGGTATCCCAACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGc                                                                                     >  2:3265848/1‑150 (MQ=255)
                                                                 cAGTACGATGTGCTCATTCTGGGTATCCCAACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGc                                                                                     >  2:639028/1‑150 (MQ=255)
                                                                      cGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTg                                                                                >  2:3381383/1‑150 (MQ=255)
                                                                      cGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTg                                                                                <  2:4348650/150‑1 (MQ=255)
                                                                            gCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCt                                                                          <  2:4784610/150‑1 (MQ=255)
                                                                             cTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTc                                                                         >  1:2969599/1‑150 (MQ=255)
                                                                             cTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTc                                                                         <  1:1964487/150‑1 (MQ=255)
                                                                             cTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTc                                                                         <  1:1343002/150‑1 (MQ=255)
                                                                                aTTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGAt                                                                      <  2:4392550/150‑1 (MQ=255)
                                                                                     gggTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGAGTGGGAAGCCGTCTGGGATCAGCTGGACGACCTGAACCTTGAAGGTAAAATAGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGATCCTCGATGCGCt                                                                 >  2:1695416/1‑150 (MQ=255)
                                                                                     gggTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCt                                                                 >  1:2158711/1‑150 (MQ=255)
                                                                                     gggTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCt                                                                 <  2:2628270/150‑1 (MQ=255)
                                                                                      ggTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTc                                                                >  1:342136/1‑150 (MQ=255)
                                                                                      ggTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTc                                                                >  2:2226046/1‑150 (MQ=255)
                                                                                       gTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCg                                                               <  1:654876/150‑1 (MQ=255)
                                                                                          tCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTa                                                            <  1:2569402/150‑1 (MQ=255)
                                                                                           cccGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTAt                                                           <  1:192679/150‑1 (MQ=255)
                                                                                            ccGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTAtg                                                          <  1:510771/150‑1 (MQ=255)
                                                                                              gACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTAtgct                                                        >  1:4622076/1‑150 (MQ=255)
                                                                                              gACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTAtgct                                                        >  2:55526/1‑150 (MQ=255)
                                                                                               aCCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTAtgctg                                                       <  2:4236658/150‑1 (MQ=255)
                                                                                                ccTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTAtgctgc                                                      >  2:1289344/1‑150 (MQ=255)
                                                                                                 cTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCa                                                     >  1:2287714/1‑150 (MQ=255)
                                                                                                 cTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCa                                                     >  1:2407629/1‑150 (MQ=255)
                                                                                                   gggATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATg                                                   <  2:1931025/150‑1 (MQ=255)
                                                                                                   gggATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATg                                                   >  1:2101902/1‑150 (MQ=255)
                                                                                                    ggATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGa                                                  >  2:1803434/1‑150 (MQ=255)
                                                                                                         ttGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAAc                                             <  2:4250434/150‑1 (MQ=255)
                                                                                                           ggTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAActc                                           >  2:216612/1‑150 (MQ=255)
                                                                                                           ggTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAActc                                           >  2:4089531/1‑150 (MQ=255)
                                                                                                           ggTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAActc                                           >  1:1156774/1‑150 (MQ=255)
                                                                                                             tGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACtctc                                         >  2:1232001/1‑150 (MQ=255)
                                                                                                             tGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACtctc                                         >  2:415730/1‑150 (MQ=255)
                                                                                                               aaaTCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGa                                       <  1:1564877/150‑1 (MQ=255)
                                                                                                                aaTCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGAc                                      >  2:499488/1‑150 (MQ=255)
                                                                                                                 aTCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGAcc                                     <  1:3767654/150‑1 (MQ=255)
                                                                                                                 aTCCAGGAAGACTGGGAAGCCGTCAGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGAcc                                     >  2:955357/1‑150 (MQ=255)
                                                                                                                     aGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAg                                 <  1:3808408/150‑1 (MQ=255)
                                                                                                                     aGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAg                                 >  1:4679641/1‑150 (MQ=255)
                                                                                                                      ggAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAgg                                >  2:380143/1‑150 (MQ=255)
                                                                                                                        aaGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCg                              <  1:2725064/150‑1 (MQ=255)
                                                                                                                         aGACTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGt                             <  1:4998402/150‑1 (MQ=255)
                                                                                                                         aGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGt                             >  1:879253/1‑150 (MQ=255)
                                                                                                                            cTGGGAAGCCGTCTGGGATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGaa                          <  2:900040/150‑1 (MQ=255)
                                                                                                                              gggAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGt                        <  1:1272487/150‑1 (MQ=255)
                                                                                                                              gggAAGCCGTCTGGGATCAGCTCGAAGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGAATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGt                        >  2:3685751/1‑150 (MQ=255)
                                                                                                                                gAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTtc                      <  1:1453404/150‑1 (MQ=255)
                                                                                                                                 aaGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTtcg                     <  1:4547972/150‑1 (MQ=255)
                                                                                                                                 aaGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTtcg                     <  1:2301731/150‑1 (MQ=255)
                                                                                                                                 aaGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTtcg                     >  1:1917426/1‑150 (MQ=255)
                                                                                                                                  aGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTtcgt                    >  1:2818894/1‑150 (MQ=255)
                                                                                                                                    ccGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTtcgtcg                  <  2:1589790/150‑1 (MQ=255)
                                                                                                                                       tCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCt               <  2:3410617/150‑1 (MQ=255)
                                                                                                                                        cTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTa              <  1:785722/150‑1 (MQ=255)
                                                                                                                                          gggATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTGTCGACCAAAGGCGTGAAGTTCGTCGGCTACt            >  2:1372232/1‑150 (MQ=255)
                                                                                                                                           ggATCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTg           <  1:1773089/150‑1 (MQ=255)
                                                                                                                                            gATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTgg          >  2:2288180/1‑150 (MQ=255)
                                                                                                                                             aTCAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGc         >  1:2062042/1‑150 (MQ=255)
                                                                                                                                               cAGCTCTACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCa       <  2:1161172/150‑1 (MQ=255)
                                                                                                                                                   tCGACGACCTGAACCTTGAAGGTAAAATTGGTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCAACgg   <  1:1262579/150‑1 (MQ=255)
                                                                                                                                                    cgacgaCCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCAACGGa  <  2:5084129/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATATTATCGGCCCAGAACTGGTGACCTTACATAACCTCAAGGACGACTCCCCGAAATTAATGGAGCAGTACGATGTGCTCATTCTGGGTATCCCGACCTGGGATTTTGGTGAAATCCAGGAAGACTGGGAAGCCGTCTGGGATCAGCTCGACGACCTGAACCTTGAAGGTAAAATTGTTGCGCTGTATGGGCTTGGCGATCAACTGGGATACGGCGAGTGGTTCCTCGATGCGCTCGGTATGCTGCATGACAAACTCTCGACCAAAGGCGTGAAGTTCGTCGGCTACTGGCCAACGGA  >  CP009273/3033281‑3033578

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: