Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,515,863 T→C intergenic (+617/+267) PA14_RS16005 → / ← PA14_RS16010 AAA family ATPase/DUF4354 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,515,8630TC99.0% 1132.9 / NA 302intergenic (+617/+267)PA14_RS16005/PA14_RS16010AAA family ATPase/DUF4354 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (137/162);  minor base G (2/1);  total (139/163)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 302 total aligned reads displayed.

GTATACAGCCCCTCCCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCTTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAA  >  NC_008463/3515715‑3516010
                                                                                                                                                     |                                                                                                                                                   
 tataCAGCCCCTCCCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCa                                                                                                                                                  <  2:8621890/151‑1 (MQ=255)
    aCAGCCCCTCCCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCAccc                                                                                                                                               >  1:7646760/1‑151 (MQ=255)
    aCAGCCCCTCCCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCAccc                                                                                                                                               <  1:11141598/151‑1 (MQ=255)
     cAGCCCCTCCCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCAcccc                                                                                                                                              <  1:4067693/151‑1 (MQ=255)
        ccccTCCCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCACCCCAAc                                                                                                                                           >  1:6865895/1‑151 (MQ=255)
         ccctccctAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAAcc                                                                                                                                          <  2:6455663/151‑1 (MQ=255)
         ccctccctAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCACACCCCaaac                                                                                                                                          >  2:4860240/1‑149 (MQ=255)
            tccctAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCctc                                                                                                                                       <  1:8417255/151‑1 (MQ=255)
            tccctAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCACACCCCAAcccac                                                                                                                                       >  2:1453470/1‑149 (MQ=255)
             ccctAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCctct                                                                                                                                      <  2:2495157/151‑1 (MQ=255)
              cctAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAACCCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCtctc                                                                                                                                     <  2:9658118/151‑1 (MQ=255)
                tAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCCCAACCCAAaccacaacc                                                                                                                                   >  1:2937018/1‑145 (MQ=255)
                 aaGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCg                                                                                                                                  <  2:6096222/151‑1 (MQ=255)
                 aaGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCACACCCCAACCCACAcccc                                                                                                                                  >  1:3498069/1‑150 (MQ=255)
                  aGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGc                                                                                                                                 <  2:10843876/151‑1 (MQ=255)
                    ccATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCAc                                                                                                                               <  1:2747264/151‑1 (MQ=255)
                    ccATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCAc                                                                                                                               <  2:6757273/151‑1 (MQ=255)
                    ccATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGCACCCCCACCCACTACCGCAc                                                                                                                               >  1:5560366/1‑151 (MQ=255)
                         gCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCACACCCCCACCCACACCCCCCCCCggg                                                                                                                          >  1:7150095/1‑151 (MQ=255)
                            cTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCACCCCAAACCACTCCCCCACCCGGGggag                                                                                                                      >  1:5947911/1‑151 (MQ=255)
                             ttCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGagag                                                                                                                      <  2:2734783/151‑1 (MQ=255)
                              tCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCACCCCAAACCACTCCCCCACGCGGGAGAgg                                                                                                                     >  2:92329/1‑151 (MQ=255)
                                cccATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCACACCCCAAACCAAACCCCCACCCGGGGGAgggg                                                                                                                   >  1:5188987/1‑151 (MQ=255)
                                 ccATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCACCCCAAACCACACCCCCAAGCGGGGGAggggg                                                                                                                  >  1:308070/1‑151 (MQ=255)
                                 ccATACCTGAACGGCACACCAATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGACCCCCAACCCGCACCCCCCCGCCGGAGAggggg                                                                                                                  >  2:6523610/1‑151 (MQ=255)
                                  cATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGCACCCCAACCCTCTCCCCCAAGCGGGAGTGGGGGc                                                                                                                 >  2:11158774/1‑151 (MQ=255)
                                     aCCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCgg                                                                                                              <  2:1009702/151‑1 (MQ=255)
                                     aCCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAAACCTCTCCCGCACGAGGGAGAGGGGGCCgg                                                                                                              >  2:8016605/1‑151 (MQ=255)
                                      ccTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCACACCCCAACCCACACCCGCACGCGGGAGAGGGGGCCggg                                                                                                             >  2:1656618/1‑151 (MQ=255)
                                        tGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCg                                                                                                           >  2:5967603/1‑151 (MQ=255)
                                        tGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGCACCCCAAACCACACCCCCCCCCCGGGGTGGGGGGCCGGCg                                                                                                           >  1:10339088/1‑151 (MQ=255)
                                                    cTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACGGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGGTTTTGGAATGTGGTGTGGGGCCGCACCCCATCCCGCTCCCGCACGGGGGGGAGGGGGCGGGGGGGGGCGTGTgggg                                                                                               >  2:5126735/1‑150 (MQ=255)
                                                     tATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCACTCCCCCACGCCGGAGAGGGGGCCCGGCGTGCCGTGTGGGtg                                                                                              >  1:7558956/1‑150 (MQ=255)
                                                       ttGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCGTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCACTCCCCCCCGCGGGTGAGGGGGGCGGGCGGGCCAGGTGGGTAGa                                                                                            >  2:5310208/1‑151 (MQ=255)
                                                        tGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAg                                                                                           <  2:7839119/151‑1 (MQ=255)
                                                         gCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCCGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAgc                                                                                          <  2:7475291/151‑1 (MQ=255)
                                                          ccACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACCCGGGGGGGGGGGCCGGGCGTGGCGTGTGGGTCGTgcg                                                                                         >  2:2431342/1‑151 (MQ=255)
                                                          ccACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGCACCCCAACCCACTCCCCCACCCGGGAGAGGGGGCCGGGCGTGGCGTGTGGGTCGAgcg                                                                                         >  1:10096323/1‑151 (MQ=255)
                                                           cACCCGCAAAACCACACCCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAgcgc                                                                                        <  1:3893127/151‑1 (MQ=255)
                                                                gCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCACTCCCGCACCCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCctcgc                                                                                   >  2:7125016/1‑151 (MQ=255)
                                                                  aaaaCCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCACAACCCAACCCTCACCCCCCCCCCGGGGAGGGGGGCCGGCGTGCCGTGTGGGTGGAGCGCcgcgcgc                                                                                 >  1:3648124/1‑149 (MQ=255)
                                                                   acaCCACAACCTCCCCCGCGGCTAACCGCAATAGCTGCGGCTGGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGCACCCCAACCCCCTCCCGCCCGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCctcgctcg                                                                                <  1:607297/149‑1 (MQ=255)
                                                                   aaaCCACACCCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCctcgctcg                                                                                <  1:6795873/151‑1 (MQ=255)
                                                                    aaCCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCgg                                                                               <  2:10383640/151‑1 (MQ=255)
                                                                     aCCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGCACCCCCACCCACACCCCCCCCCCGGGGAGGGGGGCCGGCGTGGCGTGGTGGTCGAGCCCCACCCGCggg                                                                              >  1:1813342/1‑151 (MQ=255)
                                                                     aCCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCACAACCCAAACCAAAAACCCACCCGGGGGAGGGGGCCGGGCGTGCCGTGGTGGTGGAGGCCCACGCACggg                                                                              >  1:1590399/1‑151 (MQ=255)
                                                                      ccACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGgc                                                                             <  1:6276995/151‑1 (MQ=255)
                                                                       cacaACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGCACCCCAACCCGCTCCCGCACGCGGGCGTGGGGGGCGGGGGGGCCGGGGGGGGCGGGCGCCGCGCGGGGgcg                                                                            >  2:7016036/1‑151 (MQ=255)
                                                                        acacCCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGgcgc                                                                           <  1:855432/151‑1 (MQ=255)
                                                                         caACCTCCCCCGCTGCTAACCGCAATCGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCa                                                                          <  1:1306254/151‑1 (MQ=255)
                                                                         caACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCa                                                                          <  2:974712/151‑1 (MQ=255)
                                                                         caACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCa                                                                          <  1:5417713/151‑1 (MQ=255)
                                                                         caACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTAACCCCCAACCACAACCGCACGAGGGAGATGGGGCCCGGCGTGCCTGGTGGGTAGTGCGCCACGCTCGGgcgcc                                                                          >  2:10479029/1‑150 (MQ=255)
                                                                          aaCCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAg                                                                         <  2:5022138/151‑1 (MQ=255)
                                                                                ccccGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCACCCCAAACCACGACCCCCCCCGGGGGAGGGGGCCGGGCGGGGCTTGGGGGTCGCGCGCCACGCTCGGGCGCAGCCAggg                                                                   >  1:9195457/1‑150 (MQ=255)
                                                                                ccccGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCACCCCAAACCACAACCGCACGCGGGAGAGGGGGCCGGGCGGGCCCGGTGGGTCGAGCGCCGCGCACGGGCGCCGCCAGGa                                                                   >  2:10318000/1‑151 (MQ=255)
                                                                                ccccGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCAGACCCCAAACCACACCCCCACCCGGGGGAGGGGGGCGGGGGGGGCGTGTGGGGGGCGCGCCCAGCGCGGGGGCAGCCAGGa                                                                   >  1:8281475/1‑151 (MQ=255)
                                                                                 cccGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCACAACCCACCCCACTCACGCCCGCGGGAGAGGGGGCCGGGGGTGGCGTGTGGGGGGCGCGCCGCGCACGGGCGCCGCCggggc                                                                  >  1:4122157/1‑149 (MQ=255)
                                                                                  ccGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCACCCCAACCCACAACCCCACCCGGGGGTGGGGGGCGGGGGTGCCGTGTGGGTGGGGCGCCCCGCGCGGGGGCAGCCgggtaa                                                                 >  1:7510251/1‑148 (MQ=255)
                                                                                    gctgctAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCACCCCCAACCACACCCCCACCCCGGGGGGGGGGGCGGGGGGGGCGTGGGGGGGGGGGGCCGCGCGCGGGGGCCGCCGGGGACGa                                                               >  1:4778729/1‑151 (MQ=255)
                                                                                     ctgctAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTACACCCAACCCACTCCCCCACGCGGGAGAGGGGGGCGGGCGTGCCGTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCAGGAACgtt                                                              >  2:4432162/1‑149 (MQ=255)
                                                                                       gctAGCCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGt                                                            <  1:6749182/151‑1 (MQ=255)
                                                                                         tAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCACCCCAAACCACTCCCCCACCCGGGCGAGGGGGGCGGGCGTGCCGTGTTGGTCGAGCGCCACGCGAGGGCGCCGCCAGGCACGAttggt                                                          >  1:4648659/1‑147 (MQ=255)
                                                                                         tAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCt                                                          <  2:10581337/151‑1 (MQ=255)
                                                                                         tAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGCAACCCAACCCACAACCCCACCCGGGGGGGGGGGGCGGGGGTGGCGGGTGGGTCGAGCGCCACGGGCGGGCGCCGCCGGGGACGATGtgt                                                          >  1:6453535/1‑149 (MQ=255)
                                                                                          aaCCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGACACCCCACCCACCACCCCCCCCCGGGGTGGGGGGCGGGGGGGGCCTGGGGGGGGGGGGGCCCCGCACGGGGGCGCCCGGGAAcaggtgtt                                                          >  1:4916349/1‑144 (MQ=255)
                                                                                            ccGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTAcc                                                       <  1:1932786/151‑1 (MQ=255)
                                                                                             cGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTAccc                                                      <  2:8496646/151‑1 (MQ=255)
                                                                                              gCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCACCCCAACCCAAAACCCCACCCGGGGGAGGGGGGCCGGGGGGGCGTGTGGGTGGAGCGCCGCCCGAGGGCGCAGCCAGGGACGATGtgaaacca                                                     >  1:9502729/1‑144 (MQ=255)
                                                                                               cAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTc                                                    <  1:2325494/151‑1 (MQ=255)
                                                                                               cAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTc                                                    <  1:8056455/151‑1 (MQ=255)
                                                                                               cAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTc                                                    <  2:3087230/151‑1 (MQ=255)
                                                                                                aaTAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTcc                                                   <  2:2923366/151‑1 (MQ=255)
                                                                                                 aTAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCACACCCCAAACCACTCCCCCACCCGGGAGCAGGGGCCGGGGGTTCCTTGTGGGGACAGCGCCACGCACGGGCGCAGCCGGGAACGATGTCTAACCAccc                                                  >  2:11217261/1‑151 (MQ=255)
                                                                                                   aGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCt                                                <  2:1529747/151‑1 (MQ=255)
                                                                                                          gCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGcccc                                         <  2:5500969/151‑1 (MQ=255)
                                                                                                                gtgtCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTAt                                   <  2:2063140/151‑1 (MQ=255)
                                                                                                                 tgtCGCGGATCGGTTTTTGGCATGTGGTGTTGGGCCGCACCCCAACCCAATCCCCCACGCGGGCGTGGGGGGCGGGGGTGCCGTGTGGGTCGAGCGCCAGGCTCGGGCGCAGCCAGGAACGAGGTCTACCCACCCCAAGGGCCCGGGTaag                                  >  2:3350478/1‑149 (MQ=255)
                                                                                                                 tgtCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATc                                  <  2:3989766/151‑1 (MQ=255)
                                                                                                                 tgtCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATc                                  <  2:5626588/151‑1 (MQ=255)
                                                                                                                     gcgGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGa                              <  1:4561515/151‑1 (MQ=255)
                                                                                                                      cgGATCGGTTTTTGGAATGTGGTGTTGGCCCTAACCCCCACCCGCACCCGCACCCGGGAGAGGGGGCCGGGCGTGGCTTGGGGGGCGAGCGCCGCGCTCGGGCGCCGCCGGGTACGATGTCTACCCGCCCCAAGGCCCCGGGGATCAGGaa                             >  2:4513182/1‑151 (MQ=255)
                                                                                                                      cgGATCGGTTTTTGGAATGTGGTGTTGGCCCGCACCCCAAACCACGCCCCCACGCGGGAGAGGGGGCCCGGCGTGCCGTGTGGGTGGAGCGCCGCGCTCGGGCGCAGCCAGGAACGATGTATAACCGCCCCAAGGCCCCGGGGATCGGGac                             >  2:11026829/1‑150 (MQ=255)
                                                                                                                        gATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAAtt                           <  2:4704500/151‑1 (MQ=255)
                                                                                                                            ggTTTTTGGAATGTGGTGTTGGCCCACACCCCAACCCACACCCCCACCCGGGAGAGGGGGGCCGGCGTGGCTTGTGGGTCGAGCGCCACGCTCGGGCGCAGCCTGGCAAGATGTCTACCCACCCCTAGGCCCCGGGGTTCTGGaaataggg                       >  2:2617713/1‑145 (MQ=255)
                                                                                                                             gTTTTTGGAATGTGGTGTTGGCCCTCACCCCAAACCACTACCCCACCCCGGAGAGGGGGCCGGGCCTTGCGTTTTGGTCGAGGGCCACGCACGGGGGCACACAGGAAAGATGTAAAACCACCCCGAAGGCCCGGGGGTCTTGAAATAGGcc                      >  2:1780769/1‑151 (MQ=255)
                                                                                                                               ttttGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGcccc                    <  2:8228357/151‑1 (MQ=255)
                                                                                                                                 ttGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCccca                  <  1:8863016/151‑1 (MQ=255)
                                                                                                                                  tGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCcccag                 <  1:2020832/151‑1 (MQ=255)
                                                                                                                                   ggAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCcccagc                <  1:5917183/151‑1 (MQ=255)
                                                                                                                                        gtggtgTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAgg           <  2:7226117/151‑1 (MQ=255)
                                                                                                                                         tggtgTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGc          <  1:1389933/151‑1 (MQ=255)
                                                                                                                                         tggtgTTGGCCCACACCCCAACCCACAACCCCACCCGGGAGAGGGGGCCGGGCGTGCCATGTGGGTCGAGCGCCAACCACCGGCGCAGCCAGGAACCAAGTCTACCCAACCCAAGGCCCCCGGGATCTTGAAATACGCCCCCCCCCCACGc          >  1:7357671/1‑151 (MQ=255)
                                                                                                                                         gggtgTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGc          <  1:1059256/150‑1 (MQ=255)
                                                                                                                                           gtgtTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCaa        >  1:5780517/1‑151 (MQ=255)
                                                                                                                                            tgtTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAAc       <  2:8854062/151‑1 (MQ=255)
                                                                                                                                            tgtTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCGAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAAc       <  2:11391202/151‑1 (MQ=255)
                                                                                                                                               tGGCCCTCCCCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCaa    <  1:858506/151‑1 (MQ=255)
                                                                                                                                                 ggCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCaaa   >  1:3290317/1‑151 (MQ=255)
                                                                                                                                                 ggCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCaaa   >  2:5180867/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTATACAGCCCCTCCCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCTTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAA  >  NC_008463/3515715‑3516010

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