Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,202,693 C→G 10.1% P268R (CCG→CGG)  cca → fused tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase and phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,202,6930CG10.1% 216.2 / 13.1 69P268R (CCG→CGG) ccafused tRNA nucleotidyltransferase/2',3'‑cyclic phosphodiesterase/2' nucleotidase and phosphatase
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base G (4/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CCCTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAACCCAAA  >  NC_000913/3202544‑3202840
                                                                                                                                                     |                                                                                                                                                   
cccTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGcc                                                                                                                                                     >  1:534923/1‑150 (MQ=255)
cccTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGcc                                                                                                                                                     >  2:777805/1‑150 (MQ=255)
   tGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCcgcg                                                                                                                                                  <  1:725416/150‑1 (MQ=255)
    gCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGt                                                                                                                                                 <  2:639368/150‑1 (MQ=255)
        aGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTGTGGCGGCGtcagc                                                                                                                                             <  2:856258/150‑3 (MQ=255)
        aGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCGGCCAGAGTTCTGGCGGCGtcagc                                                                                                                                             <  2:373652/150‑3 (MQ=255)
           ggCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatg                                                                                                                                          >  1:200216/1‑150 (MQ=255)
              aTCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtc                                                                                                                                       >  2:803019/1‑150 (MQ=255)
               tCCGGAAATCGATACGGGTATTAATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCAAGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtaatggtca                                                                                                                                      >  2:999628/1‑150 (MQ=255)
                ccGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtcat                                                                                                                                     >  2:139356/1‑150 (MQ=255)
                ccGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCAtcatggtcat                                                                                                                                     >  2:1222004/1‑150 (MQ=255)
                    aaaTCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGcc                                                                                                                                 <  2:661873/150‑1 (MQ=255)
                     aaTCGAAACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGccc                                                                                                                                >  2:1091614/1‑150 (MQ=255)
                           tACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCggg                                                                                                                          <  2:403652/150‑1 (MQ=255)
                           tACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCggg                                                                                                                          <  1:352341/150‑1 (MQ=255)
                            aCGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgt                                                                                                                         <  1:858559/150‑1 (MQ=255)
                             cGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGt                                                                                                                                                 >  2:787784/1‑125 (MQ=255)
                             cGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGt                                                                                                                                                 <  1:787784/125‑1 (MQ=255)
                             cGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                                                                                                                        <  2:805066/150‑1 (MQ=255)
                             cGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtg                                                                                                                        <  1:597886/150‑1 (MQ=255)
                              gggTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGgtgt                                                                                                                       <  1:1109669/150‑1 (MQ=255)
                               ggTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGtt                                                                                                                      >  2:165121/1‑150 (MQ=255)
                                  aTTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAg                                                                                                                   <  1:311922/150‑1 (MQ=255)
                                     cATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTa                                                                                                                >  2:595484/1‑150 (MQ=255)
                                      aTACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAg                                                                                                               >  2:1186414/1‑150 (MQ=255)
                                       tACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGt                                                                                                              >  1:524074/1‑150 (MQ=255)
                                        aCCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTg                                                                                                             <  2:945542/150‑1 (MQ=255)
                                        aCCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGGCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTg                                                                                                             <  1:316271/150‑1 (MQ=255)
                                          cTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGa                                                                                                           >  1:539589/1‑150 (MQ=255)
                                          cTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGa                                                                                                           >  1:544430/1‑150 (MQ=255)
                                                     tctcAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGcc                                                                                                <  2:206261/150‑1 (MQ=255)
                                                      ctcAATGGCGGCGATGCTGAGTCCACAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGGGTTAAGTTAGTGGAACAACTATGCCa                                                                                               <  2:975272/150‑1 (MQ=255)
                                                         aaTGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGtt                                                                                                                      <  2:275053/124‑1 (MQ=255)
                                                         aaTGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCAAGGTCATGGCCCGGCGGGTGtt                                                                                                                      >  1:275053/1‑124 (MQ=255)
                                                          aTGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGt                                                                                           >  1:668412/1‑150 (MQ=255)
                                                               ggcgATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTgcg                                                                                      <  2:805765/150‑1 (MQ=255)
                                                                    tGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGc                                                                                 >  2:847265/1‑150 (MQ=255)
                                                                    tGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGc                                                                                 >  2:1342997/1‑150 (MQ=255)
                                                                             cGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAAGATCTATGCCAGCGTCTGCGCGTGGCAAATGaaa                                                                        <  2:279316/150‑1 (MQ=255)
                                                                                  gTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTcgc                                                                   <  1:1342997/150‑1 (MQ=255)
                                                                                     gATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCGAGCGTCTGCGCGTGCCAAATGAAATCCGCGAt                                                                <  1:762778/150‑1 (MQ=255)
                                                                                     gATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGAt                                                                <  1:595484/150‑1 (MQ=255)
                                                                                        gTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTa                                                             <  1:1222004/150‑1 (MQ=255)
                                                                                          ccGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGc                                                           >  2:581581/1‑150 (MQ=255)
                                                                                            gTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCa                                                         <  2:534923/150‑1 (MQ=255)
                                                                                                  cAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACtgg                                                   >  2:896573/1‑150 (MQ=255)
                                                                                                           gCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGTCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGt                                          <  1:777805/150‑1 (MQ=255)
                                                                                                            ccACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGtt                                         >  1:498570/1‑150 (MQ=255)
                                                                                                             cACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGttt                                        <  1:136525/150‑1 (MQ=255)
                                                                                                                gATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCAc                                     <  1:139356/150‑1 (MQ=255)
                                                                                                                   ctcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGGGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGAt                                  <  2:763644/150‑1 (MQ=255)
                                                                                                                   ctcGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGAt                                  >  2:893710/1‑150 (MQ=255)
                                                                                                                     cGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGAtct                                <  2:200216/150‑1 (MQ=255)
                                                                                                                     cGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGAtct                                >  1:564778/1‑150 (MQ=255)
                                                                                                                       gTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCa                              >  1:1343406/1‑150 (MQ=255)
                                                                                                                        tAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCAt                             >  1:1082593/1‑150 (MQ=255)
                                                                                                                          aaGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATcc                           <  1:590795/150‑1 (MQ=255)
                                                                                                                           aGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCca                          <  1:165121/150‑1 (MQ=255)
                                                                                                                               atGACGCGGCCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACAcc                      >  1:485050/2‑150 (MQ=255)
                                                                                                                                 gACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGa                                   >  1:384128/1‑135 (MQ=255)
                                                                                                                                 gACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGa                                   <  2:384128/135‑1 (MQ=255)
                                                                                                                                  aCGCGGCCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTc                   >  1:868344/1‑150 (MQ=255)
                                                                                                                                   cgcggcCAGAGTTCTGGCGGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTcc                  >  2:1337021/1‑150 (MQ=255)
                                                                                                                                   cgcggcCAGAGTTCTGGCGGCGTCATCACGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTcc                  >  2:1058380/1‑150 (MQ=255)
                                                                                                                                     ccgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCa                >  1:154814/1‑150 (MQ=255)
                                                                                                                                      cgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCaa               <  1:803019/150‑1 (MQ=255)
                                                                                                                                      cgccAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCaa               >  1:395318/1‑150 (MQ=255)
                                                                                                                                              ctctGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAAc       <  2:524074/150‑1 (MQ=255)
                                                                                                                                              ctctGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAAc       <  1:1186414/150‑1 (MQ=255)
                                                                                                                                                   gCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACAcc                      >  1:1195606/1‑130 (MQ=255)
                                                                                                                                                   gCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACAcc                      <  2:1195606/130‑1 (MQ=255)
                                                                                                                                                   gCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAACCCaaa  >  1:631780/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCCTGCCAAGTGGCATCCGGAAATCGATACGGGTATTCATACCTTAATGACGCTCTCAATGGCGGCGATGCTGAGTCCGCAGGTCGATGTCCGTTTCGCAACTTTATGCCACGATCTCGGTAAAGGGCTGACGCCGCCAGAACTCTGGCCGCGTCATCATGGTCATGGCCCGGCGGGTGTTAAGTTAGTGGAACAACTATGCCAGCGTCTGCGCGTGCCAAATGAAATTCGCGATTTAGCCAGACTGGTGGCTGAGTTTCACGATCTCATCCACACCTTCCCAATGCTGAACCCAAA  >  NC_000913/3202544‑3202840

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: