Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,516,423 A→T 5.1% H470Q (CAT→CAA pdeA ← putative c‑di‑GMP phosphodiesterase PdeA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,516,4230AT 5.1% 256.3 / 2.8 79H470Q (CAT→CAApdeAputative c‑di‑GMP phosphodiesterase PdeA
Reads supporting (aligned to +/- strand):  ref base A (48/27);  new base T (2/2);  total (50/29)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

CATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCA  >  NC_000913/2516275‑2516571
                                                                                                                                                    |                                                                                                                                                    
catAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACat                                                                                                                                                     <  1:361658/150‑1 (MQ=255)
catAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACat                                                                                                                                                     >  2:1437795/1‑150 (MQ=255)
catAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACat                                                                                                                                                     >  2:361658/1‑150 (MQ=255)
 atAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGGCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACatg                                                                                                                                                    <  2:453362/150‑1 (MQ=255)
   aaCCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACatgat                                                                                                                                                  <  2:93929/150‑1 (MQ=255)
    aCCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACatgatg                                                                                                                                                 <  1:979606/150‑1 (MQ=255)
     ccTTCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACatgatga                                                                                                                                                <  1:672389/150‑1 (MQ=255)
       ttCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACatgatgatg                                                                                                                                              <  1:1402072/150‑1 (MQ=255)
        tCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGAtgt                                                                                                                                             >  1:504261/1‑150 (MQ=255)
        tCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGAtgt                                                                                                                                             >  2:1391580/1‑150 (MQ=255)
          aCCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTgc                                                                                                                                           >  2:1045897/1‑150 (MQ=255)
                       cGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTcc                                                                                                                              >  2:345374/1‑150 (MQ=255)
                       cGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTcc                                                                                                                              >  2:11264/1‑150 (MQ=255)
                        gAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCg                                                                                                                             <  1:1426560/150‑1 (MQ=255)
                         aaTTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACTTGATGATGTGCGCAGGATTGCTCCgc                                                                                                                            >  1:1504945/1‑150 (MQ=255)
                         aaTTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACTTGATGATGTGCGCAGGATTGCTCCgc                                                                                                                            >  1:1510060/1‑150 (MQ=255)
                             ggCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCgccagc                                                                                                                        >  1:34285/1‑150 (MQ=255)
                             ggCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCgccagc                                                                                                                        <  2:1055971/150‑1 (MQ=255)
                             ggCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTAACCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTACgccagc                                                                                                                        >  2:464356/1‑150 (MQ=255)
                                 gggCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCcagccagc                                                                                                                    <  2:504261/150‑1 (MQ=255)
                                    cGTAGAGCAGCAAATCACACTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTa                                                                                                                 >  1:1288997/1‑150 (MQ=255)
                                          gcagcaAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGcc                                                                                                           >  1:178783/1‑150 (MQ=255)
                                          gcagcaAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGcc                                                                                                           >  2:514669/1‑150 (MQ=255)
                                             gcaAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCca                                                                                                        <  2:105943/150‑1 (MQ=255)
                                              caAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCcaa                                                                                                       >  2:614057/1‑150 (MQ=255)
                                                aaTCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCcaaca                                                                                                     >  1:1225604/1‑150 (MQ=255)
                                                aaTCACCCTGATCTAACGCCCTCCGCATGGTATTACGCAATAGCACCTGTTTAGTGGTCTCCCCGCTACCCATCTCCTCTCAGCTATCCAGCGCCAGCACATGATGATGTGCCCAGGATTGCTCCGACAGCCCGATTAACTGCCCcaaca                                                                                                     >  2:1164655/1‑150 (MQ=255)
                                                     cccTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGt                                                                                                >  2:570105/1‑150 (MQ=255)
                                                     cccTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGt                                                                                                <  1:1045897/150‑1 (MQ=255)
                                                       cTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTg                                                                                              >  2:1532110/1‑150 (MQ=255)
                                                        tGATCTAACGCCGTGCGAATGGTATTCAGCAATATCACCTGTTTAGTGTTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGc                                                                                             >  1:58422/1‑150 (MQ=255)
                                                           tCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAgg                                                                                          <  2:178783/150‑1 (MQ=255)
                                                             tAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGt                                                                                        >  2:1327033/1‑150 (MQ=255)
                                                                   cGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTg                                                                                  <  1:402132/150‑1 (MQ=255)
                                                                    gTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGa                                                                                 >  1:615688/1‑150 (MQ=255)
                                                                       cGAATGGTATTCAGCAATATCACCTGTTTAGTGGTCTGCCAGGAAACCATCTCCTCTCTGCTATCAAGAGCCATAAAATGAGGATGTGCGCAGGATTTCTCCGCCATCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTTAttt                                                                              <  2:1048277/150‑1 (MQ=255)
                                                                       cGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGAttt                                                                              <  1:1201960/150‑1 (MQ=255)
                                                                        gAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACTTGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTc                                                                             <  2:1504945/150‑1 (MQ=255)
                                                                        gAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACTTGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTc                                                                             <  2:1510060/150‑1 (MQ=255)
                                                                            ggTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATc                                                                         <  2:1218020/150‑1 (MQ=255)
                                                                               aTTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCaaa                                                                      >  1:1491597/1‑150 (MQ=255)
                                                                                ttCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCaaaa                                                                     >  2:1071806/1‑150 (MQ=255)
                                                                                 tCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAAc                                                                    >  1:681766/1‑150 (MQ=255)
                                                                                  cAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACg                                                                   >  2:1325191/1‑150 (MQ=255)
                                                                                  cAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACg                                                                   >  2:363512/1‑150 (MQ=255)
                                                                                      aaTAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGcccc                                                               <  1:345374/150‑1 (MQ=255)
                                                                                       aTAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGccccc                                                              >  1:1226043/1‑150 (MQ=255)
                                                                                                   tAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCAt                                                  >  1:1513685/1‑150 (MQ=255)
                                                                                                    aGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATa                                                 >  2:372919/1‑150 (MQ=255)
                                                                                                      tGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGc                                               <  1:1325191/150‑1 (MQ=255)
                                                                                                       ggTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGcc                                              >  2:1474748/1‑150 (MQ=255)
                                                                                                       ggTCTGCCCGGAAACCATCTCCTCTCTGCTATCAAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCCGCCACCTTAACTGCCCCAACAACGGATGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGcc                                              >  1:945099/1‑150 (MQ=255)
                                                                                                        gTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGccc                                             >  1:1243380/1‑150 (MQ=255)
                                                                                                         tCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCa                                            <  1:570105/150‑1 (MQ=255)
                                                                                                          cTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCat                                           <  2:1225604/150‑1 (MQ=255)
                                                                                                            gCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCatat                                         >  2:147910/1‑150 (MQ=255)
                                                                                                             cccGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATc                                        <  1:11264/150‑1 (MQ=255)
                                                                                                              ccGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATcc                                       >  1:9794/1‑150 (MQ=255)
                                                                                                                ggAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAg                                     <  2:935170/150‑1 (MQ=255)
                                                                                                                 gAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGc                                    >  1:503146/1‑150 (MQ=255)
                                                                                                                    aCCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCg                                 <  2:584258/150‑1 (MQ=255)
                                                                                                                    aCCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCg                                 <  2:578778/150‑1 (MQ=255)
                                                                                                                     ccATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCgg                                <  1:197936/150‑1 (MQ=255)
                                                                                                                      cATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGt                               >  1:754258/1‑150 (MQ=255)
                                                                                                                         ctcctcTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTAtt                            >  2:1531878/1‑150 (MQ=255)
                                                                                                                          tcctcTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTg                           <  1:1474748/150‑1 (MQ=255)
                                                                                                                          tcctcTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTg                           >  2:614654/1‑150 (MQ=255)
                                                                                                                          tcctcTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTg                           <  2:571/150‑1 (MQ=255)
                                                                                                                           cctcTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGt                          >  1:866977/1‑150 (MQ=255)
                                                                                                                           cctcTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGt                          >  2:44467/1‑150 (MQ=255)
                                                                                                                           cctcTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGt                          >  2:1035010/1‑150 (MQ=255)
                                                                                                                            ctctctGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGtt                         >  2:268122/1‑150 (MQ=255)
                                                                                                                             tctctGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTc                        <  1:200961/150‑1 (MQ=255)
                                                                                                                                ctGCTATCCAGCGCCAGCACATGATGATGTGCGCAGTATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAg                     >  1:1330384/1‑150 (MQ=255)
                                                                                                                                ctGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAg                     <  1:933311/150‑1 (MQ=255)
                                                                                                                                   cTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGa                  >  2:1488346/1‑150 (MQ=255)
                                                                                                                                     aTCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAAt                <  1:113827/150‑1 (MQ=255)
                                                                                                                                      tCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAAtt               >  2:863964/1‑150 (MQ=255)
                                                                                                                                       ccAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAAttt              >  2:648535/1‑150 (MQ=255)
                                                                                                                                        cAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTg             >  2:627767/1‑150 (MQ=255)
                                                                                                                                         aGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGc            >  2:1289696/1‑150 (MQ=255)
                                                                                                                                          gcgcCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTccc           <  1:965456/150‑1 (MQ=255)
                                                                                                                                           cgcCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCg          <  2:503146/150‑1 (MQ=255)
                                                                                                                                            ctgcgCACATCATCATGTGTGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGa         >  1:601537/5‑150 (MQ=255)
                                                                                                                                             tgcgCACATCATCATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGAc        >  2:356402/4‑150 (MQ=255)
                                                                                                                                             ccAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGAc        >  2:1141719/1‑150 (MQ=255)
                                                                                                                                               aGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTa      <  2:34285/150‑1 (MQ=255)
                                                                                                                                                  acaTGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTc   <  1:1532110/150‑1 (MQ=255)
                                                                                                                                                   caTGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCa  >  2:439384/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATAACCTTCACCCTCTTTGTTGCGAATTGGCTGGGCGTAGAGCAGCAAATCACCCTGATCTAACGCCGTGCGAATGGTATTCAGCAATAGCACCTGTTTAGTGGTCTGCCCGGAAACCATCTCCTCTCTGCTATCCAGCGCCAGCACATGATGATGTGCGCAGGATTGCTCCGCCAGCCAGCTTAACTGCCCCAACAAGGGTTGCAGGGTTTCCTGATTTCCATCAAAACGCCCCCAGGCAGCACCATAGCCCATATCCAGCCCGGTATTGTTCCAGTGAATTTGCCGACTATTCA  >  NC_000913/2516275‑2516571

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