Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,250,702 T→C 50.5% L329P (CTG→CCG)  yeiH → UPF0324 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,250,7020TC50.5% ‑0.4 / 135.9 95L329P (CTG→CCG) yeiHUPF0324 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (20/27);  new base C (18/30);  total (38/57)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

CCTTCCTGCTGGCAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGTTAACCCGCTATCA  >  NC_000913/2250603‑2250800
                                                                                                   |                                                                                                  
ccTTCCTGCTGGCAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACt                                                                                                    >  2:2492229/1‑100 (MQ=255)
 cTTCCTGCTGGCAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTg                                                                                                   <  2:1026123/100‑1 (MQ=255)
      tgctgGCAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGc                                                                                              >  1:951923/1‑100 (MQ=255)
       gctgGCAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGcc                                                                                             <  2:1303624/100‑1 (MQ=255)
       gctgGCAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGcc                                                                                             <  2:1303420/100‑1 (MQ=255)
       gctgGCAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGcc                                                                                             >  1:1591757/1‑100 (MQ=255)
        ctgGCAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGcc                                                                                             >  2:1128127/1‑99 (MQ=255)
        ctgGCAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCt                                                                                            >  1:580733/1‑100 (MQ=255)
        ctgGCAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCt                                                                                            <  1:394565/100‑1 (MQ=255)
          ggCAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTgg                                                                                          >  1:202945/1‑100 (MQ=255)
          ggCAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTg                                                                                           <  2:1424102/99‑1 (MQ=255)
            cAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCt                                                                                        <  1:1920530/100‑1 (MQ=255)
             aattggcatggtgGCGCTTGGTCTTACCTCTCACTTCAGCGTGCTGAAAACAGCTGGGGCGAAACCTCTGCTGATTGCACTGGTACTGTGTTCCTGGCTg                                                                                       <  2:2234871/93‑1 (MQ=255)
             aaTGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTg                                                                                       <  1:771681/100‑1 (MQ=255)
               tggcgatggcgGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGAt                                                                                     <  1:1083703/100‑1 (MQ=255)
                    atggcgGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGAttgttg                                                                                <  2:545834/100‑1 (MQ=255)
                    atggcgGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGAttgttg                                                                                <  2:196380/100‑1 (MQ=255)
                     tggcgGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTtgg                                                                               >  2:2091264/1‑100 (MQ=255)
                      ggcggcGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTtggt                                                                              <  1:2026153/100‑1 (MQ=255)
                      ggcggcGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTtggt                                                                              >  1:1442416/1‑100 (MQ=255)
                         ggcgCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTtggtggt                                                                           <  2:1881714/100‑1 (MQ=255)
                          gcgcTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTtggtggtg                                                                          <  1:443676/100‑1 (MQ=255)
                             cggtgTCTGATCACTCACGTCAGCACGCTGACAAAAGCTTGGGCGAAACGTCTGCTGATGGCACTGGTACCGGTTGCCTGGCTGATTGTtggtg                                                                             <  2:2314001/90‑1 (MQ=255)
                             cTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGgtggtggtg                                                                       <  1:1184259/100‑1 (MQ=255)
                             cTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGgtggtggtg                                                                       >  1:2211504/1‑100 (MQ=255)
                              tGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGc                                                                      <  1:69881/100‑1 (MQ=255)
                              ccgatCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGc                                                                      <  1:2244479/96‑1 (MQ=255)
                               gggTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCg                                                                     <  2:2237606/100‑1 (MQ=255)
                                ggTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGa                                                                    >  1:434897/1‑100 (MQ=255)
                                ggTCTGACCACGCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCg                                                                     <  2:1094454/99‑1 (MQ=255)
                                 gTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGAt                                                                   >  1:1068729/1‑100 (MQ=255)
                                 gTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGAt                                                                   >  2:1694732/1‑100 (MQ=255)
                                   cTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTa                                                                 <  1:327113/100‑1 (MQ=255)
                                     gACCACTCACTTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAAc                                                               <  2:1241525/100‑1 (MQ=255)
                                     gACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGcagatcgg                                                               >  2:2244466/1‑93 (MQ=255)
                                     gACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAAc                                                               >  1:603596/1‑100 (MQ=255)
                                       ccACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTa                                                             >  2:1831594/1‑100 (MQ=255)
                                       ccACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTa                                                             <  2:2230884/100‑1 (MQ=255)
                                         aCTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTAtg                                                           <  1:1103637/100‑1 (MQ=255)
                                         aCTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTAtg                                                           <  1:496802/100‑1 (MQ=255)
                                            cACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGa                                                        >  1:2216048/1‑100 (MQ=255)
                                             aCGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGGTGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGAt                                                       <  1:944784/100‑1 (MQ=255)
                                             aCGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGAt                                                       <  1:2246960/100‑1 (MQ=255)
                                                tCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGAt                                                       <  2:402135/97‑1 (MQ=255)
                                                   gcgcgcTGAGAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCaaa                                                  <  2:1153790/99‑1 (MQ=255)
                                                   gcgcgcTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAg                                                 <  1:584114/100‑1 (MQ=255)
                                                    cgcgcTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGc                                                <  1:2492242/100‑1 (MQ=255)
                                                    cgcgcTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGc                                                >  1:666086/1‑100 (MQ=255)
                                                     gcgcTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCg                                               <  2:1614867/100‑1 (MQ=255)
                                                     gcgcTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCg                                               <  1:1218409/100‑1 (MQ=255)
                                                       gcTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTc                                             <  2:2441380/100‑1 (MQ=255)
                                                       gcTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTc                                             <  1:1955141/100‑1 (MQ=255)
                                                        cTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTCACTATGTGGTTCAAAACGTc                                             >  2:1617459/1‑99 (MQ=255)
                                                          gAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATc                                          <  2:760169/100‑1 (MQ=255)
                                                          gAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATc                                          >  2:1333998/1‑100 (MQ=255)
                                                          gAAAAAAGCTGGGGCGAAACCCCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCAt                                           >  1:519021/1‑99 (MQ=255)
                                                           aaaaaaGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCg                                         <  1:1308021/100‑1 (MQ=255)
                                                            aaaaaGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGc                                        <  1:190679/100‑1 (MQ=255)
                                                            aaaaaGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGc                                        >  2:1358988/1‑100 (MQ=255)
                                                             aaaaGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCa                                       <  1:2524706/100‑1 (MQ=255)
                                                             aaaaGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCa                                       <  2:2352133/100‑1 (MQ=255)
                                                             aaaaGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCa                                       <  2:2214210/100‑1 (MQ=255)
                                                               aaGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATGAACTATGTGATTCAAAGCGTCATCGCAt                                      <  2:523386/99‑1 (MQ=255)
                                                                  cTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAAc                                  >  1:868630/1‑100 (MQ=255)
                                                                   tGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAAcc                                 <  2:1188657/100‑1 (MQ=255)
                                                                   tGGGGCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAAcc                                 <  1:1568665/100‑1 (MQ=255)
                                                                     gggCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCAc                               <  1:2398748/100‑1 (MQ=255)
                                                                     gggCGAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCAc                               <  1:604716/100‑1 (MQ=255)
                                                                      ggCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACt                              >  1:1500974/1‑100 (MQ=255)
                                                                       gCGAAACCTTTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCAc                               <  2:682354/98‑1 (MQ=255)
                                                                        ctaaaCCTCTGCTGAGGGCACTGTGACCGTTCGCCTGGCTGATTGTGGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTAc                            <  2:1577326/98‑1 (MQ=255)
                                                                         gAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACa                           >  2:113902/1‑100 (MQ=255)
                                                                         gAAACCTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACa                           <  1:434943/100‑1 (MQ=255)
                                                                          aaaCCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACAt                          <  2:1053469/100‑1 (MQ=255)
                                                                             ccTCTGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCtt                       >  2:1214024/1‑100 (MQ=255)
                                                                               tctGCTGATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGc                     <  1:2228751/100‑1 (MQ=255)
                                                                                ctgctgATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTTGTGCTATTAAAAAAAAAATAAAAAGCGTCATCGCATAAACCACTACATCTTg                      <  2:578210/98‑1 (MQ=255)
                                                                                 tgctgATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTc                   >  1:1267897/1‑100 (MQ=255)
                                                                                  gctgATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTcc                  >  2:2389206/1‑100 (MQ=255)
                                                                                  gctgATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTcc                  <  2:1848255/100‑1 (MQ=255)
                                                                                  gctgATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTcc                  >  1:680878/1‑100 (MQ=255)
                                                                                  gctgATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATGAACTATGTGATTCAAAGCGGCATCGCATAAACCACTACATCTTGCTcc                  >  2:1075919/1‑100 (MQ=255)
                                                                                    tgATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCt                 >  1:400025/1‑99 (MQ=255)
                                                                                     gATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGt               <  2:1574719/100‑1 (MQ=255)
                                                                                     gATGGCACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGt               >  1:312468/1‑100 (MQ=255)
                                                                                       tGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGTTa             >  1:420076/1‑100 (MQ=255)
                                                                                          cACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGTTAAcc          >  2:899070/1‑100 (MQ=255)
                                                                                          cACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGTTAAcc          >  1:259920/1‑100 (MQ=255)
                                                                                          cACTGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGTTAAcc          >  2:1336314/1‑100 (MQ=255)
                                                                                           aCTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTTGTGCGATTAACTATGTGATTCCAAGCGTCCTCGCATAAACCACTACATCTTGCTCCTGTTAAcc          <  2:121435/99‑1 (MQ=255)
                                                                                             tGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGTTAACCCGc       >  2:2598684/1‑100 (MQ=255)
                                                                                             tGGTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGTTAACCCGc       <  2:2394222/100‑1 (MQ=255)
                                                                                              ggTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGTTAACCCGCt      >  1:123160/1‑100 (MQ=255)
                                                                                               gTACCGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGTTAACCCGCTa     >  1:555255/1‑100 (MQ=255)
                                                                                                 aCTGTTTGCCTGGCTGATTGTTGGTGGTGTTGCGATTAACTATGTTATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGTTAACCCGCTATc   <  2:1060113/100‑1 (MQ=255)
                                                                                                  cTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGTTAACCCGCTATCa  >  2:2302100/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CCTTCCTGCTGGCAATGGCGATGGCGGCGCTGGGTCTGACCACTCACGTCAGCGCGCTGAAAAAAGCTGGGGCGAAACCTCTGCTGATGGCACTGGTACTGTTTGCCTGGCTGATTGTTGGTGGTGGTGCGATTAACTATGTGATTCAAAGCGTCATCGCATAAACCACTACATCTTGCTCCTGTTAACCCGCTATCA  >  NC_000913/2250603‑2250800

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