Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 594,946 C→A 15.8% Q166H (CAG→CAT) ‡ cusR ← response regulator in two‑component regulatory system with CusS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913594,9460CA15.8% 376.9 / 51.0 138Q166H (CAG→CATcusRresponse regulator in two‑component regulatory system with CusS
Reads supporting (aligned to +/- strand):  ref base C (55/61);  new base A (5/17);  total (60/78)
Fisher's exact test for biased strand distribution p-value = 3.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 146 total aligned reads displayed.

GCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATC  >  NC_000913/594847‑595045
                                                                                                   |                                                                                                   
 cgcAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGccc                                                                                                     <  2:3461007/99‑1 (MQ=255)
     gCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATg                                                                                                <  1:1716620/100‑1 (MQ=255)
      ccGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATg                                                                                                <  2:1542171/99‑1 (MQ=255)
       cGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGc                                                                                              <  1:3534716/100‑1 (MQ=255)
         cTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATTTCCCATACCTGCGAGGCGATACGCGAGCGGGGCAGCACTTCGCCCTGATGGCga                                                                                            >  2:2669487/1‑100 (MQ=255)
         cTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCga                                                                                            <  2:2148506/100‑1 (MQ=255)
          ttCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaa                                                                                           <  1:1872934/100‑1 (MQ=255)
           tCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaag                                                                                          <  1:1052848/100‑1 (MQ=255)
           tCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaag                                                                                          <  2:1825274/100‑1 (MQ=255)
           tCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaa                                                                                           <  2:3568826/99‑1 (MQ=255)
            cACCGCCACATCAATAGCATTGGTATCGCTGTCAGAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagt                                                                                         <  1:1731164/100‑2 (MQ=255)
            cACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaagg                                                                                         >  2:1773870/1‑100 (MQ=255)
            cACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagt                                                                                         <  1:2768920/100‑2 (MQ=255)
            cACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaag                                                                                          <  1:1325066/99‑1 (MQ=255)
             aCCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagtg                                                                                        <  2:1481211/100‑3 (MQ=255)
              ccGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagtgc                                                                                       <  1:854711/100‑4 (MQ=255)
              ccGCCACATAAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACATGCGAGGCGATAAGCGAGCGGGCCAGCACTTCGCCCTGATGGCgaagga                                                                                        >  2:1845018/1‑99 (MQ=255)
               cGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagtgct                                                                                      <  2:27454/100‑5 (MQ=255)
                 ccACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagtgctgc                                                                                    <  1:2184924/100‑7 (MQ=255)
                 ccACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCATCAGGGCgaagtgctgc                                                                                    <  2:2181987/100‑7 (MQ=255)
                 ccACATCAATAGCATTGGTATCGCTGGTAAAAGTGAGGTCCCATACCTGCGAGGCGATAAGCGAGCGCGGCAGGACTTCGCCATCAGGGCgaagtg                                                                                        <  2:1497216/96‑3 (MQ=255)
                  cacaTCGATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAACACTTCGCCATCAGGGCgaagtgc                                                                                       <  2:1371292/96‑4 (MQ=255)
                  cacaTCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAAc                                                                                   <  1:2979227/100‑1 (MQ=255)
                    caTCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTc                                                                                 >  2:3116687/1‑100 (MQ=255)
                      tCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCca                                                                               <  1:2859069/100‑1 (MQ=255)
                       cAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcag                                                                              >  1:2211931/1‑100 (MQ=255)
                       cAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcag                                                                              >  2:73807/1‑100 (MQ=255)
                       cAATAGCACTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGCAGGAAGAACTCcag                                                                              >  2:7502/1‑100 (MQ=255)
                         aTAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagca                                                                            >  2:1430095/1‑100 (MQ=255)
                         aTAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagca                                                                            >  2:3490648/1‑100 (MQ=255)
                           aGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagcag                                                                           <  2:3174611/99‑1 (MQ=255)
                            gCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGGTGGCGAAGGAAGAACTCCAGCagag                                                                         >  2:3248977/1‑100 (MQ=255)
                            gCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCagag                                                                         >  1:105943/1‑100 (MQ=255)
                              aTTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTa                                                                       >  1:1528896/1‑100 (MQ=255)
                                 ggTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAAc                                                                    <  1:2405210/100‑1 (MQ=255)
                                  gTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAAct                                                                   <  1:2539956/100‑1 (MQ=255)
                                     tCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCtt                                                                <  1:2887711/100‑1 (MQ=255)
                                     tCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCtt                                                                <  1:818180/100‑1 (MQ=255)
                                        cTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTAc                                                             <  2:1925789/100‑1 (MQ=255)
                                        cTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTAc                                                             <  1:1823346/100‑1 (MQ=255)
                                         tGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACt                                                            >  2:892556/1‑100 (MQ=255)
                                          gTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTg                                                           >  2:2527020/1‑100 (MQ=255)
                                          gTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTg                                                           >  1:2388371/1‑100 (MQ=255)
                                            cAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGt                                                         >  2:1884002/1‑100 (MQ=255)
                                            cAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGt                                                         >  2:1039790/1‑100 (MQ=255)
                                              aaaTTCATGTCCCATACCTGCGAGGCGCTAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAGGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCa                                                       >  2:894629/1‑100 (MQ=255)
                                                aTTCATGTCCCATACCTGCGAGGCGATAAGCTAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCaaa                                                     <  1:1525218/100‑1 (MQ=255)
                                                   cATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGt                                                   <  2:661510/99‑1 (MQ=255)
                                                   cATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTg                                                  >  1:1200553/1‑100 (MQ=255)
                                                      gTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATg                                               <  1:1021563/100‑1 (MQ=255)
                                                      gTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATg                                               >  2:2125251/1‑100 (MQ=255)
                                                       tCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATgc                                              >  1:1185071/1‑100 (MQ=255)
                                                        cccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATgcg                                             <  1:3289403/100‑1 (MQ=255)
                                                        cccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATgcg                                             <  2:2696838/100‑1 (MQ=255)
                                                          cATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGcgcg                                           >  2:2312397/1‑100 (MQ=255)
                                                          cATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGcgcg                                           <  1:2661815/100‑1 (MQ=255)
                                                            tACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTg                                         >  1:668262/1‑100 (MQ=255)
                                                            tACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTg                                         <  2:487068/100‑1 (MQ=255)
                                                            tACCTCCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTg                                         <  2:1671502/100‑1 (MQ=255)
                                                               cTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCg                                      >  2:1595617/1‑100 (MQ=255)
                                                               cTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCg                                      >  1:735377/1‑100 (MQ=255)
                                                                  cGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTg                                   >  1:3210436/1‑100 (MQ=255)
                                                                  cGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTg                                   <  2:2473973/100‑1 (MQ=255)
                                                                   gAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGc                                  <  2:1729986/100‑1 (MQ=255)
                                                                   gAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGc                                  <  2:3252565/100‑1 (MQ=255)
                                                                     ggCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTGACTGGTCAAAGTGATGCGCGTGCCGCTGCgg                                >  1:2813435/1‑100 (MQ=255)
                                                                      gCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCggg                               <  2:534561/100‑1 (MQ=255)
                                                                      gCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCACGTGCCGCTGCggg                               <  2:2127396/100‑1 (MQ=255)
                                                                          tAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGAc                           <  2:1181910/100‑1 (MQ=255)
                                                                           nagcgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACt                          >  2:2186936/2‑100 (MQ=255)
                                                                            agcgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACtt                         >  1:1338929/1‑100 (MQ=255)
                                                                             gcgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACttt                        >  2:114402/1‑100 (MQ=255)
                                                                             gcgagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACttt                        >  2:3127591/1‑100 (MQ=255)
                                                                               gagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCTAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGc                      <  2:2683793/100‑1 (MQ=255)
                                                                               gagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGc                      <  1:653964/100‑1 (MQ=255)
                                                                                agcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCg                     <  2:2042457/100‑1 (MQ=255)
                                                                                agcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCg                     <  1:73807/100‑1 (MQ=255)
                                                                                 gcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCgg                    >  1:576064/1‑100 (MQ=255)
                                                                                  cgGGGCAGCACTTCGNCCTGNTGGCGANGGAAGAANTCCAGCAGAGTAANCTCNTTACTGGTCAAAGTGATGCGCGTGNCGCNGCGGGTGACTTTGCGGc                   >  2:2391947/1‑100 (MQ=255)
                                                                                  cgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTCAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGCCTTCGCGGc                   >  1:1396462/1‑100 (MQ=255)
                                                                                   ggggCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCt                  <  2:14702/100‑1 (MQ=255)
                                                                                    gggCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTg                 <  1:1626439/100‑1 (MQ=255)
                                                                                    gggCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTg                 <  2:3387393/100‑1 (MQ=255)
                                                                                    gggCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTg                 >  2:3161822/1‑100 (MQ=255)
                                                                                    gggCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTg                 >  1:3608678/1‑100 (MQ=255)
                                                                                     ggCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGa                >  2:611155/1‑100 (MQ=255)
                                                                                      ttcttccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTg                 >  2:3522917/7‑98 (MQ=38)
                                                                                      gcagcaCTTCTCCCTGATGGCGAAGGAAGAACGCCAGCAGAGTAAACTCTTTACTGCTCAAAGTGATGCGCGTGCCTCTGCGGGTGACTTTGCAGc                   <  2:3392199/96‑1 (MQ=255)
                                                                                        cttccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACga             >  1:2957543/5‑100 (MQ=255)
                                                                                         ttccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACgag            >  1:3193825/4‑100 (MQ=255)
                                                                                         gcaCTTCGCCCTGATGGCGACGGAAGAACTCCAGCAGAGTAACCTCTTTACTGGTCAAAGTGCTGCGCGTGCCGCTGCGGGCGACTCTGCGGc                   >  2:546778/1‑93 (MQ=255)
                                                                                          tccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGTGTGACTTTGCGGCTGACgaga           >  1:943726/3‑100 (MQ=255)
                                                                                             ttCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCg        >  1:676183/1‑100 (MQ=255)
                                                                                               cGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGTGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGAc      >  2:2776219/1‑100 (MQ=255)
                                                                                               cGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGAc      >  2:2141072/1‑100 (MQ=255)
                                                                                                gCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGAcc     <  1:141412/100‑1 (MQ=255)
                                                                                                 cccTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCa    >  1:2434032/1‑100 (MQ=255)
                                                                                                  ccTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCa    <  2:1814126/99‑1 (MQ=255)
                                                                                                   cTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATc  <  1:2305976/100‑1 (MQ=255)
                                                                                                   cTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATc  <  1:2305910/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATC  >  NC_000913/594847‑595045

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