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.98.9% 1782.7 / 11.1 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 . (170/196);  minor base T (4/0);  total (174/198)
Fisher's exact test for biased strand distribution p-value = 4.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 400 total aligned reads displayed.

GAAAATTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTATCC  >  NC_008463/5513907‑5514201
                                                                                                                                                     |                                                                                                                                                 
 aaaaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAgg                                                                                                                                                  >  2:8004111/1‑150 (MQ=255)
  aaaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggg                                                                                                                                                 >  2:162603/1‑150 (MQ=255)
   aaTTATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCATGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAgggg                                                                                                                                                <  1:4518918/150‑1 (MQ=255)
     ttATCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggtt                                                                                                                                              >  2:8527295/1‑148 (MQ=255)
       aTCTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAggggttag                                                                                                                                            >  2:1809142/1‑146 (MQ=255)
         cTAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAgcg                                                                                                                                         >  2:7514395/1‑150 (MQ=255)
          tAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTTGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCATGGGGTTAgcgc                                                                                                                                        >  2:215084/1‑150 (MQ=255)
          tAGGCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCGGATGTCCCAGGTTCAAGTCCCGGTGTCGCCACCATATTTTTCAAGGGGTTAgcgc                                                                                                                                        >  2:2738120/1‑150 (MQ=255)
             gCCCAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAg                                                                                                                                     >  2:1314067/1‑150 (MQ=255)
               ccAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCt                                                                                                                                   <  2:7750693/150‑1 (MQ=255)
               ccAGGGGTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCt                                                                                                                                   <  2:1970016/150‑1 (MQ=255)
                 aGGGGTTGACAGCCCCGCCCGATATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGAAGCATGCATAATTCAGATGTCCCAGGTTGAAGTCCAGGTGTAGCAACAATATTTTTCTAGGGGTTAGCGCTAGCTaa                                                                                                                                 >  1:1884326/1‑150 (MQ=255)
                  ggggTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAAc                                                                                                                                <  2:5502972/150‑1 (MQ=255)
                  ggggTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAAc                                                                                                                                >  2:3829147/1‑150 (MQ=255)
                  ggggTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAAc                                                                                                                                >  1:6185053/1‑150 (MQ=255)
                   gggTTGACAGCCCCGCCCGAGATACGGAAATTGCGCGCCACTAGGCTACATATCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAAcc                                                                                                                               <  2:1481781/150‑1 (MQ=255)
                     gTTGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAAcccc                                                                                                                             <  1:10251643/150‑1 (MQ=255)
                       tGACAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCtt                                                                                                                           >  1:1001781/1‑150 (MQ=255)
                          cAGCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCttttt                                                                                                                        <  2:7079431/150‑1 (MQ=255)
                            gCCCCGCCCGAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgt                                                                                                                      >  2:1719741/1‑150 (MQ=255)
                              cccgcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgttt                                                                                                                    <  1:135432/150‑1 (MQ=255)
                               ccgcccgAGATACGGAAAATGCGCGCCACGCGGCGACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGACCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCCAGCTAACCCCTttttgtttt                                                                                                                   <  2:2603135/150‑1 (MQ=255)
                                cgcccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTttttgttttg                                                                                                                  <  1:9799823/150‑1 (MQ=255)
                                  cccgAGATACGTAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGt                                                                                                                <  2:2930677/150‑1 (MQ=255)
                                  cccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGt                                                                                                                >  2:3570048/1‑150 (MQ=255)
                                   ccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTc                                                                                                               >  2:3722050/1‑150 (MQ=255)
                                   ccgAGATACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTc                                                                                                               >  2:507260/1‑150 (MQ=255)
                                      agaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGc                                                                                                            <  2:1911402/150‑1 (MQ=255)
                                      agaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGc                                                                                                            >  2:10492550/1‑150 (MQ=255)
                                      agaTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGc                                                                                                            >  1:1668632/1‑150 (MQ=255)
                                        aTACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTCACCCCTTTTTGTTTTGGTCGGCAg                                                                                                          >  1:2802620/1‑150 (MQ=255)
                                         tACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAga                                                                                                         <  1:3399095/150‑1 (MQ=255)
                                         tACGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAga                                                                                                         <  1:1075571/150‑1 (MQ=255)
                                           cGGAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgact                                                                                                       >  1:6807976/1‑150 (MQ=255)
                                            ggAAAATGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgacta                                                                                                      >  1:3257146/1‑150 (MQ=255)
                                                aaTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaa                                                                                                  <  1:2124395/150‑1 (MQ=255)
                                                 aTGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaac                                                                                                 <  1:9271096/150‑1 (MQ=255)
                                                  tGCGCGCCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaacc                                                                                                >  2:7248370/1‑150 (MQ=255)
                                                     gcgcCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgac                                                                                             <  1:7522027/150‑1 (MQ=255)
                                                       gcCACTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgacta                                                                                           <  1:3133588/150‑1 (MQ=255)
                                                          aCTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactacaa                                                                                        >  2:4899633/1‑150 (MQ=255)
                                                          aCTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactacaa                                                                                        <  2:6270492/150‑1 (MQ=255)
                                                           cTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactacaaa                                                                                       <  1:3706820/150‑1 (MQ=255)
                                                           cTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactacaaa                                                                                       >  2:9217811/1‑150 (MQ=255)
                                                           cTCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactacaaa                                                                                       >  1:6439072/1‑150 (MQ=255)
                                                            tCGGCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAgactacaaaccgactacaaac                                                                                      <  1:8747789/150‑1 (MQ=255)
                                                              ggCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGactacaaaccgactacaaaccg                                                                                    >  1:1253383/1‑150 (MQ=255)
                                                              ggCTACAGAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCAGATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGactacaaaccgactacaaaccg                                                                                    >  2:9324203/1‑150 (MQ=255)
                                                               gCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGt                                                                                   <  1:2898972/150‑1 (MQ=255)
                                                               gCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGt                                                                                   <  2:8283764/150‑1 (MQ=255)
                                                               gCTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGt                                                                                   >  1:5012607/1‑150 (MQ=255)
                                                                cTACATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTc                                                                                  >  1:1552650/1‑150 (MQ=255)
                                                                  aCATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTa                                                                                >  2:6426393/1‑150 (MQ=255)
                                                                   cATAGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTAc                                                                               <  1:7473839/150‑1 (MQ=255)
                                                                      aGCTCAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGcc                                                                            >  1:2461941/1‑150 (MQ=255)
                                                                          cAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCt                                                                        <  2:7398197/150‑1 (MQ=255)
                                                                          cAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCt                                                                        <  1:3431699/150‑1 (MQ=255)
                                                                          cAGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCt                                                                        <  2:9781416/150‑1 (MQ=255)
                                                                           aGTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCttt                                                                      >  2:10193352/1‑150 (MQ=255)
                                                                            gTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCttt                                                                      <  1:4409968/150‑1 (MQ=255)
                                                                            gTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCttt                                                                      <  1:3538715/150‑1 (MQ=255)
                                                                            gTCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCATTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCttt                                                                      <  2:6679196/150‑1 (MQ=255)
                                                                             tCGGTTAGAGCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTa                                                                     <  1:8092994/150‑1 (MQ=255)
                                                                                   agagCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAtt                                                               >  2:4553148/1‑150 (MQ=255)
                                                                                   agagCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAtt                                                               <  2:9931899/150‑1 (MQ=255)
                                                                                   agagCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAtt                                                               <  2:1229485/150‑1 (MQ=255)
                                                                                    gagCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAttt                                                              <  1:5543289/150‑1 (MQ=255)
                                                                                     agCGCAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTATTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCAACGCCCCCTTTACCGATTTa                                                             >  1:387543/1‑150 (MQ=255)
                                                                                      gcgcAGCATTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCTTATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTAc                                                            <  2:9622549/150‑1 (MQ=255)
                                                                                         cagcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGt                                                         <  2:4436853/150‑1 (MQ=255)
                                                                                          agcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGtt                                                        >  2:1124100/1‑150 (MQ=255)
                                                                                           gcaTTCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTc                                                       >  1:520070/1‑150 (MQ=255)
                                                                                              ttCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCa                                                    >  1:10538657/1‑150 (MQ=255)
                                                                                               tCATAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACCCCCCCTTTACCGATTTACGGTTCGCAt                                                   <  1:8932278/150‑1 (MQ=255)
                                                                                                 aTAATGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGAGTTACGGTTCGCATGa                                                 >  1:8299671/1‑150 (MQ=255)
                                                                                                   aaTGCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGaaa                                               <  1:7704003/150‑1 (MQ=255)
                                                                                                      gCTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAg                                            >  1:7392607/1‑150 (MQ=255)
                                                                                                       cTGATGTCCCTGGTTCAAGTCCCGGTGTAGCCACCTTATTTTTCAATGGGTTATCGCAATCTAACACCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTTCGCACCATTTACCGATTTTCGGTTCGCATGTAAAAGc                                           <  2:8785370/150‑1 (MQ=255)
                                                                                                       cTGATGTCCCAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGc                                           <  1:854721/150‑1 (MQ=255)
                                                                                                              cccAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCAt                                    <  1:7211382/150‑1 (MQ=255)
                                                                                                              cccAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCAt                                    <  1:8045357/150‑1 (MQ=255)
                                                                                                               ccAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATc                                   >  2:4041173/1‑150 (MQ=255)
                                                                                                                cAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCt                                  >  2:3427534/1‑150 (MQ=255)
                                                                                                                cAGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCATTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCt                                  <  2:5846237/150‑1 (MQ=255)
                                                                                                                 aGGTTCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTg                                 >  1:5780870/1‑150 (MQ=255)
                                                                                                                    ttCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTCCCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCgg                              <  2:5292708/150‑1 (MQ=255)
                                                                                                                    ttCAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCgg                              <  2:3120141/150‑1 (MQ=255)
                                                                                                                      cAAGTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCgggg                            >  1:6277044/1‑150 (MQ=255)
                                                                                                                         gTCCCGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGctc                         >  2:5659418/1‑150 (MQ=255)
                                                                                                                           cccTGTGTAGCCACCATATTTTTCAAGGGTTTAGCTCAAGCTAACCCCATATTTTTTTGGTCGGCTGACTACAAACCGAATACACACCGTCTCCTCAACCTTTAACGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTa                       <  2:860299/150‑1 (MQ=255)
                                                                                                                            ccGGTGTAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTAc                      >  1:9825617/1‑150 (MQ=255)
                                                                                                                               gtgtAGCCACCATATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGAt                   >  2:2349087/1‑150 (MQ=255)
                                                                                                                                      caccaTATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTc            <  1:4253551/150‑1 (MQ=255)
                                                                                                                                        ccaTATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCAc          <  1:9897078/150‑1 (MQ=255)
                                                                                                                                        ccaTATTTTTCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCAc          <  1:6060401/150‑1 (MQ=255)
                                                                                                                                            ctttttCATGGGGTTAGCGCACGCTACCCCCTTTTTGTTGTGGTCGGCCGACTACAACCCGACTACACACCGTCTCCGCCCCCTTTACCGAGTTACGGGGCGCATGAAAAAGCCCCGCAGCTGCGGGGCGCTACGAGCAACGTCACCGCt      >  2:5910243/2‑150 (MQ=255)
                                                                                                                                             tttttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTa     >  1:1061654/1‑150 (MQ=255)
                                                                                                                                             tttttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTa     >  1:956475/1‑150 (MQ=255)
                                                                                                                                              ttttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTATACCGATTTACGGTTCGCATGAAAAAGCCCCGCATCTGCGGGGCTCTACGATCAACGTCACCGCTAt    >  2:2078374/1‑150 (MQ=255)
                                                                                                                                                ttCAAGGGGTTAGCGCAAGCTAACCCCTTTTTGTTTTGGTCGGCAGACTACAAACCGACTACAAACCGTCTACGCCCCCTTTACCGATTTACGGTTCGCATGAAAAAGCACCGCATCTGCGGGGCTCTAAGATCAACGTCACCGCTATcc  >  2:9438021/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: