Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,514,060 (G)5→4 intergenic (+15/+136) PA14_RS25260 → / ← PA14_RS25265 tRNA‑Met/HigA family addiction module antitoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084635,514,0560G.99.4% 1692.3 / 1.8 346intergenic (+11/+140)PA14_RS25260/PA14_RS25265tRNA‑Met/HigA family addiction module antitoxin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (185/159);  minor base C (0/2);  total (185/161)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 376 total aligned reads displayed.

GAAAATTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATCC  >  NC_008463/5513907‑5514201
                                                                                                                                                     |                                                                                                                                                 
gAAAATTCTCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAg                                                                                                                                                   <  2:1874177/150‑1 (MQ=255)
gAAAATTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAg                                                                                                                                                   <  1:10642280/150‑1 (MQ=255)
 aaaaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAgg                                                                                                                                                  <  1:3021395/150‑1 (MQ=255)
   aaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGACAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAgggg                                                                                                                                                >  1:5165028/1‑150 (MQ=255)
   aaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAgggg                                                                                                                                                >  2:3484303/1‑150 (MQ=255)
   aaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAgggg                                                                                                                                                <  2:733043/150‑1 (MQ=255)
   aaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAgggg                                                                                                                                                >  2:3851040/1‑150 (MQ=255)
      tATCGAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggtta                                                                                                                                             <  1:10395517/150‑4 (MQ=255)
        tCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAgc                                                                                                                                          <  1:6715151/150‑1 (MQ=255)
            ggCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCaa                                                                                                                                      <  2:1063193/150‑1 (MQ=255)
            ggCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCaa                                                                                                                                      <  1:8936444/150‑1 (MQ=255)
            ggCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCaa                                                                                                                                      <  1:2810741/150‑1 (MQ=255)
             gCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAg                                                                                                                                     >  2:9123107/1‑150 (MQ=255)
             gCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAg                                                                                                                                     <  1:8639023/150‑1 (MQ=255)
              cccAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGc                                                                                                                                    >  2:3706323/1‑150 (MQ=255)
                  ggggTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAAc                                                                                                                                >  2:5986958/1‑150 (MQ=255)
                  ggggTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGGTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCTAGCTAAc                                                                                                                                >  2:839342/1‑150 (MQ=255)
                     gTTGACAGCCCCGCCCGAGATACGGAAAATGGGCGCCACTCGGCTACATAGCTCAGTCGGTTAGCGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAAccgc                                                                                                                             >  2:312041/1‑148 (MQ=255)
                         aCAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCtttt                                                                                                                         >  2:9396080/1‑150 (MQ=255)
                          cAGCCCCGCCCGAGATACGGAACATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCttttt                                                                                                                        <  1:6598767/150‑1 (MQ=255)
                          cAGCCCCGCCCGAGATACGGAAACTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCttttt                                                                                                                        <  2:5165281/150‑1 (MQ=255)
                          cAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCttttt                                                                                                                        <  2:1502414/150‑1 (MQ=255)
                          cAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCttttt                                                                                                                        >  2:8795894/1‑150 (MQ=255)
                          cAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCttttt                                                                                                                        <  1:5861182/150‑1 (MQ=255)
                           aGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttg                                                                                                                       >  1:4000696/1‑150 (MQ=255)
                            gCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgt                                                                                                                      <  1:74803/150‑1 (MQ=255)
                              cccgcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgttt                                                                                                                    <  1:7410550/150‑1 (MQ=255)
                               ccgcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgtttt                                                                                                                   <  2:7826706/150‑1 (MQ=255)
                                cgcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgttttg                                                                                                                  <  2:7511571/150‑1 (MQ=255)
                                 gcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTgg                                                                                                                 >  1:8280180/1‑150 (MQ=255)
                                     gagaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCgg                                                                                                             >  1:10756684/1‑150 (MQ=255)
                                     gagaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCgg                                                                                                             <  2:2429550/150‑1 (MQ=255)
                                         tACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAga                                                                                                         >  1:1943639/1‑150 (MQ=255)
                                         tACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAga                                                                                                         >  2:4610035/1‑150 (MQ=255)
                                         tACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTGTTTGTGTTGGTCGGCAga                                                                                                         >  1:4222493/1‑150 (MQ=255)
                                          aCGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgac                                                                                                        <  2:7626346/150‑1 (MQ=255)
                                            ggAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGAGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgacta                                                                                                      >  2:6839029/1‑150 (MQ=255)
                                             gAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactac                                                                                                     >  2:10759879/1‑150 (MQ=255)
                                              aaaaTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCCAGCTAACCCCGTTTTGTTTTGGTCGGCAgactaca                                                                                                    >  1:4828342/1‑150 (MQ=255)
                                              aaaaTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGCGTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTATTTGTTTTGGTCGGCAgactaca                                                                                                    <  2:8231197/150‑1 (MQ=255)
                                                aaTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaa                                                                                                  >  1:367207/1‑150 (MQ=255)
                                                  tGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCGAACCCCTTTTTGTTTTGGTCGGCAgactacaaacc                                                                                                >  2:6972071/1‑150 (MQ=255)
                                                           cTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactacaaa                                                                                       <  1:9329792/150‑1 (MQ=255)
                                                            tCGGCTACATAGCTCAGTCTGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgtctacaaaccgactacaaac                                                                                      <  2:7186886/150‑1 (MQ=255)
                                                                cTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTc                                                                                  >  1:6412395/1‑150 (MQ=255)
                                                                   cATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGTTCGGCAGACTACAAACCGACTACAAACCGTCTAc                                                                               >  2:5119111/1‑150 (MQ=255)
                                                                   cATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTAc                                                                               >  2:10522426/1‑150 (MQ=255)
                                                                    aTAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCACGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACg                                                                              <  2:1296994/150‑1 (MQ=255)
                                                                    aTAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAATCCGTCTACg                                                                              >  1:1563381/1‑150 (MQ=255)
                                                                     tAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGc                                                                             <  1:7192487/150‑1 (MQ=255)
                                                                        cTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGcccc                                                                          >  2:8029962/1‑150 (MQ=255)
                                                                          cAGTCGGTTAGAGCTCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCt                                                                        <  1:8147009/150‑1 (MQ=255)
                                                                           aGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCtt                                                                       >  2:1977355/1‑150 (MQ=255)
                                                                           aGTCGGTTAGAGCGAAGCTTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGACACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTTCAAACCGTCTACGCCCCCtt                                                                       <  2:4028174/150‑1 (MQ=255)
                                                                            gTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCttt                                                                      >  1:9092771/1‑150 (MQ=255)
                                                                            gTCGGTTAGAGCGCAGCATTCATAATGCTGAAGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTACCCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCttt                                                                      >  2:7871046/1‑150 (MQ=255)
                                                                             tCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTa                                                                     <  2:9046264/150‑1 (MQ=255)
                                                                              cGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTAc                                                                    >  2:10412449/1‑150 (MQ=255)
                                                                                gTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCg                                                                  >  2:9191658/1‑150 (MQ=255)
                                                                                  tAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAt                                                                <  1:2281596/150‑1 (MQ=255)
                                                                                  tAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAt                                                                <  1:5661847/150‑1 (MQ=255)
                                                                                   agagCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAtt                                                               <  2:9808571/150‑1 (MQ=255)
                                                                                   agagCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAtt                                                               >  1:4630438/1‑150 (MQ=255)
                                                                                   agagCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAtt                                                               >  1:7889689/1‑150 (MQ=255)
                                                                                   agagCGCAGCATTCATAATGCTGAGGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAtt                                                               <  2:5426660/150‑1 (MQ=255)
                                                                                    gagCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTCACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGCCTACAAACCGTCTACGCCCCCTTTACCGAgtt                                                              >  1:9528254/1‑150 (MQ=255)
                                                                                     agCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTa                                                             >  2:8554493/1‑150 (MQ=255)
                                                                                        gaagcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGATAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACgg                                                          <  2:5850864/148‑1 (MQ=255)
                                                                                                 aTAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGa                                                 >  2:3755753/1‑150 (MQ=255)
                                                                                                   aaTGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCTCCCTTTACCGATTTACGGTTCGCATGaaa                                               <  1:3840402/150‑1 (MQ=255)
                                                                                                           tGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCg                                       >  1:6075186/1‑150 (MQ=255)
                                                                                                            gTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGc                                      <  1:2661076/150‑1 (MQ=255)
                                                                                                             tCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCa                                     <  2:7487436/150‑1 (MQ=255)
                                                                                                              cccAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCAt                                    <  2:4891874/150‑1 (MQ=255)
                                                                                                              cccAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCAt                                    >  1:7004819/1‑150 (MQ=255)
                                                                                                               ccAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATc                                   <  2:8633776/150‑1 (MQ=255)
                                                                                                                 aGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTg                                 >  2:4185060/1‑150 (MQ=255)
                                                                                                                  ggTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGc                                >  1:5283162/1‑150 (MQ=255)
                                                                                                                   gTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCg                               <  1:10181642/150‑1 (MQ=255)
                                                                                                                    ttCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCgg                              <  2:6926249/150‑1 (MQ=255)
                                                                                                                      caacccccGCTGTAGCCACCATATTTTTCACCGGGTTAGCGCAAGCTAACCCCTTTGTGTTTGCGTCGGCAGCCTACAAACCGCCTACCCACCGTCTACGCCCCCGGTACCGACTGACGGTTCGCATGAAAAAGCCCCGCATCTGCgggg                            <  2:2760263/145‑1 (MQ=255)
                                                                                                                      cAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCgggg                            <  1:5472965/150‑1 (MQ=255)
                                                                                                                      cAAGTCCCGGTGTAGCCACCATAGTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCgggg                            >  2:2993097/1‑150 (MQ=255)
                                                                                                                         gTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGctc                         >  1:6163541/1‑150 (MQ=255)
                                                                                                                         gTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGctc                         <  1:555313/150‑1 (MQ=255)
                                                                                                                         gTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGctc                         >  2:6444298/1‑150 (MQ=255)
                                                                                                                         gTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACAACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGctc                         >  2:253280/1‑150 (MQ=255)
                                                                                                                             cGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACg                     >  2:1773935/1‑150 (MQ=255)
                                                                                                                               gtgtAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACACACCGTCTACGCCCCCTTTACCGAATTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGAt                   <  2:1003348/150‑1 (MQ=255)
                                                                                                                                  tAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCaa                >  1:1650894/1‑150 (MQ=255)
                                                                                                                                    gCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGctc                         <  1:3783973/139‑1 (MQ=255)
                                                                                                                                       accaTATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCa           <  1:2515312/150‑1 (MQ=255)
                                                                                                                                       accaTATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCa           >  1:10073776/1‑150 (MQ=255)
                                                                                                                                        ccaTATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCAc          <  1:4582557/150‑1 (MQ=255)
                                                                                                                                        ccaTATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCAc          <  1:8534709/150‑1 (MQ=255)
                                                                                                                                            atTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCt      >  2:5921942/1‑150 (MQ=255)
                                                                                                                                              ttttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTAt    >  1:9728228/1‑150 (MQ=255)
                                                                                                                                               tttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATc   >  1:3183904/1‑150 (MQ=255)
                                                                                                                                                ttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATcc  <  1:7054963/150‑1 (MQ=255)
                                                                                                                                                ttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATcc  >  1:3390967/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GAAAATTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATCC  >  NC_008463/5513907‑5514201

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: