Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 332,891 G→A 6.2% R174H (CGC→CAC)  pdeL → DNA‑binding transcriptional activator/c‑di‑GMP phosphodiesterase PdeL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913332,8910GA 6.2% 199.5 / 2.2 65R174H (CGC→CAC) pdeLDNA‑binding transcriptional activator/c‑di‑GMP phosphodiesterase PdeL
Reads supporting (aligned to +/- strand):  ref base G (29/31);  new base A (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

TGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTA  >  NC_000913/332743‑333040
                                                                                                                                                    |                                                                                                                                                      
tGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCCTCCGGTCTTATTGTCATAATGACCCGc                                                                                                                                                       >  1:834358/1‑150 (MQ=255)
  gATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCa                                                                                                                                                     >  2:659494/1‑150 (MQ=255)
   aTCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCaa                                                                                                                                                    >  1:89101/1‑150 (MQ=255)
    tCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAAc                                                                                                                                                   <  2:420723/150‑1 (MQ=255)
    tCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAAc                                                                                                                                                   >  1:1087055/1‑150 (MQ=255)
     ccAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACt                                                                                                                                                  >  1:1502422/1‑150 (MQ=255)
      cAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACtg                                                                                                                                                 >  1:723276/1‑150 (MQ=255)
       aaCCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACtga                                                                                                                                                <  1:1376000/150‑1 (MQ=255)
         ccGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACtgatg                                                                                                                                              <  2:709333/150‑1 (MQ=255)
         ccGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACtgatg                                                                                                                                              <  1:824175/150‑1 (MQ=255)
             ttttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAAc                                                                                                                                          <  1:1415111/150‑1 (MQ=255)
              tttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACa                                                                                                                                         <  2:1588876/150‑1 (MQ=255)
               ttCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAg                                                                                                                                        >  2:1019015/1‑150 (MQ=255)
                      cgcAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCgg                                                                                                                                 <  1:127775/150‑1 (MQ=255)
                              ggCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTg                                                                                                                         >  1:278581/1‑150 (MQ=255)
                                       aCGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTa                                                                                                                <  1:526063/150‑1 (MQ=255)
                                          ggCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTaaaa                                                                                                             >  2:1463578/1‑150 (MQ=255)
                                          ggCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAAATGATGAAACAGACTGCGGATATTCTGATGCCGGTaaaa                                                                                                             >  2:1615509/1‑150 (MQ=255)
                                             tgtgAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACAt                                                                                                          >  1:67365/1‑150 (MQ=255)
                                                gAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTtg                                                                                                       >  1:842932/1‑150 (MQ=255)
                                                 cggTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACGGGATCAGGTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTtgc                                                                                                      <  1:264985/149‑1 (MQ=255)
                                                  ggTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTtgct                                                                                                     >  1:898328/1‑150 (MQ=255)
                                                     gCTTGTCCGCTGGGAACATCCACAAACGGTAATTATCCCACCGGATCAGTTTATTCCTATGCCGGAGTCATCCGGTCTTATTGTCATAATGACCCTCGAACTGATGAAAAAGACTTCGTATATTCTGATGCCGGTAAAACATTTGCTGcc                                                                                                  <  2:194994/150‑1 (MQ=255)
                                                     gCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGcc                                                                                                  <  1:1019015/150‑1 (MQ=255)
                                                      cTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCg                                                                                                 >  2:31056/1‑150 (MQ=255)
                                                            cGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCACCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAAt                                                                                           >  1:460996/1‑150 (MQ=255)
                                                            cGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCACCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAAt                                                                                           >  1:466377/1‑150 (MQ=255)
                                                               tGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTc                                                                                        <  2:723276/150‑1 (MQ=255)
                                                                gggAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTcc                                                                                       <  2:842932/150‑1 (MQ=255)
                                                                 ggAACATCCACAAACGGGAATTATCCCACCGGACCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCa                                                                                      <  2:1536769/150‑1 (MQ=255)
                                                                         cacaAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCa                                                                              <  1:1463578/150‑1 (MQ=255)
                                                                         cacaAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCa                                                                              >  2:962971/1‑150 (MQ=255)
                                                                             aaCGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGAATCaa                                                                          >  2:135690/1‑150 (MQ=255)
                                                                               cGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACg                                                                        >  1:1399359/1‑150 (MQ=255)
                                                                               cGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACg                                                                        <  1:879880/150‑1 (MQ=255)
                                                                               cGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACg                                                                        >  2:1312601/1‑150 (MQ=255)
                                                                                  gAATTATCCCACCGGATCAGTTTATTCCTCTGGCGTAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGtct                                                                     <  1:606145/150‑1 (MQ=255)
                                                                                         cccACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCggg                                                             <  1:1387365/150‑1 (MQ=255)
                                                                                         cccACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCggg                                                             <  2:278581/150‑1 (MQ=255)
                                                                                            aCCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGgttg                                                          <  2:898328/150‑1 (MQ=255)
                                                                                                gATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGtttt                                                      <  2:778103/150‑1 (MQ=255)
                                                                                                         aTTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCg                                             >  2:293794/1‑150 (MQ=255)
                                                                                                         aTTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCg                                             >  1:117286/1‑150 (MQ=255)
                                                                                                         aTTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCg                                             >  2:969666/1‑150 (MQ=255)
                                                                                                           tCCTCTGGCGTCGGCATCCGGTCTTATTGTCAAAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCggg                                           <  2:1502422/150‑1 (MQ=255)
                                                                                                            ccTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGa                                          >  2:973165/1‑150 (MQ=255)
                                                                                                            ccTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGa                                          >  1:1207525/1‑150 (MQ=255)
                                                                                                                 ggCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGa                                     <  1:678119/150‑1 (MQ=255)
                                                                                                                   cGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGaaa                                   <  2:89101/150‑1 (MQ=255)
                                                                                                                        gcATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAagag                              <  1:300013/149‑1 (MQ=255)
                                                                                                                         cATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCggg                                           >  2:365068/1‑136 (MQ=255)
                                                                                                                         cATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCggg                                           <  1:365068/136‑1 (MQ=255)
                                                                                                                           tCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTgg                                                  <  2:253837/127‑1 (MQ=255)
                                                                                                                           tCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTgg                                                  >  1:253837/1‑127 (MQ=255)
                                                                                                                              ggTCTTATTGTCATAATGACCCGCCAACTGTTGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGCCAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAACAGAGTGTCTg                        >  1:1031950/1‑150 (MQ=255)
                                                                                                                                   tATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCt                   >  1:1158480/1‑150 (MQ=255)
                                                                                                                                    aTTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTg                  <  1:1301908/150‑1 (MQ=255)
                                                                                                                                        tCATTATGACCCGCCAACTTATGAAACAGTCTTCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCTGCGGGTTGTTTTTTGGCAGAGGGATTTGAAAAAGAGTGTCTGAACCTGGTta              <  2:1087055/150‑1 (MQ=255)
                                                                                                                                           tAATGAACCGCGAACTGATGAATCAGACTGCGGATATTCTGATGCCGGTAAAAGATTTGCTGCCGCACAATTTCCATATTCGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGACTGTCTGAACCTGGTtaata           <  2:149398/150‑1 (MQ=255)
                                                                                                                                              tGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTGTCTTCTGGTTTTTTTTTTGGCAGCGGGATTTGAAAAAGATTTTCTGAACCTTGGTATTaaa         <  2:834358/150‑1 (MQ=255)
                                                                                                                                                aCCCACCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTa      <  2:466377/150‑1 (MQ=255)
                                                                                                                                                aCCCACCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTa      <  2:460996/150‑1 (MQ=255)
                                                                                                                                                 cccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAg     <  1:31056/150‑1 (MQ=255)
                                                                                                                                                 cccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAg     >  2:112160/1‑150 (MQ=255)
                                                                                                                                                 cccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAg     >  2:1605725/1‑150 (MQ=255)
                                                                                                                                                    gCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTa  <  1:1273131/150‑1 (MQ=255)
                                                                                                                                                    gCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTa  >  1:1268237/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TGGATCCAACCGGTTTTCTGCGCGCAGACTGGCGTACTGACGGGCTGTGAGGTGCTTGTCCGCTGGGAACATCCACAAACGGGAATTATCCCACCGGATCAGTTTATTCCTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTA  >  NC_000913/332743‑333040

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: