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.5% 1798.8 / 0.4 372intergenic (+11/+140)PA14_RS25260/PA14_RS25265tRNA‑Met/HigA family addiction module antitoxin
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base . (189/178);  minor base C (0/2);  total (191/181)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 411 total aligned reads displayed.

GAAAATTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATCC  >  NC_008463/5513907‑5514201
                                                                                                                                                     |                                                                                                                                                 
gAAAATTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAg                                                                                                                                                   >  2:1365599/1‑150 (MQ=255)
  aaaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggg                                                                                                                                                 >  2:3627/1‑150 (MQ=255)
   aaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAgggg                                                                                                                                                >  2:8795605/1‑150 (MQ=255)
    aTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCACGTCCCGGTGGAGCCACCATATTTTTCAAggggt                                                                                                                                               >  2:5171089/1‑149 (MQ=255)
    aTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggt                                                                                                                                               <  1:2400029/150‑2 (MQ=255)
    aTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggt                                                                                                                                               <  1:6468263/150‑2 (MQ=255)
    aTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggt                                                                                                                                               >  2:6685646/1‑149 (MQ=255)
     ttATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggtt                                                                                                                                              <  2:628345/150‑3 (MQ=255)
     ttATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggtt                                                                                                                                              <  1:5595958/150‑3 (MQ=255)
       aTCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggttag                                                                                                                                            >  1:5991734/1‑146 (MQ=255)
        tCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAgc                                                                                                                                          <  1:9786582/150‑1 (MQ=255)
          tAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAgcgc                                                                                                                                        >  2:2896927/1‑150 (MQ=255)
          tAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAgcgc                                                                                                                                        >  1:5012442/1‑150 (MQ=255)
           aGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCa                                                                                                                                       >  2:8531555/1‑150 (MQ=255)
             gCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAg                                                                                                                                     >  2:2246008/1‑150 (MQ=255)
              cccAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGc                                                                                                                                    <  1:9894310/150‑1 (MQ=255)
               ccAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCt                                                                                                                                   <  1:3976176/150‑1 (MQ=255)
                 aGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTGCAAGTCCCGGTGTAGCCACCAGAGTTTTCATGGGGTTAGCGCAAGCTat                                                                                                                                 >  1:6296714/1‑149 (MQ=255)
                  ggggTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAAc                                                                                                                                <  2:4399168/150‑1 (MQ=255)
                      gtGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCt                                                                                                                            <  1:4784762/149‑1 (MQ=255)
                       tGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCtt                                                                                                                           <  2:8010895/150‑1 (MQ=255)
                        gACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCTAGGGGTTAGCGCAAGCTAACCCCttt                                                                                                                          >  2:3119844/1‑150 (MQ=255)
                         aCAGCCCCGCCCGAGCTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCGAACCCCtttt                                                                                                                         <  2:5043713/150‑1 (MQ=255)
                         aCAGCCCCGCCCGAGATACGGAAACTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCGAACCCCtttt                                                                                                                         <  2:6842505/150‑1 (MQ=255)
                         aCAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCtttt                                                                                                                         >  2:9522171/1‑150 (MQ=255)
                           aGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttg                                                                                                                       <  1:8529168/150‑1 (MQ=255)
                            gCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgt                                                                                                                      >  1:7278291/1‑150 (MQ=255)
                             ccccGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgtt                                                                                                                     <  1:273218/150‑1 (MQ=255)
                               ccgcccgAGATACGGAAACTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCAttttgtttt                                                                                                                   <  2:9877538/150‑1 (MQ=255)
                               ccgcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgtttt                                                                                                                   <  1:10800150/150‑1 (MQ=255)
                               ccgcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgtttt                                                                                                                   <  2:5367055/150‑1 (MQ=255)
                               ccgcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgtttt                                                                                                                   <  2:8286107/150‑1 (MQ=255)
                                cgcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgttttg                                                                                                                  <  1:6083788/150‑1 (MQ=255)
                                cgcccgAGATACGGAAAATGCGCGCCACACGGCTACAGAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCAttttgttttg                                                                                                                  <  1:690722/150‑1 (MQ=255)
                                  cccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGt                                                                                                                >  1:3102293/1‑150 (MQ=255)
                                  cccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCATGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGt                                                                                                                >  2:9253431/1‑150 (MQ=255)
                                   ccgAGATACGGAAAATGCTCTCCACTCGGCTACATAGCTCAGTCGGTTAGAGGGATGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTATCGCTTGATAACCCCTTTTTTTTTTGGTc                                                                                                               <  2:5212897/150‑1 (MQ=255)
                                   ccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTc                                                                                                               <  2:478822/150‑1 (MQ=255)
                                   ccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTc                                                                                                               >  2:6479253/1‑150 (MQ=255)
                                    cgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCg                                                                                                              <  1:8742714/150‑1 (MQ=255)
                                     gagaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCgg                                                                                                             <  1:9244475/150‑1 (MQ=255)
                                      agaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGc                                                                                                            >  1:10003939/1‑150 (MQ=255)
                                      agaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGc                                                                                                            >  1:1241043/1‑150 (MQ=255)
                                       gaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCa                                                                                                           >  2:2391442/1‑150 (MQ=255)
                                        aTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATACTGCTGATGTCCCAGGTTCACGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCATTTTGTTTTGGTCGGCAg                                                                                                          <  2:3980457/150‑1 (MQ=255)
                                                       gcCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgacta                                                                                           >  1:595653/1‑150 (MQ=255)
                                                       gcCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgacta                                                                                           >  2:552510/1‑150 (MQ=255)
                                                        ccACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactac                                                                                          <  1:10092415/150‑1 (MQ=255)
                                                         cACTCGGCTACATAGCTCAGTCGGATAGAGCGCAGCATTCATAATGCTGCTGTCCCAGGTTCACGTCCCGCTGTCGCCACCAGATGTTTCAAGGGGTTAGCGCAAGCTAACCCCATTTTGTTTTGGTCGGCAgactacaaaccgactaca                                                                                         <  2:1074153/150‑1 (MQ=255)
                                                          aCTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactacaa                                                                                        >  1:10531431/1‑150 (MQ=255)
                                                           cTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactacaaa                                                                                       <  2:10440456/150‑1 (MQ=255)
                                                               gCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGt                                                                                   >  2:7140320/1‑150 (MQ=255)
                                                               gCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGt                                                                                   <  1:5210642/150‑1 (MQ=255)
                                                                cTACATAGCTCAGTCGTTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTATTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTc                                                                                  <  1:2683051/150‑1 (MQ=255)
                                                                 tACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCt                                                                                 <  1:7472225/150‑1 (MQ=255)
                                                                  aCATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTa                                                                                <  1:11038065/150‑1 (MQ=255)
                                                                  aCATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTa                                                                                <  1:8659866/150‑1 (MQ=255)
                                                                        cTCAGTCGGTTAGAGCGCAGCATTAATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTACAAGGGGTTAGCGCAAGCTAACCCCTATTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGcacc                                                                          >  2:1116195/1‑150 (MQ=255)
                                                                            gTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCttt                                                                      <  1:3512029/150‑1 (MQ=255)
                                                                             tCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCATTa                                                                     >  1:5318160/1‑150 (MQ=255)
                                                                                 ttAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGa                                                                 <  1:429063/150‑1 (MQ=255)
                                                                                     agCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTGTTGGTCGGCAGACTACACACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTa                                                             >  2:6933651/1‑150 (MQ=255)
                                                                                      gcgcAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTAc                                                            <  1:2088013/150‑1 (MQ=255)
                                                                                       cgcAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACg                                                           >  2:8891195/1‑150 (MQ=255)
                                                                                          agcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGtt                                                        <  2:2300094/150‑1 (MQ=255)
                                                                                          agcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCGTTACCGATTTACGGtt                                                        >  1:1669050/1‑150 (MQ=255)
                                                                                           gcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAGGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTc                                                       >  2:8481896/1‑150 (MQ=255)
                                                                                             aTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGc                                                     >  1:2962959/1‑150 (MQ=255)
                                                                                             aTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGc                                                     >  1:7042435/1‑150 (MQ=255)
                                                                                               tCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCAt                                                   <  2:587343/150‑1 (MQ=255)
                                                                                                 aTAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGa                                                 >  2:4213697/1‑150 (MQ=255)
                                                                                                  tAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGaa                                                <  1:973621/150‑1 (MQ=255)
                                                                                                     tGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGaaaaa                                             >  2:7557621/1‑150 (MQ=255)
                                                                                                     tGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGaaaaa                                             >  1:4521266/1‑150 (MQ=255)
                                                                                                     tGCAGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGaaaaa                                             <  1:6663316/150‑1 (MQ=255)
                                                                                                      gCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAg                                            >  1:10223990/1‑150 (MQ=255)
                                                                                                       cTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGc                                           <  1:4607955/150‑1 (MQ=255)
                                                                                                         gATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGccc                                         >  2:143834/1‑150 (MQ=255)
                                                                                                          aTGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGcccc                                        <  1:841422/150‑1 (MQ=255)
                                                                                                          aTGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGcccc                                        <  2:10219629/150‑1 (MQ=255)
                                                                                                          aTGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGcccc                                        >  1:1816832/1‑150 (MQ=255)
                                                                                                           tGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCg                                       >  2:7307087/1‑150 (MQ=255)
                                                                                                            gTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGc                                      <  2:193931/150‑1 (MQ=255)
                                                                                                              cccAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCAt                                    >  1:7917580/1‑150 (MQ=255)
                                                                                                                       aaGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGc                           >  1:8222481/1‑150 (MQ=255)
                                                                                                                          tCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGctct                        <  1:763957/150‑1 (MQ=255)
                                                                                                                           cccGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTa                       <  1:9015547/150‑1 (MQ=255)
                                                                                                                            ccGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCTCATGAAAAAGCCCCGCATCTGCGGGGCTCTAc                      >  2:10307919/1‑150 (MQ=255)
                                                                                                                              ggTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGTTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGa                    <  1:1524784/150‑1 (MQ=255)
                                                                                                                                tgtAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATc                  >  1:5780271/1‑150 (MQ=255)
                                                                                                                                 gtAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCa                 >  2:7086728/1‑150 (MQ=255)
                                                                                                                                   aGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAAc               >  1:9358272/1‑150 (MQ=255)
                                                                                                                                     ccaccaTATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGt             >  2:2070252/1‑150 (MQ=255)
                                                                                                                                        ccaTATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCAc          <  2:4717310/150‑1 (MQ=255)
                                                                                                                                           tatTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGc       <  2:9448654/150‑1 (MQ=255)
                                                                                                                                           tatTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGc       >  1:5436102/1‑150 (MQ=255)
                                                                                                                                               tttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATc   >  2:5761582/1‑150 (MQ=255)
                                                                                                                                               tttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACAGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCAACTGAGGGGCTCTACGATCAACGTAACCGCTATc   >  1:10448513/1‑150 (MQ=255)
                                                                                                                                                ttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATcc  <  2:3414464/150‑1 (MQ=255)
                                                                                                                                                ttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTGTGGTCGGCAGACTACAAACCGTCTACAATCCGTCTGCGCCAGCTTTACCGATTTACGGTTCGCTTGACAAAGCCCCGCTTCTGCGGGGCTCTACGATCACCGTCACCGCTATcc  >  2:6125430/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: