Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,009,187 G→T 88.7% intergenic (+90/‑232) mocA → / → ygfK CTP:molybdopterin cytidylyltransferase/putative Fe‑S subunit oxidoreductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,009,1870GT88.7% 674.6 / 83.2 254intergenic (+90/‑232)mocA/ygfKCTP:molybdopterin cytidylyltransferase/putative Fe‑S subunit oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base G (12/17);  new base T (119/106);  total (131/123)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

AGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGC  >  CP009273/3009040‑3009336
                                                                                                                                                   |                                                                                                                                                     
 gATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaaaagataaa                                                                                                                                                    <  1:2768041/150‑1 (MQ=255)
 gATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTgagaaaagataaa                                                                                                                                                    >  1:2616743/1‑150 (MQ=255)
    atTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAAt                                                                                                                                                 <  1:4548206/150‑1 (MQ=255)
      tGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTg                                                                                                                                               >  2:1157620/1‑150 (MQ=255)
      tGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTg                                                                                                                                               >  1:974851/1‑150 (MQ=255)
      tGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTg                                                                                                                                               >  1:4416375/1‑150 (MQ=255)
      tGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTg                                                                                                                                               <  1:513710/150‑1 (MQ=255)
        aTACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGAt                                                                                                                                             <  1:2247815/150‑1 (MQ=255)
         tACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATa                                                                                                                                            <  1:3465269/150‑1 (MQ=255)
          aCTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTATTTAAACGATTTCGTTCAAATATTTAATCCTTCCGGAAATGATTATCATAATGATAAGGAATAACTTCTTTTGAGAAAAGATAAAAATTGATaa                                                                                                                                           <  2:3516266/150‑1 (MQ=255)
          aCTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATaa                                                                                                                                           >  1:3046339/1‑150 (MQ=255)
           cTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATaaa                                                                                                                                          <  2:4385619/150‑1 (MQ=255)
           cTCCGGACGACTTTATTACAGCGAAGGAAAGGTAAACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATaaa                                                                                                                                          <  1:366938/150‑1 (MQ=255)
              cGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATgg                                                                                                                                       >  2:4508926/1‑150 (MQ=255)
                  cgacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTg                                                                                                                                   <  2:1391629/150‑1 (MQ=255)
                   gacTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTgg                                                                                                                                  >  2:2598346/1‑150 (MQ=255)
                     cTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGACACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCa                                                                                                                                >  2:3992811/1‑150 (MQ=255)
                         attaCAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCaataa                                                                                                                            >  2:4487447/1‑150 (MQ=255)
                          ttaCAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCaataat                                                                                                                           <  1:4412350/150‑1 (MQ=255)
                              aGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAAtt                                                                                                                       >  1:4236789/1‑150 (MQ=255)
                              aGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAAtt                                                                                                                       >  2:618596/1‑150 (MQ=255)
                              aGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAAtt                                                                                                                       >  1:3396513/1‑150 (MQ=255)
                                 gAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTc                                                                                                                    >  2:376724/1‑150 (MQ=255)
                                   aGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAg                                                                                                                  <  1:3164590/150‑1 (MQ=255)
                                     gAAAGGTATAATGAAATTTATAAAAATTATTTAAACGATTGCGTTCAAATATTTAATCGTTCATGAACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTAAgta                                                                                                                <  2:104381/150‑3 (MQ=255)
                                      aaaGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCCTTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAg                                                                                                               <  1:1723375/150‑1 (MQ=255)
                                       aaGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGt                                                                                                              >  1:3325261/1‑150 (MQ=255)
                                          gTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAt                                                                                                           >  1:4443206/1‑150 (MQ=255)
                                          gTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAt                                                                                                           >  1:556181/1‑150 (MQ=255)
                                          gTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAt                                                                                                           <  1:1465632/150‑1 (MQ=255)
                                             taCTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAtttt                                                                                                        <  1:1114434/150‑1 (MQ=255)
                                             taCTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAtttt                                                                                                        <  1:2082669/150‑1 (MQ=255)
                                              aCTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTa                                                                                                       >  2:1636666/1‑150 (MQ=255)
                                               cTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTaa                                                                                                      >  2:4073539/1‑150 (MQ=255)
                                                       tAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCAtt                                                                                              >  1:2104708/1‑150 (MQ=255)
                                                        aGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCAttt                                                                                             <  1:3755412/150‑1 (MQ=255)
                                                         gAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTa                                                                                            >  2:1730671/1‑150 (MQ=255)
                                                         gAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTa                                                                                            >  2:1729475/1‑150 (MQ=255)
                                                           aaaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGc                                                                                          <  2:3046339/150‑1 (MQ=255)
                                                           aaaCGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGc                                                                                          >  2:3892658/1‑150 (MQ=255)
                                                                  gTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCtgtg                                                                                   <  2:1118958/150‑1 (MQ=255)
                                                                   ttAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGa                                                                                  >  1:3797324/1‑150 (MQ=255)
                                                                   ttAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGAACCTGTGa                                                                                  >  1:4345781/1‑150 (MQ=255)
                                                                    tAAACGATTGCTTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAg                                                                                 <  2:2371073/150‑1 (MQ=255)
                                                                    tAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAg                                                                                 <  1:1325700/150‑1 (MQ=255)
                                                                      aaCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGtt                                                                               >  1:1893590/1‑150 (MQ=255)
                                                                      aaCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGtt                                                                               <  1:3516252/150‑1 (MQ=255)
                                                                      aaCGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGtt                                                                               <  2:2421543/150‑1 (MQ=255)
                                                                        cGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttca                                                                             >  2:2166597/1‑150 (MQ=255)
                                                                           ttGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGttcattc                                                                          <  2:3780091/150‑1 (MQ=255)
                                                                                 tCAAATATTTAATCCTTCCGGCACTCATTATCATAATTATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACaaaa                                                                    <  2:202429/150‑1 (MQ=255)
                                                                                       atTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCt                                                              >  1:3110442/1‑150 (MQ=255)
                                                                                        tttAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACATGTGAGTTCATTCACAAAAAGCCCtt                                                             <  2:295558/150‑1 (MQ=255)
                                                                                        tttAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACATGTGAGTTCATTCACAAAAAGCCCtt                                                             <  2:3414918/150‑1 (MQ=255)
                                                                                         ttAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTa                                                            <  1:1951583/150‑1 (MQ=255)
                                                                                           aaTCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAg                                                          <  2:1103419/150‑1 (MQ=255)
                                                                                            aTCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGt                                                         >  1:4626471/1‑150 (MQ=255)
                                                                                            aTCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGt                                                         <  2:4135438/150‑1 (MQ=255)
                                                                                             tCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGtt                                                        <  2:4536689/150‑1 (MQ=255)
                                                                                             tCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGtt                                                        >  2:1468138/1‑150 (MQ=255)
                                                                                             tCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATGTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGtt                                                        >  2:1931752/1‑150 (MQ=255)
                                                                                              ccTTCCGGCACTCATTATCATAATGATATGGAATCTCTGATTTTGAGTCAAGTTAATAATTTTTATATGGGCTGGCTTTAATAATTTTCATGAGTGATTTTAATTCTAATTTAGCACCTGTGAGTTCATTCACAAAAAGCGATTAAGttt                                                       <  2:846867/150‑1 (MQ=255)
                                                                                              ccTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGttt                                                       >  1:1169248/1‑150 (MQ=255)
                                                                                              ccTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGttt                                                       <  2:1552872/150‑1 (MQ=255)
                                                                                              ccTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGttt                                                       <  1:2034535/150‑1 (MQ=255)
                                                                                               cTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTc                                                      >  2:3616910/1‑150 (MQ=255)
                                                                                               cTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTc                                                      <  2:3942991/150‑1 (MQ=255)
                                                                                                 tCCGGCACTCATTATCATAATGATAAGGAATCACTTATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGa                                                    <  2:3285300/150‑1 (MQ=255)
                                                                                                 tCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGa                                                    >  2:4222576/1‑150 (MQ=255)
                                                                                                 tCCGGCACTCATTATCATAATGATAAGGAATAACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGa                                                    >  2:2479339/1‑150 (MQ=255)
                                                                                                  ccGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGAGTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAg                                                   >  1:4589125/1‑150 (MQ=255)
                                                                                                   cGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgc                                                  >  2:1231966/1‑150 (MQ=255)
                                                                                                      cACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAgccgc                                               >  1:2759187/1‑150 (MQ=255)
                                                                                                        cTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCt                                             <  2:2666542/150‑1 (MQ=255)
                                                                                                         tCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTc                                            <  2:2573404/150‑1 (MQ=255)
                                                                                                         tCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTc                                            <  2:4441359/150‑1 (MQ=255)
                                                                                                              aTCATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtt                                       <  2:463662/150‑1 (MQ=255)
                                                                                                               taaTAATGATAAGGATTCACTTATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAttt                                      <  1:2708199/148‑1 (MQ=255)
                                                                                                                cATAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAtttt                                     <  1:3644924/150‑1 (MQ=255)
                                                                                                                 aTAATGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAttttt                                    >  2:2796299/1‑150 (MQ=255)
                                                                                                                   aaTGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACAttttttt                                  <  1:1636666/150‑1 (MQ=255)
                                                                                                                     tGATAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTtat                                <  2:3644521/150‑1 (MQ=255)
                                                                                                                        tAAGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATAtt                             <  1:2587079/150‑1 (MQ=255)
                                                                                                                         aaGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATAttt                            >  1:2797026/1‑150 (MQ=255)
                                                                                                                          aGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTc                           >  2:31422/1‑150 (MQ=255)
                                                                                                                          aGGAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTc                           >  1:1823350/1‑150 (MQ=255)
                                                                                                                            gAATCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTccc                         >  2:2987500/1‑150 (MQ=255)
                                                                                                                             aaTCACTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCg                        <  1:2152882/150‑1 (MQ=255)
                                                                                                                                 aCTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCaa                    >  1:3876532/1‑150 (MQ=255)
                                                                                                                                  cTGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCaaa                   >  1:2498805/1‑150 (MQ=255)
                                                                                                                                   tGATTTTGAGAAAAGATAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACATGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAAc                  >  2:2331622/1‑150 (MQ=255)
                                                                                                                                         tGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCa            <  2:3056903/150‑1 (MQ=255)
                                                                                                                                         tGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCa            >  2:2499483/1‑150 (MQ=255)
                                                                                                                                          gagaaaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCaa           >  1:3778878/1‑150 (MQ=255)
                                                                                                                                            gaaaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAaga         >  2:553953/1‑150 (MQ=255)
                                                                                                                                             aaaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAagag        >  1:4167534/1‑150 (MQ=255)
                                                                                                                                              aaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAgt       <  2:974851/150‑1 (MQ=255)
                                                                                                                                              aaagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAgt       >  1:4110939/1‑150 (MQ=255)
                                                                                                                                               aagataaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAgtg      <  1:1220248/150‑1 (MQ=255)
                                                                                                                                                   gaaaaATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGc  >  1:2726531/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGATATTGATACTCCGGACGACTTTATTACAGCGAAGGAAAGGTATACTGAAATTTAGAAAACGTAGTTAAACGATTGCGTTCAAATATTTAATCCTTCCGGCACTCATTATCATAATGATAAGGAATCACTGATTTTGAGAAAAGAGAAAAATTGATAAATGGGCTGGCAATAATAATTTTCAGGAGTGATTTTAATTCTCATTTAGCACCTGTGAGTTCATTCACAAAAAGCCCTTAAGTTTCGAGCCGCCTCACATTTTTTTATATTTCCCGCCAAACCTGGCAAGAGTGGTGC  >  CP009273/3009040‑3009336

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: