Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,957,549 2 bp→AT intergenic (‑371/‑300) PA14_RS22700 ← / → PA14_RS22705 16S ribosomal RNA/exonuclease SbcCD subunit D C‑terminal domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,957,5490GA100.0% 475.9 / NA 237intergenic (‑371/‑301)PA14_RS22700/PA14_RS2270516S ribosomal RNA/exonuclease SbcCD subunit D C‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (103/134);  total (103/134)
*NC_0084634,957,5500CT99.6% 471.9 / NA 237intergenic (‑372/‑300)PA14_RS22700/PA14_RS2270516S ribosomal RNA/exonuclease SbcCD subunit D C‑terminal domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (103/133);  minor base G (0/1);  total (103/134)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

TGCGCCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGATGAAGCTTCGTTTCGAGCCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACT  >  NC_008463/4957455‑4957693
                                                                                              ||                                                                                                                                               
 gcgcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTc                                                                                          <  2:3410011/150‑1 (MQ=9)
 gcgcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTc                                                                                          >  1:7063698/1‑150 (MQ=9)
   gcCGACCTTCACCCCATCTCCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAAGTCGCTGTCACGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCt                                                                                        <  1:10255983/150‑1 (MQ=11)
   gcCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCt                                                                                        <  2:8531083/150‑1 (MQ=12)
    ccGACCTTCAAACCATCTTCCGATCATCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTg                                                                                       >  2:103616/1‑150 (MQ=255)
    ccGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTg                                                                                       <  1:5738520/150‑1 (MQ=12)
    ccGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTg                                                                                       <  2:1213328/150‑1 (MQ=12)
     cGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTgg                                                                                      <  2:3945312/150‑1 (MQ=12)
     cGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTgg                                                                                      <  1:5240686/150‑1 (MQ=12)
     cGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTgg                                                                                      <  1:6764829/150‑1 (MQ=12)
      gACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACTACCTGCCGGAGCGCTTCACCACTCCTggg                                                                                     <  1:7292705/150‑1 (MQ=255)
      gACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTggg                                                                                     >  2:1009493/1‑150 (MQ=12)
       aCCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGga                                                                                    <  2:4745786/150‑1 (MQ=12)
          ttCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagag                                                                                 <  1:4713425/150‑1 (MQ=12)
          ttCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagag                                                                                 <  1:4328987/150‑1 (MQ=12)
           tCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGgagaga                                                                                >  2:1420170/1‑150 (MQ=12)
              aaCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATc                                                                             >  1:4549622/1‑150 (MQ=12)
              aaCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATc                                                                             >  2:6160478/1‑150 (MQ=12)
              aaCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATc                                                                             >  1:9757351/1‑150 (MQ=12)
               aCCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATcc                                                                            >  2:4151708/1‑150 (MQ=12)
                 cATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGt                                                                          <  2:925155/150‑1 (MQ=12)
                  aTCTTCCGATCAGCGGGAGGCGCATTCTACAGAGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCTGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTa                                                                         >  2:389572/1‑150 (MQ=11)
                    cTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTaaa                                                                       <  1:2228558/150‑1 (MQ=12)
                      tCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATc                                                                     <  1:9361278/150‑1 (MQ=12)
                      tCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATc                                                                     <  2:949204/150‑1 (MQ=12)
                         gATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATACGTAAATCCtt                                                                  >  2:2205163/1‑150 (MQ=255)
                          aTCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTg                                                                 >  2:3431361/1‑150 (MQ=12)
                           tCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGCACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGa                                                                <  2:6568611/150‑1 (MQ=255)
                            cAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAt                                                               >  1:1857168/1‑150 (MQ=12)
                             aGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAtt                                                              <  2:5108981/150‑1 (MQ=12)
                              gCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAttt                                                             >  2:9165033/1‑150 (MQ=12)
                                gggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTc                                                           <  1:10450704/150‑1 (MQ=12)
                                gggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTc                                                           <  2:10634500/150‑1 (MQ=12)
                                gggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTc                                                           >  2:7222977/1‑150 (MQ=12)
                                gggAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTc                                                           >  2:3092239/1‑150 (MQ=12)
                                  gAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCaa                                                         >  2:7510831/1‑150 (MQ=12)
                                  gAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCaa                                                         >  1:107244/1‑150 (MQ=12)
                                    ggCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGc                                                       <  1:7070008/150‑1 (MQ=12)
                                     gcgcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAAACCGTAAATCCTTGATTTTCAAGCa                                                      >  2:964969/1‑150 (MQ=255)
                                      cgcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAt                                                     >  1:11040121/1‑150 (MQ=12)
                                       gcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAtt                                                    >  2:7766438/1‑150 (MQ=12)
                                       gcATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAtt                                                    >  2:9917177/1‑150 (MQ=12)
                                        cATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCAttt                                                   >  2:5610720/1‑150 (MQ=12)
                                        cATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAAGCCTTGATTTTCAAGCAttt                                                   >  1:9448114/1‑150 (MQ=255)
                                         aTTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTa                                                  >  2:4322621/1‑150 (MQ=12)
                                          ttCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTAc                                                 <  1:794278/150‑1 (MQ=12)
                                           tCTACAGCGTTCACATTCGCTGACAATCACCTCCGTGAATCTTCGTATCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGAATTTCAAGCATTTAct                                                <  2:1791537/150‑1 (MQ=11)
                                           tCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTAct                                                <  1:5724849/150‑1 (MQ=12)
                                            cTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTActc                                               <  1:9815555/150‑1 (MQ=12)
                                            cTACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTActc                                               <  1:9820422/150‑1 (MQ=12)
                                             tACAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTActct                                              <  1:3203577/150‑1 (MQ=12)
                                             tACAGCGGTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTActct                                              <  2:9051380/150‑1 (MQ=255)
                                               cAGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtc                                            <  2:6788682/150‑1 (MQ=12)
                                                aGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcg                                           <  1:8214836/150‑1 (MQ=12)
                                                aGCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcg                                           <  2:716414/150‑1 (MQ=12)
                                                aGCGTTCAAATTCGCTGTCAAGCACCACGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAACTCCTTGATTTTCAAGCATTTACTCTtcg                                           <  2:4225004/150‑1 (MQ=11)
                                                 gCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgt                                          >  1:2681628/1‑150 (MQ=12)
                                                 gCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgt                                          >  2:600870/1‑150 (MQ=12)
                                                 gCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgt                                          <  1:2976810/150‑1 (MQ=12)
                                                 gCGTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgt                                          >  2:1979717/1‑150 (MQ=12)
                                                   gTTCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTAGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcg                                        >  1:10050335/1‑150 (MQ=255)
                                                    ttCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgt                                       <  1:1748470/150‑1 (MQ=12)
                                                    ttCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgt                                       >  1:411737/1‑150 (MQ=12)
                                                    ttCAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgt                                       <  2:9113543/150‑1 (MQ=12)
                                                      cAAATTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTtcgtcgtcg                                     >  1:3763243/1‑150 (MQ=12)
                                                       aaaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGc                                    <  1:9834952/150‑1 (MQ=12)
                                                        aaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGTATCCGTAATTCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGcc                                   >  2:5397987/1‑150 (MQ=11)
                                                        aaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGcc                                   >  1:5974582/1‑150 (MQ=12)
                                                        aaTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGcc                                   >  1:1348457/1‑150 (MQ=12)
                                                         aTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCg                                  <  1:4930662/150‑1 (MQ=12)
                                                         aTTCGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCg                                  >  2:754285/1‑150 (MQ=12)
                                                            cGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGa                               <  2:6957626/150‑1 (MQ=12)
                                                            cGCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGa                               >  1:7253849/1‑150 (MQ=12)
                                                             gCTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAg                              <  1:10649771/150‑1 (MQ=12)
                                                              cTGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgt                             >  1:3589490/1‑150 (MQ=12)
                                                               tGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtg                            >  1:2480040/1‑150 (MQ=12)
                                                               tGTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtg                            >  1:10836201/1‑150 (MQ=12)
                                                                gTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtgg                           <  1:4633220/150‑1 (MQ=12)
                                                                gTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtgg                           >  2:10815774/1‑150 (MQ=12)
                                                                gTCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtgg                           >  1:6406646/1‑150 (MQ=12)
                                                                 tCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtggt                          <  2:4595529/150‑1 (MQ=12)
                                                                 tCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtggt                          >  1:3703860/1‑150 (MQ=12)
                                                                 tCAAGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAgtggt                          <  1:9008982/150‑1 (MQ=12)
                                                                   aaGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTgc                        >  1:8658928/1‑150 (MQ=12)
                                                                   aaGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTgc                        <  1:7274117/150‑1 (MQ=12)
                                                                    aGCACCTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTgcg                       <  2:2666566/150‑1 (MQ=12)
                                                                        ccTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCagt                   >  2:9366656/1‑148 (MQ=12)
                                                                        ccTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCAtt                   >  2:3543665/1‑150 (MQ=12)
                                                                         cTCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATta                  <  2:2775183/150‑1 (MQ=12)
                                                                          tCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATtat                 <  1:11006750/150‑1 (MQ=12)
                                                                          tCGGTGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATtat                 <  2:811009/150‑1 (MQ=12)
                                                                            ggtGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAg               >  1:1192444/3‑150 (MQ=12)
                                                                             gtGAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAgg              <  1:4875938/149‑1 (MQ=17)
                                                                               gAAGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGga            >  2:10261759/1‑150 (MQ=17)
                                                                                aaGCTTCGTTTCGAATCCCTTTCGGGACCGCACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgag           <  2:6550475/150‑1 (MQ=12)
                                                                                aaGCTTCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgag           <  1:8377154/150‑1 (MQ=17)
                                                                                 aGCTTCGTTACGAATCCCTTTCTGGACCGAACCGCTGCGTCGAACAACCTGCCGGAGCTCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGgaga          <  1:1385954/150‑1 (MQ=11)
                                                                                    ttCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAg       <  1:4192866/150‑1 (MQ=17)
                                                                                     tCGTTTCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGa      <  2:3989551/150‑1 (MQ=17)
                                                                                        tttCGAATCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAAc   <  1:1785597/150‑1 (MQ=17)
                                                                                              ||                                                                                                                                               
TGCGCCGACCTTCAAACCATCTTCCGATCAGCGGGAGGCGCATTCTACAGCGTTCAAATTCGCTGTCAAGCACCTCGATGAAGCTTCGTTTCGAGCCCCTTTCGGGACCGAACCGCTGCGGCGAACAACCTGCCGGAGCGCTTCACCACTCCTGGGAGAGAATCCGTAAATCCTTGATTTTCAAGCATTTACTCTTCGTCGTCGCCGGGAGTGGTGCGCATTATAGGGAGATAGAAACT  >  NC_008463/4957455‑4957693

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: