Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,714,073 C→T P141L (CCG→CTG)  yiaC → GNAT family putative N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,714,0730CT95.6% 288.6 / 2.5 91P141L (CCG→CTG) yiaCGNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (51/36);  total (52/39)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-01

CAACCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAG  >  NC_000913/3713976‑3714169
                                                                                                 |                                                                                                
cAACCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTTCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGctggctgg                                                                                                >  2:1417137/1‑100 (MQ=255)
 aaCCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTggt                                                                                               <  2:2039875/100‑1 (MQ=255)
 aaCCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTggt                                                                                               <  2:1471372/100‑1 (MQ=255)
  aCCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTggtg                                                                                              <  2:461414/100‑1 (MQ=255)
  aCCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTggtg                                                                                              <  1:1478233/100‑1 (MQ=255)
   ccGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCtggtgg                                                                                             <  2:2101924/100‑1 (MQ=255)
   ccGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCtggtgg                                                                                             <  1:2318747/100‑1 (MQ=255)
   ccGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCtggtgg                                                                                             >  2:336194/1‑100 (MQ=255)
   ccGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCggtgg                                                                                             >  1:110254/1‑100 (MQ=255)
    cGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTggtggt                                                                                            >  2:429643/1‑100 (MQ=255)
    cGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTggtggt                                                                                            <  2:247224/100‑1 (MQ=255)
    cGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTggtggt                                                                                            <  1:792184/100‑1 (MQ=255)
     ggCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGtt                                                                                           >  2:2249358/1‑100 (MQ=255)
      gCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTc                                                                                          <  2:1222219/100‑1 (MQ=255)
       cGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCa                                                                                         >  2:1207462/1‑100 (MQ=255)
       cGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCa                                                                                         <  1:1924142/100‑1 (MQ=255)
       cGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCa                                                                                         <  2:448730/100‑1 (MQ=255)
            aaTTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGc                                                                                    >  1:2116388/1‑100 (MQ=255)
             aTTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCt                                                                                   >  1:1921108/1‑100 (MQ=255)
             aTTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCt                                                                                   >  1:299241/1‑100 (MQ=255)
              tttttACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTg                                                                                  <  1:2326045/100‑1 (MQ=255)
               ttttACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGt                                                                                 <  1:1474842/100‑1 (MQ=255)
                tttACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTa                                                                                <  2:2817164/100‑1 (MQ=255)
                  tACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAg                                                                              <  1:24341/100‑1 (MQ=255)
                  tACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAg                                                                              <  2:110254/100‑1 (MQ=255)
                    ccAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAgtg                                                                            >  2:2799809/1‑100 (MQ=255)
                       ggCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCa                                                                         >  1:573103/1‑100 (MQ=255)
                       ggCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACACACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCa                                                                         >  2:877194/1‑100 (MQ=255)
                         cacaGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGt                                                                       >  1:682862/1‑100 (MQ=255)
                           caGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTc                                                                     >  2:2289203/1‑100 (MQ=255)
                             gggTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCg                                                                   <  1:1818498/100‑1 (MQ=255)
                                ttttCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGAc                                                                <  1:2799830/100‑1 (MQ=255)
                                ttttCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGAc                                                                >  2:1705136/1‑100 (MQ=255)
                                ttttCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGAc                                                                <  2:289569/100‑1 (MQ=255)
                                  ttCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCg                                                              >  2:2085395/1‑100 (MQ=255)
                                    cacaTTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACAGGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGtt                                                            >  1:1668731/1‑100 (MQ=255)
                                     acaTTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTg                                                           >  2:2175459/1‑100 (MQ=255)
                                     acaTTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTg                                                           <  1:349514/100‑1 (MQ=255)
                                     acaTTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTg                                                           >  1:471761/1‑100 (MQ=255)
                                     acaTCGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTg                                                           <  2:2548494/100‑1 (MQ=255)
                                      caTTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGt                                                          >  2:1818524/1‑100 (MQ=255)
                                       aTTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTa                                                         >  2:1097252/1‑100 (MQ=255)
                                       aTTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTa                                                         >  2:1059903/1‑100 (MQ=255)
                                       aTTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTa                                                         >  2:2711366/1‑100 (MQ=255)
                                          gTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTAttt                                                      >  2:2402576/1‑100 (MQ=255)
                                            cGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTAttttt                                                    >  2:1532194/1‑100 (MQ=255)
                                             gATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTc                                                   <  2:682860/100‑1 (MQ=255)
                                             gATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTc                                                   >  2:878834/1‑100 (MQ=255)
                                             gATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTc                                                   >  1:1184712/1‑100 (MQ=255)
                                               ttGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCa                                                 <  1:839492/100‑1 (MQ=255)
                                                 gcgcATGGCAGGATGAAACCCAACTACCCGCATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGc                                               <  1:2388589/100‑1 (MQ=255)
                                                 gcgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGc                                               >  1:1223059/1‑100 (MQ=255)
                                                 gcgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGc                                               >  1:2418927/1‑100 (MQ=255)
                                                  cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGcc                                              >  1:1394177/1‑100 (MQ=255)
                                                  cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGcc                                              >  1:877573/1‑100 (MQ=255)
                                                  cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGcc                                              >  2:2477073/1‑100 (MQ=255)
                                                     aTGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATg                                           <  2:1534635/100‑1 (MQ=255)
                                                      tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGc                                          <  1:99083/100‑1 (MQ=255)
                                                      tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGc                                          <  2:2079308/100‑1 (MQ=255)
                                                        gCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCa                                        <  1:1554962/100‑1 (MQ=255)
                                                          aGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAgc                                      <  1:2394565/100‑1 (MQ=255)
                                                          aGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAgc                                      >  1:1984640/1‑100 (MQ=255)
                                                            gATGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAgcgc                                    <  1:878835/100‑1 (MQ=255)
                                                             aTGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGcc                                   >  2:2777323/1‑100 (MQ=255)
                                                             aTGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGcc                                   >  2:678423/1‑100 (MQ=255)
                                                             aTGAAACCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGcc                                   >  2:678842/1‑100 (MQ=255)
                                                                aaaCCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTa                                >  1:298705/1‑100 (MQ=255)
                                                                 aaCCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTAt                               >  1:2178246/1‑100 (MQ=255)
                                                                 aaCCCAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTAt                               <  2:104598/100‑1 (MQ=255)
                                                                   cccAACTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTa                             >  2:2066894/1‑100 (MQ=255)
                                                                       aCTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCg                         <  2:299240/100‑1 (MQ=255)
                                                                        cTACCCACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCg                         >  2:1398133/1‑99 (MQ=255)
                                                                           cccACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGcaca                     <  2:573101/100‑1 (MQ=255)
                                                                           cccACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGcaca                     >  2:1420269/1‑100 (MQ=255)
                                                                           cccACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGcac                      >  1:1676664/1‑99 (MQ=255)
                                                                            ccACATGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAg                    <  2:378294/100‑1 (MQ=255)
                                                                             cacaTGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGc                   >  2:2241909/1‑100 (MQ=255)
                                                                             cacaTGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGc                   <  1:2679297/100‑1 (MQ=255)
                                                                             cacaTGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGc                   >  1:2096860/1‑100 (MQ=255)
                                                                             cacaTGGATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGc                   >  2:1439063/1‑100 (MQ=255)
                                                                                 tGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCAt                <  2:1983760/99‑1 (MQ=255)
                                                                                  ggATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTc              >  2:1280372/1‑100 (MQ=255)
                                                                                  ggATTATGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTc              <  1:1911553/100‑1 (MQ=255)
                                                                                       aTGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAg         >  2:2529242/1‑100 (MQ=255)
                                                                                       aTGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAg         >  1:565859/1‑100 (MQ=255)
                                                                                       aTGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAg         <  2:1102179/100‑1 (MQ=255)
                                                                                       aTGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAg         <  1:1916763/100‑1 (MQ=255)
                                                                                        tGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTTTTTTGTTTTTAGCTCCGTATTACCCGCACAGCTAtt               >  2:1573917/1‑93 (MQ=255)
                                                                                        tGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGc        >  2:1510704/1‑100 (MQ=255)
                                                                                        tGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGc        >  1:2612372/1‑100 (MQ=255)
                                                                                        tGAGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGc        >  1:2731919/1‑100 (MQ=255)
                                                                                          aGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCcg      >  1:1699667/1‑100 (MQ=255)
                                                                                          aGCTGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCcg      <  1:2249371/100‑1 (MQ=255)
                                                                                           gctggctgGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCcgc     >  1:1175669/1‑100 (MQ=255)
                                                                                           gctggctgGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCcgc     >  2:1450365/1‑100 (MQ=255)
                                                                                            ctggctggTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCcgcg    <  1:2510064/100‑1 (MQ=255)
                                                                                            ctggctggTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCcgcg    <  1:2777344/100‑1 (MQ=255)
                                                                                             tGGCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGctcgcg   <  2:2221385/100‑6 (MQ=255)
                                                                                              ggCTGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAg  >  2:202489/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CAACCGGCGATAAATTTTTACCAGGCACAGGGTTTTCACATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCCGCGAG  >  NC_000913/3713976‑3714169

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