Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,804 G→T 17.5% intergenic (+431/+63) icd → / ← ymfD isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase/e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,8040GT17.5% 276.9 / 46.5 109intergenic (+431/+63)icd/ymfDisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase/e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (42/48);  new base T (10/9);  total (52/57)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

Only 100 of 112 total aligned reads displayed.

AACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAA  >  NC_000913/1196708‑1196903
                                                                                                |                                                                                                   
aaCTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAgcgc                                                                                                  >  2:1032807/1‑100 (MQ=255)
 actactCACCCGGGAGTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACAAAAAAAAAAAAAAAAAGGCCGCGAgcg                                                                                                   <  2:2087173/98‑1 (MQ=255)
   tactCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGc                                                                                               >  2:890811/1‑100 (MQ=255)
     ctCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGcct                                                                                             >  1:200567/1‑99 (MQ=255)
     cgcACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGcct                                                                                             >  1:2373552/3‑99 (MQ=255)
       cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCt                                                                                           >  1:2237344/1‑100 (MQ=255)
       cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCt                                                                                           <  2:2528058/100‑1 (MQ=255)
        aCCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGcctttt                                                                                          <  1:548461/100‑5 (MQ=255)
         cccGGTACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGccttttt                                                                                         <  2:359966/100‑6 (MQ=255)
         cccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAACAAAGGCCGCGAGCGCGGCCCCTTc                                                                                         >  2:1008837/1‑100 (MQ=255)
         cccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGccttttt                                                                                         <  1:2258678/100‑6 (MQ=255)
         cccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCtt                                                                                          >  1:515791/1‑99 (MQ=255)
         cccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTc                                                                                         <  2:1210929/100‑1 (MQ=255)
         cccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTc                                                                                         <  1:692368/100‑1 (MQ=255)
          ccGGGACTCTCCAGGGTACAGCCAACAGGCATTGGTTGGAATCACCTTAGCGTTCAGGTACATGCGCAATGTAAAAAAGGCCGCGCTCGCGGcctttttt                                                                                        <  2:1721814/100‑7 (MQ=255)
          ccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGcctttttt                                                                                        <  1:2135657/100‑7 (MQ=255)
          ccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTca                                                                                        <  2:1664889/100‑1 (MQ=255)
           cGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGcctttttta                                                                                       <  2:2373537/100‑8 (MQ=255)
           cGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTcac                                                                                       <  1:2340403/100‑1 (MQ=255)
            gggACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGccttttttac                                                                                      <  1:1267409/100‑9 (MQ=255)
            gggACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGccttttttac                                                                                      <  1:1150647/100‑9 (MQ=255)
            gggACTCGCCAGGGGACAGCCAACAGGCATTGGGAGCAATCACCTTCGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGccgtttttac                                                                                      <  2:1341685/100‑9 (MQ=255)
              gACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCacata                                                                                    <  2:177979/100‑1 (MQ=255)
              gACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATATAAAAAAGACCGCGAGCGCGGCCCCTTCacata                                                                                    >  2:935559/1‑100 (MQ=255)
                cTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCacataca                                                                                  >  1:1661603/1‑100 (MQ=255)
                cTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCacataca                                                                                  >  2:290379/1‑100 (MQ=255)
                 tCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCTCTTctca                                                                                      <  2:1320757/95‑1 (MQ=255)
                       ggggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAg                                                                           >  2:750015/1‑100 (MQ=255)
                         ggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTa                                                                         <  2:358260/100‑1 (MQ=255)
                          gACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTAc                                                                        <  1:771516/100‑1 (MQ=255)
                          gACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTAc                                                                        <  2:2071199/100‑1 (MQ=255)
                            cAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTg                                                                      >  1:1653159/1‑100 (MQ=255)
                             aGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTga                                                                     <  2:1197628/100‑1 (MQ=255)
                             aGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTga                                                                     <  2:362611/100‑1 (MQ=255)
                              gCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTgag                                                                    >  1:2380860/1‑100 (MQ=255)
                              gCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTgag                                                                    <  2:780963/100‑1 (MQ=255)
                              gCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTgag                                                                    >  1:1342531/1‑100 (MQ=255)
                              gCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTgag                                                                    >  1:1163254/1‑100 (MQ=255)
                               ccAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTgaga                                                                   >  2:2000919/1‑100 (MQ=255)
                               ccAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTgaga                                                                   <  2:699069/100‑1 (MQ=255)
                                cAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGAc                                                                  >  1:281735/1‑100 (MQ=255)
                                  aCAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTg                                                                >  1:992875/1‑100 (MQ=255)
                                  aCAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTg                                                                <  1:537614/100‑1 (MQ=255)
                                   cAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGt                                                               <  2:1461735/100‑1 (MQ=255)
                                     ggCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGttt                                                             <  1:2009072/100‑1 (MQ=255)
                                      gCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGttt                                                             <  2:317544/99‑1 (MQ=255)
                                        aTTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAAc                                                          >  2:1293570/1‑100 (MQ=255)
                                         ttGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAAcc                                                         <  1:1240927/100‑1 (MQ=255)
                                         ttGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAAcc                                                         <  2:902821/100‑1 (MQ=255)
                                           gggTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTACCCTa                                                       >  2:1183055/1‑100 (MQ=255)
                                             gTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAgg                                                     >  1:303102/1‑100 (MQ=255)
                                             gTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAgg                                                     >  1:2051098/1‑100 (MQ=255)
                                                 aaTCACCTTAGCGTTCAGGTACATGCGGAATGTAAAACAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAat                                                 >  2:1366785/1‑100 (MQ=255)
                                                 aaTCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAat                                                 >  1:404474/1‑100 (MQ=255)
                                                  aTCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAatt                                                <  2:499560/100‑1 (MQ=255)
                                                    cACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAattat                                              <  1:173746/100‑1 (MQ=255)
                                                      ccTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAAttatta                                            >  1:2542214/1‑100 (MQ=255)
                                                        ttAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATc                                          <  1:1552317/100‑1 (MQ=255)
                                                        ttAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATc                                          >  1:1213575/1‑100 (MQ=255)
                                                            cGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGtat                                      <  2:1163246/100‑1 (MQ=255)
                                                            cGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGtat                                      >  1:1829341/1‑100 (MQ=255)
                                                            cGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGtat                                      <  1:2285367/100‑1 (MQ=255)
                                                              ttCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGtatta                                    <  2:1960650/100‑1 (MQ=255)
                                                                  ggTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATAtt                                >  1:686722/1‑100 (MQ=255)
                                                                    tACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGc                              <  2:1143762/100‑1 (MQ=255)
                                                                      cATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCAt                            >  1:2491128/1‑100 (MQ=255)
                                                                       aTGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGGATTATTATCGTATTATATTGCAtt                           >  2:377677/1‑100 (MQ=255)
                                                                          cGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACa                        <  1:643620/100‑1 (MQ=255)
                                                                             aaTGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAAc                      <  2:2212170/99‑1 (MQ=255)
                                                                             aaTGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACt                     >  1:1890512/1‑100 (MQ=255)
                                                                              aTGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACTTCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATCACAACTc                    >  2:1374595/1‑100 (MQ=255)
                                                                              aTGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTc                    <  2:323378/100‑1 (MQ=255)
                                                                              aTGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTc                    <  1:637881/100‑1 (MQ=255)
                                                                              aTGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTc                    <  1:303356/100‑1 (MQ=255)
                                                                               atgtgaaggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTcc                   >  1:1115405/10‑100 (MQ=255)
                                                                                  tgaaggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCaa                >  2:1573419/7‑100 (MQ=255)
                                                                                  aaaaaaGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCaa                >  1:2119135/1‑100 (MQ=255)
                                                                                   aaaaaGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAAc               >  1:205038/1‑100 (MQ=255)
                                                                                   aaaaaGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAAc               >  1:513209/1‑100 (MQ=255)
                                                                                   aaaaaGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAAc               <  1:635830/100‑1 (MQ=255)
                                                                                   aaaaaGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAAc               <  2:1407114/100‑1 (MQ=255)
                                                                                     aggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGCCTGTTTAACCTAGGGAATTATTATAGTATTATATTGCATTACAACTCCAAACAg             >  2:1910465/4‑100 (MQ=255)
                                                                                     aggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAg             >  1:1410733/4‑100 (MQ=255)
                                                                                      aaGGCCGCGAGCGCGGCACCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCCACAg             >  1:1264357/1‑99 (MQ=255)
                                                                                       gggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTg           >  2:1510523/2‑100 (MQ=255)
                                                                                       gggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTg           >  1:959766/2‑100 (MQ=255)
                                                                                       aGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTg           <  2:1318797/100‑1 (MQ=255)
                                                                                        ggCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGa          <  2:176720/100‑1 (MQ=255)
                                                                                         gCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGAc         <  2:2195050/100‑1 (MQ=255)
                                                                                         gCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGAc         >  1:2045566/1‑100 (MQ=255)
                                                                                         gCCACGAGCGCGGCCCCTTCACATACATCTTTAGTACTGTGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTg           >  2:1755340/1‑98 (MQ=255)
                                                                                          ccGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCCACAGTGACa        >  1:2443289/1‑100 (MQ=255)
                                                                                          ccGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACa        <  2:2476640/100‑1 (MQ=255)
                                                                                          ccGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACa        <  1:1278308/100‑1 (MQ=255)
                                                                                          ccGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAAATCCCAACCGTGACa        >  1:212024/1‑100 (MQ=255)
                                                                                           cgcgAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACat       >  2:30341/1‑100 (MQ=255)
                                                                                            gcgAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACata      <  1:1877094/100‑1 (MQ=255)
                                                                                             cgAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACatat     <  1:833680/100‑1 (MQ=255)
                                                                                              gAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTTTTTAACCTAGTGAATTATTATCGTATTATATTTCATTACATCTTCCAACAGTGACATATg    <  2:559590/100‑1 (MQ=255)
                                                                                                gcgcGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGaa  >  1:2099468/1‑100 (MQ=255)
                                                                                                |                                                                                                   
AACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAA  >  NC_000913/1196708‑1196903

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