Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,870,998 (T)6→5 100% coding (614/1476 nt) yeaI → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,870,9930T.55.0% 11.4 / 254.9 120coding (609/1476 nt)yeaIputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (23/30);  new base . (28/38);  total (52/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

Only 100 of 135 total aligned reads displayed.

ATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGA  >  NC_000913/1870895‑1871090
                                                                                                  |                                                                                                 
aTTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtt                                                                                                  <  1:1925235/100‑1 (MQ=255)
  ttGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttt                                                                                                >  2:2532557/1‑100 (MQ=255)
  ttGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttt                                                                                                <  1:441084/100‑1 (MQ=255)
  ttGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttt                                                                                                <  1:1853897/100‑1 (MQ=255)
   tGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCttttt                                                                                               >  1:2530025/1‑100 (MQ=255)
    gTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCttttta                                                                                              <  1:2558563/100‑2 (MQ=255)
    gTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCttttta                                                                                              >  1:3109753/1‑99 (MQ=255)
    gTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCttttta                                                                                              >  2:3970529/1‑99 (MQ=255)
     tttAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttttag                                                                                             <  2:277181/100‑3 (MQ=255)
     tttAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttttag                                                                                             <  1:3528221/100‑3 (MQ=255)
     tttAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttttag                                                                                             <  1:1645414/100‑3 (MQ=255)
      ttAACCTCTCTGGCGCTCTTTTGTTATGGAAAAGATAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAg                                                                                            <  2:1478702/100‑1 (MQ=255)
      ttAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttttagt                                                                                            <  1:2580626/100‑4 (MQ=255)
        aaCCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACCTCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGttt                                                                                         >  2:1335700/1‑100 (MQ=255)
        aaCCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGttt                                                                                         >  1:1694062/1‑100 (MQ=255)
         aCCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGtttt                                                                                        <  1:4129275/100‑1 (MQ=255)
          ccTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCtttttttttttt                                                                                        >  2:2341616/1‑100 (MQ=255)
          ccTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGtttt                                                                                        >  2:2341573/1‑100 (MQ=255)
          ccTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGttttt                                                                                       <  2:2925013/100‑1 (MQ=255)
          ccTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGttttt                                                                                       >  1:3911677/1‑100 (MQ=255)
           cTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTc                                                                                      >  2:131088/1‑100 (MQ=255)
            tCGCTGGCTCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTcc                                                                                     <  2:623871/100‑1 (MQ=255)
            tCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTc                                                                                      <  2:622202/100‑1 (MQ=255)
            tCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTcc                                                                                     <  1:3517392/100‑1 (MQ=255)
                tGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCt                                                                                  <  2:1666781/100‑1 (MQ=255)
                 ggCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCtt                                                                                 <  1:408312/100‑1 (MQ=255)
                 ggCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTc                                                                                >  1:1015885/1‑100 (MQ=255)
                   cgcTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTc                                                                                >  1:3084201/1‑99 (MQ=255)
                   cgcTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTg                                                                              <  1:2467208/100‑1 (MQ=255)
                   cgcTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTg                                                                              >  2:2670984/1‑100 (MQ=255)
                    gcTCTTTTGTTATGGAAAAGACAACATCCTTGACAAAAAAAAAAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTgg                                                                             <  2:3335729/100‑1 (MQ=255)
                     ctctTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGc                                                                            >  2:3753301/1‑100 (MQ=255)
                     ctctTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGc                                                                            >  2:1908445/1‑100 (MQ=255)
                     ctctTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGc                                                                            <  1:1179411/100‑1 (MQ=255)
                      tctTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGc                                                                            <  1:2826891/100‑1 (MQ=255)
                      tctTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGca                                                                           >  2:741316/1‑100 (MQ=255)
                          ttGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGcacac                                                                        <  2:323933/100‑1 (MQ=255)
                          ttGTTATGGAAAAGACAACATCCTTGAAAAAAAAAAAAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGca                                                                           <  2:1954656/96‑1 (MQ=255)
                               aTGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTg                                                                  <  1:192606/100‑1 (MQ=255)
                               aTGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATACTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTg                                                                  >  2:1614771/1‑100 (MQ=255)
                                tGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTg                                                                  >  2:870377/1‑100 (MQ=255)
                                 ggAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGa                                                                 >  1:4037091/1‑100 (MQ=255)
                                 ggAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAg                                                                >  1:3473827/1‑100 (MQ=255)
                                  gAAAAGACAACATCCTTGACAACAATAAGAACAAAACGGGACTCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTCGCACACAATCagag                                                                <  2:189237/100‑1 (MQ=255)
                                  gAAAAGACAACATCCTTGACAACAATAAGAAAAAAACTGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGc                                                               >  2:2201089/1‑100 (MQ=255)
                                  gAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGc                                                               <  1:3495005/100‑1 (MQ=255)
                                   acaaGACAACATCCTTGAAAACAATAAGAAAAAAACGGGCATCCTGTTGCTGGCGCTGATCCCTTTTTTAGCTTTTCCCCTTCTGGCACACAATCTgcg                                                                <  2:189270/97‑3 (MQ=255)
                                   aaaaGAAAACATCCTTGACAACAATAAGAAAAAAACAGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGc                                                               <  2:282968/100‑1 (MQ=255)
                                      aGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGtt                                                           <  2:2976411/100‑1 (MQ=255)
                                       gACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGtt                                                           >  2:2345096/1‑100 (MQ=255)
                                       gACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGtt                                                           <  2:3164276/100‑1 (MQ=255)
                                        acaacaTCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTta                                                          <  1:1360315/100‑1 (MQ=255)
                                          aacaTCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATaa                                                       <  1:131089/100‑1 (MQ=255)
                                          aacaTCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATaa                                                       >  2:1117152/1‑100 (MQ=255)
                                           acaTCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGGGCAGTTATAAc                                                      <  2:223009/100‑1 (MQ=255)
                                             aTCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGc                                                    <  1:3180962/100‑1 (MQ=255)
                                             aTCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGc                                                    <  2:34751/100‑1 (MQ=255)
                                              tCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCTCCTTTTGGCACACAATCTTAGCAGTtatta                                                       >  2:202316/1‑95 (MQ=255)
                                               ccTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTg                                                  <  1:391919/100‑1 (MQ=255)
                                                cTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTg                                                  <  1:151814/100‑1 (MQ=255)
                                                     caacaaTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTAt                                             >  1:3167405/1‑100 (MQ=255)
                                                        caaTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGACCAGTTATAACGCTCACTATTCt                                          >  2:2562840/1‑100 (MQ=255)
                                                        caaTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCtt                                         >  1:180186/1‑100 (MQ=255)
                                                        caaTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCtt                                         <  2:3110279/100‑1 (MQ=255)
                                                        caaTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCtt                                         >  1:128867/1‑100 (MQ=255)
                                                         aaTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCCTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCtt                                         >  1:510820/1‑100 (MQ=255)
                                                          aTAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCttt                                        <  1:823018/100‑1 (MQ=255)
                                                           tAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGt                                      <  1:1080649/100‑1 (MQ=255)
                                                           tAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGt                                      <  1:2136345/100‑1 (MQ=255)
                                                           tAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGt                                      <  1:2566514/100‑1 (MQ=255)
                                                           tAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGt                                      <  1:4015853/100‑1 (MQ=255)
                                                            aaGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGt                                      >  2:1544763/1‑100 (MQ=255)
                                                            aaGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGt                                      <  1:870383/100‑1 (MQ=255)
                                                             aGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTa                                     >  1:2635280/1‑100 (MQ=255)
                                                             aGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTAt                                    >  2:2103839/1‑100 (MQ=255)
                                                              gAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTAt                                    >  1:312086/1‑100 (MQ=255)
                                                              gAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATg                                   <  2:1092841/100‑1 (MQ=255)
                                                              gAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATg                                   >  1:3876463/1‑100 (MQ=255)
                                                               aaaaaaaCGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATg                                   <  1:1719877/100‑1 (MQ=255)
                                                               aaaaaaaCGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTTAGCAGTTATAACGCTGACTATTCTTTGTATgt                                  <  2:1059264/100‑1 (MQ=255)
                                                                 aaaaaCGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATgtcg                                >  1:1815245/1‑100 (MQ=255)
                                                                   aaaCGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATgtcgt                               <  1:126374/100‑1 (MQ=255)
                                                                   aaaCGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATgtcgtc                              <  1:3244175/100‑1 (MQ=255)
                                                                    aaCGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATgtcgtc                              <  1:3363451/100‑1 (MQ=255)
                                                                      cGGGAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGtcgtcg                             <  2:2017124/98‑1 (MQ=255)
                                                                       gggAATCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGAtt                          >  1:2083147/1‑100 (MQ=255)
                                                                        ggAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGAtt                          <  1:2345109/100‑1 (MQ=255)
                                                                          aaTCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTAc                        <  2:3318347/100‑1 (MQ=255)
                                                                            tCCTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGt                     <  2:2979221/100‑1 (MQ=255)
                                                                              cTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTc                    <  1:3994758/100‑1 (MQ=255)
                                                                              cTGTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTcc                   <  2:180185/100‑1 (MQ=255)
                                                                               tGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTcc                   <  2:3167386/100‑1 (MQ=255)
                                                                                gTTGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCgg                 <  2:3766326/100‑1 (MQ=255)
                                                                                  tGCTGGCGCTGATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGAc               <  2:3731665/100‑1 (MQ=255)
                                                                                   gCTGGCGCTGATCCC‑TTTTTAGTTTTTTCTCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACa              >  2:3405385/1‑100 (MQ=255)
                                                                                         gcTGATCCCTTTTTTAGTTTTTCCCCTTCTGACACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCa          <  2:1344353/99‑1 (MQ=255)
                                                                                            gATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACt      >  1:2865218/1‑100 (MQ=255)
                                                                                            gATCCCTCTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACt      >  1:2166154/1‑100 (MQ=255)
                                                                                            gATCCC‑TTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTg     >  2:3482403/1‑100 (MQ=255)
                                                                                                ccTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGa  >  1:108742/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
ATTTGTTTAACCTCGCTGGCGCTCTTTTGTTATGGAAAAGACAACATCCTTGACAACAATAAGAAAAAAACGGGAATCCTGTTGCTGGCGCTGATCCCTTTTTTAGTTTTTCCCCTTCTGGCACACAATCTGAGCAGTTATAACGCTGACTATTCTTTGTATGTCGTCGATTACTGTCCGGACAACCATACTGCGA  >  NC_000913/1870895‑1871090

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: ‑