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.1% 1267.6 / NA 328intergenic (+617/+267)PA14_RS16005/PA14_RS16010AAA family ATPase/DUF4354 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (131/194);  minor base A (3/0);  total (134/194)
Fisher's exact test for biased strand distribution p-value = 6.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 328 total aligned reads displayed.

GTATACAGCCCCTCCCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCTTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAA  >  NC_008463/3515715‑3516009
                                                                                                                                                     |                                                                                                                                                  
   taCAGCCCCTCCCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCAcc                                                                                                                                               >  2:10741547/1‑151 (MQ=255)
    aCAGCCCCTCCCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCAccc                                                                                                                                              <  2:5492530/151‑1 (MQ=255)
    aCAGCCCCTCCCTAAGCCATAGCCCTTCCCCAGACCTGAACGGCACACCTAGTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCAccc                                                                                                                                              <  2:6597300/151‑1 (MQ=255)
      aGCCCCTCCCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCACAccccc                                                                                                                                            >  2:1521867/1‑150 (MQ=255)
         ccctccctAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCTCACCCCCAcc                                                                                                                                         >  1:7429643/1‑151 (MQ=255)
         ccctccctAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAAcc                                                                                                                                         <  1:6184459/151‑1 (MQ=255)
          cctccctAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGAACCCCAAccc                                                                                                                                        >  1:1384717/1‑151 (MQ=255)
               ctAAGCCAGAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTcc                                                                                                                                   <  2:1283443/151‑1 (MQ=255)
                 aaGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGCACCCCCACCCTCTcccc                                                                                                                                 >  1:9674156/1‑150 (MQ=255)
                    ccATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCAc                                                                                                                              <  2:2282225/151‑1 (MQ=255)
                      aTAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCAcgc                                                                                                                            <  1:10023120/151‑1 (MQ=255)
                       tAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCCCCCCCACCCACTCCCGCAcgcg                                                                                                                           >  2:8948453/1‑151 (MQ=255)
                         gCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCggg                                                                                                                         <  1:9005242/151‑1 (MQ=255)
                         gCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCggg                                                                                                                         <  2:2620477/151‑1 (MQ=255)
                          cccTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGga                                                                                                                        <  1:2205637/151‑1 (MQ=255)
                              tCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCACACCCCCACCCGGGCGAgg                                                                                                                    >  1:6449381/1‑151 (MQ=255)
                              tCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAAACCTATCCCGCACGCGGGTGAgg                                                                                                                    >  2:5268226/1‑151 (MQ=255)
                               ccccATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCTCCCCCCAAACCACACCCCCACCCGGGAGAggg                                                                                                                   >  1:4612633/1‑151 (MQ=255)
                                cccATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCCCCCCAACCCACTCCCGCACGCGGGAGAgggg                                                                                                                  >  2:5444632/1‑151 (MQ=255)
                                 ccATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAggggg                                                                                                                 <  1:9155754/151‑1 (MQ=255)
                                   aTACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGcc                                                                                                               <  2:7162449/151‑1 (MQ=255)
                                    tACCTGATCGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTATTGGAATGTGGTGTTGGCCCGCCACCCAAACCGAGACCGCCCGCGGGTGTGGGGGCCg                                                                                                              >  1:2608114/1‑151 (MQ=255)
                                     aCCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCgg                                                                                                             <  2:9565722/151‑1 (MQ=255)
                                            cGGCACACCTATTCCCACCCGCAAAACCACACCCTCCCCCGCTGCTAACCGCAATCGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCCCGCGGGAGAGGGGGCCGGGCGTGcc                                                                                                      <  1:5943178/151‑1 (MQ=255)
                                             ggCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCt                                                                                                     <  1:6226378/151‑1 (MQ=255)
                                              gCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCtt                                                                                                    <  1:3192860/151‑1 (MQ=255)
                                               cacaCCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTtg                                                                                                   <  2:9119483/151‑1 (MQ=255)
                                                acacCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGGCACAAACCCAACCCCGCCCGCCCCCGGGGGTGGGGGCCCGGCGTGccgggt                                                                                                  >  1:4163465/1‑147 (MQ=255)
                                                   ccTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCACACCCCAACCCACACCCGCCCGCGGGGGAGGGGGGCGGGCGTGCCGTGTggg                                                                                               >  1:10335726/1‑151 (MQ=255)
                                                     tATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTGTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCGTGTGGGTc                                                                                             >  2:4611112/1‑151 (MQ=255)
                                                      aTTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCg                                                                                            <  2:8742659/151‑1 (MQ=255)
                                                       ttGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGa                                                                                           <  1:4913087/151‑1 (MQ=255)
                                                       ttGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCACACCCCAACCCCCGCACGCCCGCGGGAGAGGGGGCCGGGCGTGCCTggggggggga                                                                                           >  2:1275710/1‑147 (MQ=255)
                                                        tGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAg                                                                                          <  2:4456133/151‑1 (MQ=255)
                                                         gCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCTCACCCCAAACCAATCCCGCACCCGGGAGAGGGGGCCGGGGGTGCCGTGTGGGGAGAgc                                                                                         >  2:738970/1‑151 (MQ=255)
                                                         gCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAgc                                                                                         <  2:9293933/151‑1 (MQ=255)
                                                         gCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAgc                                                                                         <  1:8116357/151‑1 (MQ=255)
                                                          ccACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCTCACCCCAACCCTCTCCCGCCCCCGGGCGGGGGGGCCGGGCGTGCCTTGTGGGTCGAgcg                                                                                        >  2:6950751/1‑151 (MQ=255)
                                                           cACCCGCAACACCACAACCTCCCCCGCTGCCACCCGCAATAGCTGCGGCTAGAGGGACGCGGAGCGGTTGTTGGAAGGAGGTGGAGGCCCGCCCCCCCCCCCACTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCATGGGGGTCGAgcgc                                                                                       <  2:11190414/151‑1 (MQ=255)
                                                             cccGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCct                                                                                     <  1:5512548/151‑1 (MQ=255)
                                                              ccGCAAAACCACCACCTCCCCCGCGGCTCACCGCAAGAGCTGCGGCTTGAGTGTCGCGGATCGGATTTTGGAATGTGGTGTATGCCCTCACCCCACCCCACTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCctc                                                                                    <  2:4245686/151‑1 (MQ=255)
                                                               cGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCctcg                                                                                   <  1:3342625/151‑1 (MQ=255)
                                                                   aaaCCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCctcgctcg                                                                               <  2:10354583/151‑1 (MQ=255)
                                                                    aaCCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCCCCCCCACCCGCAGCCGCCCGCGGGGGTGGGGGGCGGGCGGGGCGTGGGGGGCGGGCGCCTCGCTCgg                                                                              >  2:11409738/1‑151 (MQ=255)
                                                                    aaCCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGCACCCCCACCCTCTCCCCCACGCGGGCGCGGGGGGCGGGCGGGCCGTGTTGGTGGTGCGCCGCGCACgg                                                                              >  2:6422740/1‑151 (MQ=255)
                                                                    aaCCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCACACCCCCACCCACAACCCCACGCCGGCGTGGGGGGCCGGCGTGCCATGTGGGTGGGGCGCcacgcaagg                                                                              >  1:1764247/1‑147 (MQ=255)
                                                                    aaCCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCAAAACCCAAACCACACCCCCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGGCGAGCGCCTCGCTAgg                                                                              >  2:2788710/1‑151 (MQ=255)
                                                                      ccACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGGCCCACCCCCACCCACAACCCCCCCCCGGGGAGGGGGGCGGGGGGGCCGGGTGGGGGGGGCGCCCCGCGCgggg                                                                            >  1:4062095/1‑150 (MQ=255)
                                                                         caACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCCCGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCa                                                                         <  1:5500463/151‑1 (MQ=255)
                                                                             cTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGCACCCCAACCCACTACCCCACGCGGGAGAGGGGGCCGGGCGTGCCATGTGGGTCGAGCGCCTCGCTCGGGCGCCGCCt                                                                     >  2:9400485/1‑151 (MQ=255)
                                                                             cTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCACACCCCAAACCACTACCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGccg                                                                     >  2:2545033/1‑150 (MQ=255)
                                                                              tCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTg                                                                    <  1:10580839/151‑1 (MQ=255)
                                                                              tCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTg                                                                    <  1:2334026/151‑1 (MQ=255)
                                                                                  ccGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCACTACCGCACCCGGGAGAGGGGGCCCGGCGGGCCGTGGGGGTCGAGCGCCACGCTCGGGCGCAGCCTGGTAc                                                                >  2:11738334/1‑151 (MQ=255)
                                                                                     ctgctAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGGCGCACCCCAACCCACAACCCCACCCCGGGGAGGGGGGCGGGCGGGGCGTGGGGGGGGAGCGCCGCGGGCGGGCGCGGCCGGGacacat                                                             >  1:296302/1‑146 (MQ=255)
                                                                                      tgctAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATg                                                            <  2:10807975/151‑1 (MQ=255)
                                                                                      tgctAACCGCAATAGCTGCGGATTGAGTGTCGCGGATCGGTTTTTGGTATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATg                                                            <  2:4892151/151‑1 (MQ=255)
                                                                                        ctAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCACCCCAACCCACACCCCCCCCCGGGGGTGGGGGGCGGGGGTGGCGGGGGGGGCGGGCGCCCCGCGCGGGGGCAGCCAGGGACGATGta                                                          >  1:8525676/1‑150 (MQ=255)
                                                                                        ctAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCGCAACCCAACCCACAACCCCACGCGGGAGAGGGGGGCGGGGGTGGCTTGTTGGTAGAGCGCCACGCACGGGGGCAACCAGGAACGAttta                                                          >  1:8052841/1‑148 (MQ=255)
                                                                                          aaCCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCACACCCCAAACCACAACCGCACGCGGGGGGTGGGGGCGGGCGTGGCTTGTGGGTGGAGCGCCCCGCAAGGGCGCCGCCAGGGACGAtttata                                                        >  2:7789169/1‑146 (MQ=255)
                                                                                           aCCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTAc                                                       <  1:8759128/151‑1 (MQ=255)
                                                                                            ccGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTAcc                                                      <  2:3792269/151‑1 (MQ=255)
                                                                                             cGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCACACCCCAAACCAAAACCCCACCCGGGGGAGGGGGGCGGGGGGTGCGTGTGGGGGGAGCGCCCCCCGCGGGCGCAGCCCGGGAAGATGTGTAccc                                                     >  1:7750931/1‑151 (MQ=255)
                                                                                             cGCAAGAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTAccc                                                     <  2:570973/151‑1 (MQ=255)
                                                                                              gCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCCCACCCCCAACCACACCCCCCCCCGGGGGAGGGGGGCGGGCGGGCCGTGTGGGTGGAGCGCCACGCGCGGGGGCCCCCAGGGAAGATGTGTAacca                                                    >  1:2532852/1‑150 (MQ=255)
                                                                                              gCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTGGGCCCTCACCCCACCCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCt                                                    <  1:7051829/151‑1 (MQ=255)
                                                                                               cAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGGCCACACCCCAACCCACTCCCGCACCCGGGCGTGGGGGGCGGGCGGGCCGTGTGGGTGGCGCGCCGCGCTCGGGGGCAGCCGGGGACGATGTATAaccac                                                   >  2:7987830/1‑149 (MQ=255)
                                                                                                 aTAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCGCAACCCAACCCACACCCCCCCGCCGGCTATGGGGCGCGGCGTGCCATGGGGGTCGAGAGCCACGCACGGGCGCAGCCAGGAACGAAGACTAcccaacc                                                 >  2:6094313/1‑147 (MQ=255)
                                                                                                   aGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCt                                               <  1:5645729/151‑1 (MQ=255)
                                                                                                     cTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAg                                             <  2:3100379/151‑1 (MQ=255)
                                                                                                      tGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAgg                                            <  1:1426173/151‑1 (MQ=255)
                                                                                                        cGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGcc                                          <  1:11212410/151‑1 (MQ=255)
                                                                                                         ggCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGccc                                         <  1:10741547/151‑1 (MQ=255)
                                                                                                          gCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTGGGCCCTCACCCCACCCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGcccc                                        <  2:311989/151‑1 (MQ=255)
                                                                                                           cTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCg                                       <  1:5467760/151‑1 (MQ=255)
                                                                                                                gtgtCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTAt                                  <  2:10604045/151‑1 (MQ=255)
                                                                                                                gtgtCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTAt                                  <  2:4727255/151‑1 (MQ=255)
                                                                                                                 tgtCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATc                                 <  1:11022089/151‑1 (MQ=255)
                                                                                                                  gtCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCTCCCCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCGACCCTCCCCTAGGCCCCGGGTATCt                                <  2:11336183/151‑1 (MQ=255)
                                                                                                                    cgcgGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCCgg                              <  1:2843830/151‑1 (MQ=255)
                                                                                                                     gcgGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGa                             <  1:9854476/151‑1 (MQ=255)
                                                                                                                     gcgGATCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGa                             <  2:2869602/151‑1 (MQ=255)
                                                                                                                         aTCGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTc                         <  1:5369486/151‑1 (MQ=255)
                                                                                                                          tCGGTTTTTGGAATGTGGTGTTGGCCCACACCCCAAACCACTACCCCACGCGGGAGAGGGGGGCCGGCGGGCCGTGGGGGTGGAGCGCCACGCACGGGCGCAACCATGAAAGATTTAAAACCACCCCAACGCCCCCGGTATCAGGACTTCg                        >  2:6796621/1‑151 (MQ=255)
                                                                                                                           cGGTTTTTGGACTGTGGTGTTGGCCCGCCCCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCCGGGGGGGCGTGTGGGGGGCGCGCCTCGCGCGGGCGCAGCCGGGGACGGTGGCGACCCGCCCCGAGGCCCCGGGGATCGGGACTTCgg                       >  2:8162039/1‑151 (MQ=255)
                                                                                                                           cGGTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCgg                       <  2:2820324/151‑1 (MQ=255)
                                                                                                                            ggTTTTTGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGc                      <  1:11409738/151‑1 (MQ=255)
                                                                                                                              tttttGGAATGTGGTGTTGGCCCTAACCCCCACCCACTACCCCCCCCCGGATAGGGGGGCGGGCGTGCCATTTGGGTCGCGCGCCACCCGCGGGCGCACCCAGGAACCAATTCTAACCAACCCATCGCCCCGGGTAAATTGAAATCGGccc                    >  2:175492/1‑151 (MQ=255)
                                                                                                                                tttGGAATGTGGTGTTGGCCCAAACCCCAAACCAAACCCCCACCCGGGTGTGGGGGCCGGGCGTGCCGTGTGGGTGGGGCGCCTCGCGCCGGCGCAGCCAGGGACGATGTCTACCCAACCCAAGGCCCCCGGGATCGTGAATTAGGccccc                  >  2:6254742/1‑151 (MQ=255)
                                                                                                                                 ttGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCccca                 <  1:10476896/151‑1 (MQ=255)
                                                                                                                                 ttGGAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCccca                 <  1:7329876/151‑1 (MQ=255)
                                                                                                                                  tGGAATGTGGTGTTGGCCCTCACCCCCACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCcccag                <  1:6857403/151‑1 (MQ=255)
                                                                                                                                  gggACTGTGGTGTTGGCCCTCACCCCCCCCCGCTCCCGCCCGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGAGGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCcccag                <  2:11539699/150‑1 (MQ=255)
                                                                                                                                   ggAATGTGGTGTTGGCCCTCACCCCACCCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCcccagc               <  1:5226262/151‑1 (MQ=255)
                                                                                                                                   ggAATGTGGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCcccagc               <  1:7833315/151‑1 (MQ=255)
                                                                                                                                       tgtgGTGTTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCcccagcccag           <  1:357505/151‑1 (MQ=255)
                                                                                                                                           gtgtTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCaa       <  1:9535847/151‑1 (MQ=255)
                                                                                                                                             gtTGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAAcc     <  2:1621100/151‑1 (MQ=255)
                                                                                                                                              gtGGCCCGCACCCCCACCCGCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGCGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCGCGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCa    <  2:3309760/150‑1 (MQ=255)
                                                                                                                                               tGGCCCTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCaa   >  2:10120788/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GTATACAGCCCCTCCCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCTTCACCCCAACCCTCTCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCCTGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAA  >  NC_008463/3515715‑3516009

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: