Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 828,374 G→T 61.7% V11V (GTC→GTA)  ybiA ← DUF1768 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273828,3740GT61.7% 133.0 / 275.2 201V11V (GTC→GTA) ybiADUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base G (38/39);  new base T (64/60);  total (102/99)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 202 total aligned reads displayed.

TTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTTC  >  CP009273/828227‑828523
                                                                                                                                                   |                                                                                                                                                     
tttttCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTAc                                                                                                                                                     >  2:1145595/1‑150 (MQ=255)
    tcatcaAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACgtgt                                                                                                                                                 >  2:1218001/1‑150 (MQ=255)
      atcaAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTg                                                                                                                                               >  1:812139/1‑150 (MQ=255)
      atcaAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTg                                                                                                                                               >  2:2638394/1‑150 (MQ=255)
        caAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGa                                                                                                                                             >  1:3062995/1‑150 (MQ=255)
         aaGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGAt                                                                                                                                            >  2:3718930/1‑150 (MQ=255)
          aGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATc                                                                                                                                           <  2:4091426/150‑1 (MQ=255)
             aaTTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATccgt                                                                                                                                        <  1:1007053/150‑3 (MQ=255)
             aaTTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCtt                                                                                                                                        >  1:3092481/1‑150 (MQ=255)
              aTTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCttt                                                                                                                                       <  2:358163/150‑1 (MQ=255)
               tttttGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTtg                                                                                                                                      <  1:2592264/150‑1 (MQ=255)
                   tgtgCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCt                                                                                                                                  <  2:3677424/150‑1 (MQ=255)
                   tgtgCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCt                                                                                                                                  >  2:169317/1‑150 (MQ=255)
                   tgtgCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCt                                                                                                                                  <  1:4808663/150‑1 (MQ=255)
                    gtgCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTc                                                                                                                                 <  1:3509105/150‑1 (MQ=255)
                    gtgCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTc                                                                                                                                 >  2:3973321/1‑150 (MQ=255)
                        cTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAAc                                                                                                                             >  2:2912460/1‑150 (MQ=255)
                         tGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACg                                                                                                                            <  2:1214048/150‑1 (MQ=255)
                           aaaaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACggg                                                                                                                          >  2:2493934/1‑150 (MQ=255)
                           aaaaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACggg                                                                                                                          >  2:4326611/1‑150 (MQ=255)
                             aaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCa                                                                                                                        >  2:2677446/1‑150 (MQ=255)
                             aaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCa                                                                                                                        >  1:2391960/1‑150 (MQ=255)
                             aaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCa                                                                                                                        >  2:2677477/1‑150 (MQ=255)
                             aaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCa                                                                                                                        >  2:4250388/1‑150 (MQ=255)
                             aaTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCa                                                                                                                        <  1:495269/150‑1 (MQ=255)
                              aTAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCAt                                                                                                                       >  2:2339938/1‑150 (MQ=255)
                                  tgtTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTaaa                                                                                                                   >  1:2498659/1‑150 (MQ=255)
                                  tgtTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTaaa                                                                                                                   >  1:479069/1‑150 (MQ=255)
                                  tgtTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTaaa                                                                                                                   <  1:540942/150‑1 (MQ=255)
                                       tGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCa                                                                                                              <  2:2384243/150‑1 (MQ=255)
                                       tGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCa                                                                                                              >  1:2728326/1‑150 (MQ=255)
                                         aGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCAtt                                                                                                            >  1:2170212/1‑150 (MQ=255)
                                          ggTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTa                                                                                                           <  1:2048394/150‑1 (MQ=255)
                                             aGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGt                                                                                                        >  1:1632240/1‑150 (MQ=255)
                                             aGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGt                                                                                                        >  1:1143544/1‑150 (MQ=255)
                                                 cAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggt                                                                                                    <  2:323864/150‑1 (MQ=255)
                                                  aGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTtggtg                                                                                                   >  2:708636/1‑150 (MQ=255)
                                                     ttttCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGtt                                                                                                >  1:611054/1‑150 (MQ=255)
                                                      tttCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTa                                                                                               <  2:3386754/150‑1 (MQ=255)
                                                       ttCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAt                                                                                              <  1:3821236/150‑1 (MQ=255)
                                                        tCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAtt                                                                                             >  1:4620092/1‑150 (MQ=255)
                                                        tCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTAtt                                                                                             <  1:2262034/150‑1 (MQ=255)
                                                         cccGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTg                                                                                            <  2:3800919/150‑1 (MQ=255)
                                                            gTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGa                                                                                         >  1:4598664/1‑150 (MQ=255)
                                                            gTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGa                                                                                         >  2:2615656/1‑150 (MQ=255)
                                                                   ttGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCtt                                                                                  <  1:1050947/150‑1 (MQ=255)
                                                                   ttGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCtt                                                                                  <  2:1477870/150‑1 (MQ=255)
                                                                     gATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTGTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTc                                                                                <  2:4641459/150‑1 (MQ=255)
                                                                      aTGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTcc                                                                               <  2:1356716/150‑1 (MQ=255)
                                                                       tGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCca                                                                              >  1:2655819/1‑150 (MQ=255)
                                                                            cAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaac                                                                         <  1:3082834/150‑1 (MQ=255)
                                                                              ggcggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacag                                                                       >  2:2001785/1‑150 (MQ=255)
                                                                               gcggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacaga                                                                      <  1:2770426/150‑1 (MQ=255)
                                                                                 ggcAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCcagaaacagaaa                                                                    >  1:4180002/1‑150 (MQ=255)
                                                                                      aGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCa                                                               >  1:3031086/1‑150 (MQ=255)
                                                                                       gTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCAt                                                              <  1:729051/150‑1 (MQ=255)
                                                                                       gTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCAt                                                              <  2:2735023/150‑1 (MQ=255)
                                                                                       gTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCAt                                                              <  1:1348024/150‑1 (MQ=255)
                                                                                        ttAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCAtt                                                             <  1:4736073/150‑1 (MQ=255)
                                                                                        ttAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCAtt                                                             >  1:4017206/1‑150 (MQ=255)
                                                                                        ttAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCAtt                                                             <  1:2351139/150‑1 (MQ=255)
                                                                                           gAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTAcc                                                          >  1:1602478/1‑150 (MQ=255)
                                                                                           gAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTAcc                                                          <  1:3777127/150‑1 (MQ=255)
                                                                                            aaaaaTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTAccc                                                         >  2:1317298/1‑150 (MQ=255)
                                                                                            aaaaaTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTAccc                                                         <  1:159361/150‑1 (MQ=255)
                                                                                              aaaTCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCt                                                       >  2:1503365/1‑150 (MQ=255)
                                                                                                   aCCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGaaaa                                                  <  2:2552751/150‑1 (MQ=255)
                                                                                                    ccGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAAc                                                 >  2:2181199/1‑150 (MQ=255)
                                                                                                    ccGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAAc                                                 >  2:3037222/1‑150 (MQ=255)
                                                                                                     cGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAAcc                                                <  2:2755807/150‑1 (MQ=255)
                                                                                                      gTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCg                                               <  2:1217241/150‑1 (MQ=255)
                                                                                                       tAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGa                                              <  2:3166057/150‑1 (MQ=255)
                                                                                                          tCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGaaaa                                           >  1:444096/1‑150 (MQ=255)
                                                                                                            aTCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAAt                                         <  1:1365007/150‑1 (MQ=255)
                                                                                                            aTCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAAt                                         <  2:2757258/150‑1 (MQ=255)
                                                                                                             tCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATg                                        <  1:2941697/150‑1 (MQ=255)
                                                                                                             tCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATg                                        <  2:2385258/150‑1 (MQ=255)
                                                                                                                cTGGTGCTGTAGACATATATGATGGTGTCTTGCAGTACGGGTTGGATCCAGTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCca                                     <  2:4382720/150‑1 (MQ=255)
                                                                                                                cTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCca                                     <  2:699830/150‑1 (MQ=255)
                                                                                                                cTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCca                                     <  1:3712592/150‑1 (MQ=255)
                                                                                                                 tGGTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCcac                                    >  1:1190076/1‑150 (MQ=255)
                                                                                                                  ggTGCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCcaca                                   <  1:1307967/150‑1 (MQ=255)
                                                                                                                     gCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAata                                <  1:3326211/150‑1 (MQ=255)
                                                                                                                     gCTGTAGAAATTTATGATGGTGTCTTGCATTACGTGTTGGATCCTGTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAata                                <  1:2800164/150‑1 (MQ=255)
                                                                                                                       tGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATAtt                              <  1:367521/150‑1 (MQ=255)
                                                                                                                            aaaTTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTg                         >  1:3852227/1‑150 (MQ=255)
                                                                                                                            aaaTTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTg                         >  1:899640/1‑150 (MQ=255)
                                                                                                                                 tATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTtat                    >  2:3941144/1‑150 (MQ=255)
                                                                                                                                 tATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTtat                    <  2:574017/150‑1 (MQ=255)
                                                                                                                                 tATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTtat                    >  2:3396990/1‑150 (MQ=255)
                                                                                                                                  atgatgGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTtata                   <  2:2787734/150‑1 (MQ=255)
                                                                                                                                     atgGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAg                >  2:2465418/1‑150 (MQ=255)
                                                                                                                                      tgGTGTCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGt               >  1:2283071/1‑150 (MQ=255)
                                                                                                                                        gtgtCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAt             >  1:2361372/1‑150 (MQ=255)
                                                                                                                                        gtgtCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAt             <  1:1800161/150‑1 (MQ=255)
                                                                                                                                        gtgtCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAt             >  1:1299466/1‑150 (MQ=255)
                                                                                                                                         tgtCTTGCATTACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTAtt            >  1:3628996/1‑150 (MQ=255)
                                                                                                                                            cTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCa         >  1:4254714/1‑150 (MQ=255)
                                                                                                                                                 aTGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGttt    >  2:2926831/1‑150 (MQ=255)
                                                                                                                                                 aTGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGttt    >  1:3259962/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTTTTCATCAAGGAATTTTTGTGCCTGAAAATAGTGTTCTGAGGTAGGCCAGGTTTTCCCGTCAACCTTGATGGGCCAGGCGGCAAAGTTAGAAAAATCACCGTAGTCATCGCTGGTGCTGTAGAAATTTATGATGGTGTCTTGCATGACGTGTTGGATCCTTTGTGCTCGAACGGGCATTAAACCGCATTATGTTGGTGGTTATTGCGAGCCGCTTTCCAGAAACAGAAAAACCATTACCCCTGAAAACCGAAAAATGCCACAATATTGGCTGTTTATACAGTATTTCAGGTTTTC  >  CP009273/828227‑828523

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: ‑