Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,644,364 Δ1 bp intergenic (‑115/‑120) sohB ← / → PA14_RS16560 protease SohB/histidine phosphatase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,644,3610G.99.8% 2502.6 / ‑3.8 506intergenic (‑112/‑123)sohB/PA14_RS16560protease SohB/histidine phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (224/281);  minor base C (0/1);  total (224/282)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 527 total aligned reads displayed.

CTGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGATTCG  >  NC_008463/3644212‑3644506
                                                                                                                                                     |                                                                                                                                                 
cTGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAg                                                                                                                                                   <  2:3380098/150‑1 (MQ=255)
 tGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCTgg                                                                                                                                                  >  2:6976695/1‑150 (MQ=255)
    ccAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggca                                                                                                                                               <  2:7759577/150‑3 (MQ=255)
    ccAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggca                                                                                                                                               <  2:9092293/150‑3 (MQ=255)
    ccAGAAAGCCCGAATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGAGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAgggca                                                                                                                                               <  1:6181761/150‑3 (MQ=255)
       gAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGc                                                                                                                                           >  2:10435363/1‑150 (MQ=255)
         aaGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGccc                                                                                                                                         >  1:1022244/1‑150 (MQ=255)
          aGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCg                                                                                                                                        <  2:3954070/150‑1 (MQ=255)
           gCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCgg                                                                                                                                       >  1:9147691/1‑150 (MQ=255)
           gCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCgg                                                                                                                                       <  2:4835120/150‑1 (MQ=255)
            cccGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGa                                                                                                                                      >  2:10560799/1‑150 (MQ=255)
             ccGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCAGACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGag                                                                                                                                     <  2:2833903/150‑1 (MQ=255)
              cGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGagc                                                                                                                                    <  2:7510870/150‑1 (MQ=255)
               gCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGagca                                                                                                                                   <  1:1292900/150‑1 (MQ=255)
               gCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGagca                                                                                                                                   >  1:5035818/1‑150 (MQ=255)
                 atatTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCCGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGagcagc                                                                                                                                 <  2:4229377/150‑1 (MQ=255)
                 atatTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGagcagc                                                                                                                                 <  1:8058305/150‑1 (MQ=255)
                 atatTCGGCAACGAACTCCACACCCCCCCCTCTCCACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGagcagc                                                                                                                                 <  1:2961741/150‑1 (MQ=255)
                   atTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCg                                                                                                                               <  2:7681521/150‑1 (MQ=255)
                     tCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGAt                                                                                                                             <  1:9860504/150‑1 (MQ=255)
                       ggCAATGAACTCCACACCACCCCCTCTCGACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATaa                                                                                                                           <  2:2928384/150‑1 (MQ=255)
                       ggCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATaa                                                                                                                           >  1:9352626/1‑150 (MQ=255)
                       ggCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATaa                                                                                                                           <  2:2330146/150‑1 (MQ=255)
                       ggCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATaa                                                                                                                           <  1:5796936/150‑1 (MQ=255)
                       ggCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATaa                                                                                                                           >  2:4390208/1‑150 (MQ=255)
                        gCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATaaa                                                                                                                          >  2:5638818/1‑150 (MQ=255)
                        gCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATaaa                                                                                                                          <  2:6231636/150‑1 (MQ=255)
                        gCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATaaa                                                                                                                          <  1:2259191/150‑1 (MQ=255)
                         cAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAg                                                                                                                         >  1:1562890/1‑150 (MQ=255)
                         cAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAg                                                                                                                         >  1:3437695/1‑150 (MQ=255)
                          aaTGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAgg                                                                                                                        <  1:8873841/150‑1 (MQ=255)
                           aTGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGc                                                                                                                       <  1:4329142/150‑1 (MQ=255)
                            tGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCa                                                                                                                      >  2:1604822/1‑150 (MQ=255)
                              aaCTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGc                                                                                                                    >  1:4757536/1‑150 (MQ=255)
                               aCTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCg                                                                                                                   <  2:4051225/150‑1 (MQ=255)
                               aCTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCg                                                                                                                   >  1:10704354/1‑150 (MQ=255)
                                  ccACACACCCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAAt                                                                                                                <  1:2400759/150‑1 (MQ=255)
                                  ccACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAAt                                                                                                                <  1:5828610/150‑1 (MQ=255)
                                    acacaACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGa                                                                                                              <  1:7459870/150‑1 (MQ=255)
                                      acaACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACa                                                                                                            <  1:2090324/150‑1 (MQ=255)
                                       caACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACAt                                                                                                           >  1:6699170/1‑150 (MQ=255)
                                        aaCCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACAtt                                                                                                          >  2:3560469/1‑150 (MQ=255)
                                        aaCCCCCTCTCTACTGGTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCCCCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACAtt                                                                                                          <  2:5792125/150‑1 (MQ=255)
                                         aCCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTc                                                                                                         >  1:5240380/1‑150 (MQ=255)
                                         aCCCCCGCTCGACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTc                                                                                                         <  2:4303335/150‑1 (MQ=255)
                                                ctctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAg                                                                                                  <  2:4411612/150‑1 (MQ=255)
                                                ctctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAg                                                                                                  <  2:6367203/150‑1 (MQ=255)
                                                ctctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAg                                                                                                  <  2:3462285/150‑1 (MQ=255)
                                                ctctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAg                                                                                                  <  1:3174812/150‑1 (MQ=255)
                                                 tctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGc                                                                                                 >  2:5290716/1‑150 (MQ=255)
                                                 tctACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGc                                                                                                 >  1:4088080/1‑150 (MQ=255)
                                                        ttGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtt                                                                                          >  1:9654655/1‑150 (MQ=255)
                                                         tGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGttt                                                                                         <  1:2816334/150‑1 (MQ=255)
                                                          gACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtttt                                                                                        <  2:216675/150‑1 (MQ=255)
                                                          gACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtttt                                                                                        <  2:3185480/150‑1 (MQ=255)
                                                           aCGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGttttt                                                                                       <  2:1499475/150‑1 (MQ=255)
                                                            cGCACCCGTCTACGGGCGACACGATTGCGCACGCATACCCAGGCCACCGTCGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtttttt                                                                                      <  2:10105340/150‑1 (MQ=255)
                                                            cGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtttttt                                                                                      >  1:10377309/1‑150 (MQ=255)
                                                            cGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtttttt                                                                                      <  1:1666761/150‑1 (MQ=255)
                                                            cGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtttttt                                                                                      >  1:605586/1‑150 (MQ=255)
                                                            cGCACCCGTCTACGGGCGACACGATTGCGCAAGCAGACCCAGGCCACCGTCGGGCGGCGAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGtttttt                                                                                      <  2:8534285/150‑1 (MQ=255)
                                                              cACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCg                                                                                    <  1:1960390/150‑1 (MQ=255)
                                                                cccGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACAGTCAAACAAGCGTATGTGTTTTCGtt                                                                                  >  1:7750920/1‑150 (MQ=255)
                                                                  cGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACCAGCGTATGTTTTTTCGTTGa                                                                                <  2:3123702/150‑1 (MQ=255)
                                                                     cTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAg                                                                             <  2:6719944/150‑1 (MQ=255)
                                                                         gggCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCt                                                                         >  1:8240818/1‑150 (MQ=255)
                                                                           gCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTgg                                                                       <  2:9427030/150‑1 (MQ=255)
                                                                            cGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTggg                                                                      <  2:9196434/150‑1 (MQ=255)
                                                                             gACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGt                                                                     >  1:10035834/1‑150 (MQ=255)
                                                                                   aTTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCt                                                               <  2:4737823/150‑1 (MQ=255)
                                                                                   aTTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACAGTCAAACAAGCGTAGGTTTTTTCGTTGACAGCCCTGGGTGATCCt                                                               >  2:7144998/1‑150 (MQ=255)
                                                                                     tGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCa                                                             >  2:1271674/1‑150 (MQ=255)
                                                                                       cgcAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATa                                                           <  2:8186429/150‑1 (MQ=255)
                                                                                        gcAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCCGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCAtcc                                                          <  1:3819356/150‑3 (MQ=255)
                                                                                          aaGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATAccc                                                        >  2:8873327/1‑150 (MQ=255)
                                                                                           aGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCt                                                       <  1:5642443/150‑1 (MQ=255)
                                                                                               tACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTcgcg                                                   <  1:8525613/150‑1 (MQ=255)
                                                                                                 cccAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGTCAGCCCTGGGTGATCCTCATACCCTCGCGGa                                                 >  2:10399681/1‑150 (MQ=255)
                                                                                                       ccACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTCTGTTTTTTCGTTGACAGCCCTGGGTGATCCGCATACCCTCGCGGAACTTTc                                           >  2:2466142/1‑150 (MQ=255)
                                                                                                       ccACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTc                                           <  1:7842413/150‑1 (MQ=255)
                                                                                                          ccGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCAcc                                        >  1:7197261/1‑150 (MQ=255)
                                                                                                           cGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCAccc                                       <  2:1442165/150‑1 (MQ=255)
                                                                                                             aCGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACccgc                                     <  1:3611169/150‑1 (MQ=255)
                                                                                                                ggcggcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGa                                  <  1:2616703/150‑1 (MQ=255)
                                                                                                                  cggcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAgg                                <  1:6560955/150‑1 (MQ=255)
                                                                                                                   ggcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGt                               >  2:3663292/1‑150 (MQ=255)
                                                                                                                   ggcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGt                               >  1:9515767/1‑150 (MQ=255)
                                                                                                                    gcCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTc                              <  2:1880732/150‑1 (MQ=255)
                                                                                                                      cAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCg                            >  2:9320172/1‑150 (MQ=255)
                                                                                                                         cGATGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACAGTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCGGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGccc                         >  2:6580881/1‑150 (MQ=255)
                                                                                                                           aTGCCACTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGt                       >  2:7322875/1‑150 (MQ=255)
                                                                                                                                aCTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGgca                  <  1:903578/150‑1 (MQ=255)
                                                                                                                                aCTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGATTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGgca                  <  2:5184656/150‑1 (MQ=255)
                                                                                                                                 cTGCCTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGgcag                 >  1:5396533/1‑150 (MQ=255)
                                                                                                                                    ccTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCAt              >  1:10507663/1‑150 (MQ=255)
                                                                                                                                     cTGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATc             >  1:4610924/1‑150 (MQ=255)
                                                                                                                                      tGGAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCAGCt            >  2:9846251/1‑150 (MQ=255)
                                                                                                                                       ggAGAGCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTa           <  1:4525378/150‑1 (MQ=255)
                                                                                                                                            gCCTTTGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGa      >  2:4981246/1‑150 (MQ=255)
                                                                                                                                                ttGCAGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCAGACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGATTCg  <  2:10283373/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CTGGCCAGAAAGCCCGCATATTCGGCAATGAACTCCACACAACCCCCTCTCTACTGTTGACGCACCCGTCTACGGGCGACACGATTGCGCAAGCATACCCAGGCCACCGACGGGCGGCCAGCGATGCCACTGCCTGGAGAGCCTTTGCAGGGGCACGCCCGGAGCAGCCGATAAAGGCAGCGAATGACATTCAAACAAGCGTATGTTTTTTCGTTGACAGCCCTGGGTGATCCTCATACCCTCGCGGAACTTTCACCCGCCGAGGTCCGCCCGTGGGCAGCATCTACCTGATTCG  >  NC_008463/3644212‑3644506

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