Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 594,948 G→C 15.1% Q166E (CAG→GAG) ‡ 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,9480GC15.1% 400.2 / 54.8 145Q166E (CAG→GAG) cusRresponse regulator in two‑component regulatory system with CusS
Reads supporting (aligned to +/- strand):  ref base G (58/64);  new base C (5/17);  total (63/82)
Fisher's exact test for biased strand distribution p-value = 3.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

CAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAA  >  NC_000913/594850‑595047
                                                                                                   |                                                                                                   
  gCCGCTTTACCGCCACATCAATATCATGGGTATCGCTGTCAAAATGCATGTCCCATACCTGCGTGGCGATAAGCGATC‑GGTGCAGCACGTCTCCCTGATg                                                                                                  >  2:3273412/1‑100 (MQ=255)
  gCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATg                                                                                                  <  1:1716620/100‑1 (MQ=255)
   ccGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATgg                                                                                                 <  1:2315103/100‑1 (MQ=255)
   ccGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATg                                                                                                  <  2:1542171/99‑1 (MQ=255)
      cTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATTTCCCATACCTGCGAGGCGATACGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCga                                                                                              >  2:2669487/1‑100 (MQ=255)
      cTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCga                                                                                              >  2:1872925/1‑100 (MQ=255)
       ttCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCgaa                                                                                             <  1:652945/100‑1 (MQ=255)
       ttCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCgaa                                                                                             <  1:1872934/100‑1 (MQ=255)
       ttCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCgaa                                                                                             <  1:2843367/100‑1 (MQ=255)
        tCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCgaag                                                                                            <  2:1825274/100‑1 (MQ=255)
        tCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCgaag                                                                                            <  1:1052848/100‑1 (MQ=255)
        tCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCgaa                                                                                             <  2:3568826/99‑1 (MQ=255)
         cACCGCCACATCAATAGCATTGGTATCGCTGTCAGAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCATCAGGGCgaagt                                                                                           <  1:1731164/100‑2 (MQ=255)
         cACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCgaagg                                                                                           >  2:1773870/1‑100 (MQ=255)
         cACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCATCAGGGCgaagt                                                                                           <  2:2395815/100‑2 (MQ=255)
         cACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCATCAGGGCgaagt                                                                                           <  1:2768920/100‑2 (MQ=255)
          aCCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCATCAGGGCgaagtg                                                                                          <  2:1481211/100‑3 (MQ=255)
          aCCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCATCAGGGCgaagt                                                                                           <  2:1437025/99‑2 (MQ=255)
           ccGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGGCAGCACTTCGCCATCAGGGCgaagtg                                                                                          <  1:724397/100‑3 (MQ=255)
           ccGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCgaaggaa                                                                                         >  1:3252589/1‑100 (MQ=255)
           ccGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCATCAGGGCgaagtgc                                                                                         <  1:3582855/100‑4 (MQ=255)
            cGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCgaaggaag                                                                                        <  1:1427869/100‑1 (MQ=255)
            cGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCATCAGGGCgaagtgct                                                                                        <  2:27454/100‑5 (MQ=255)
            cGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCATCAGGGCgaagtgct                                                                                        <  2:3237657/100‑5 (MQ=255)
             gCCACTTCAATAGCATTGGTAGCGCTGAATAAATTCAGTTGCCATACCTGCGAGGCGATAAGCGAGT‑GGGGCAGCACTTCGCCCTGATGGCGAAGgaaga                                                                                       <  2:2189393/100‑1 (MQ=255)
              ccACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCATCAGGGCgaagtgctgc                                                                                      <  2:2181987/100‑7 (MQ=255)
               cacaTCGATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAACACTTCGCCATCAGGGCgaagtgc                                                                                         <  2:1371292/96‑4 (MQ=255)
               cacaTCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAAc                                                                                     <  1:2979227/100‑1 (MQ=255)
               cacaTCAATAGCATTGATATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCATCAGGGCgaagtgctgcc                                                                                     <  2:2133913/100‑8 (MQ=38)
                acaTCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACt                                                                                    >  2:1823337/1‑100 (MQ=255)
                   tCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGTGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGCAGGAAGAACTCca                                                                                 >  1:1611036/1‑100 (MQ=255)
                    cAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcag                                                                                >  1:2211931/1‑100 (MQ=255)
                    cAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcag                                                                                >  2:73807/1‑100 (MQ=255)
                    cAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcag                                                                                <  2:165411/100‑1 (MQ=255)
                    cAATAGCACTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGCAGGAAGAACTCcag                                                                                >  2:7502/1‑100 (MQ=255)
                      aTAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagca                                                                              >  2:3490648/1‑100 (MQ=255)
                      aTAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagca                                                                              <  1:1371726/100‑1 (MQ=255)
                        aGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagcag                                                                             <  2:3174611/99‑1 (MQ=255)
                         gCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCagag                                                                           >  1:105943/1‑100 (MQ=255)
                           aTTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTa                                                                         >  1:1528896/1‑100 (MQ=255)
                            tttgTATCGCAGTCAAAATGAATGTACCATACCTGCGAAGTGATAAGCTAGT‑GGTGCAGCACTTCGCCCTAATGGCAAAGGAAGATCTCcagc                                                                               <  2:1274369/90‑1 (MQ=255)
                              ggTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAAc                                                                      <  1:2405210/100‑1 (MQ=255)
                               gTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAAct                                                                     <  1:78084/100‑1 (MQ=255)
                               gTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAAct                                                                     <  1:2539956/100‑1 (MQ=255)
                                  tCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCtt                                                                  <  1:818180/100‑1 (MQ=255)
                                  tCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCtt                                                                  <  1:2887711/100‑1 (MQ=255)
                                     cTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTAc                                                               <  2:1925789/100‑1 (MQ=255)
                                      tGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACt                                                              >  2:892556/1‑100 (MQ=255)
                                        tCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTgg                                                            <  1:1892301/100‑1 (MQ=255)
                                         cAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGt                                                           <  2:2148624/100‑1 (MQ=255)
                                         cAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGt                                                           >  2:1884002/1‑100 (MQ=255)
                                         cAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGt                                                           <  2:1113026/100‑1 (MQ=255)
                                          aaaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTc                                                          >  1:917798/1‑100 (MQ=255)
                                          aaaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTc                                                          >  2:2566694/1‑100 (MQ=255)
                                           aaaTTCATGTCCCATACCTGCGAGGCGCTAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAGGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCa                                                         >  2:894629/1‑100 (MQ=255)
                                           aaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCa                                                         <  2:1741109/100‑1 (MQ=255)
                                                cATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGt                                                     <  2:661510/99‑1 (MQ=255)
                                                cATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTg                                                    >  1:1200553/1‑100 (MQ=255)
                                                 aTGTCCCATGCCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGa                                                   <  2:2058705/100‑1 (MQ=255)
                                                 aTGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGa                                                   <  2:2712229/100‑1 (MQ=255)
                                                   gTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATg                                                 <  1:1021563/100‑1 (MQ=255)
                                                   gTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATg                                                 >  2:2125251/1‑100 (MQ=255)
                                                     cccATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATgcg                                               <  1:3289403/100‑1 (MQ=255)
                                                       cATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGcgcg                                             >  2:2312397/1‑100 (MQ=255)
                                                       cATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGcgcg                                             <  1:1776942/100‑1 (MQ=255)
                                                         tACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTg                                           <  2:487068/100‑1 (MQ=255)
                                                         tACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTg                                           >  1:668262/1‑100 (MQ=255)
                                                         tACCTCCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTg                                           <  2:1671502/100‑1 (MQ=255)
                                                           ccTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGcc                                         <  1:2592717/100‑1 (MQ=255)
                                                            cTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCg                                        >  1:735377/1‑100 (MQ=255)
                                                             tGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGc                                       >  1:2414457/1‑100 (MQ=255)
                                                              gCGAGGCGATAAGCGAGC‑GGTGCAGCACTACGCCCTGATGGCGAAGGGAGAACTCCAGCAGAGTAAACTCTCTACTGGTCAAAGTGATGCGCGTCCCGCt                                      >  2:2791108/1‑100 (MQ=255)
                                                               cGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTg                                     <  2:2473973/100‑1 (MQ=255)
                                                               cGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTg                                     >  1:3210436/1‑100 (MQ=255)
                                                                gAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGc                                    <  2:3252565/100‑1 (MQ=255)
                                                                  ggCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTGACTGGTCAAAGTGATGCGCGTGCCGCTGCgg                                  >  1:2813435/1‑100 (MQ=255)
                                                                   gCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCggg                                 <  2:534561/100‑1 (MQ=255)
                                                                     gATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTg                               <  1:3028885/100‑1 (MQ=255)
                                                                       tAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGAc                             >  2:1504150/1‑100 (MQ=255)
                                                                         agcgagc‑gGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACtt                           >  1:1338929/1‑100 (MQ=255)
                                                                          gcgagc‑gGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACttt                          >  2:3127591/1‑100 (MQ=255)
                                                                          gcgagc‑gGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACttt                          >  2:114402/1‑100 (MQ=255)
                                                                            gagc‑gGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCTAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGc                        <  2:2683793/100‑1 (MQ=255)
                                                                            gagc‑gGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGc                        >  1:2235129/1‑100 (MQ=255)
                                                                            gagc‑gGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGc                        <  1:653964/100‑1 (MQ=255)
                                                                               c‑gGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCgg                      <  1:2892182/99‑1 (MQ=255)
                                                                                 ggggCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCt                    <  2:14702/100‑1 (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)
                                                                                       ttccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACgag              >  1:1263102/4‑100 (MQ=255)
                                                                                       ttccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACgag              >  1:3193825/4‑100 (MQ=255)
                                                                                       gcaCTTCGCCCTGATGGCGACGGAAGAACTCCAGCAGAGTAACCTCTTTACTGGTCAAAGTGCTGCGCGTGCCGCTGCGGGCGACTCTGCGGc                     >  2:546778/1‑93 (MQ=255)
                                                                                        tccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGTGTGACTTTGCGGCTGACgaga             >  1:943726/3‑100 (MQ=255)
                                                                                             cGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGAc        >  2:2141072/1‑100 (MQ=255)
                                                                                              gCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGAcc       <  1:141412/100‑1 (MQ=255)
                                                                                                 cTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATc    <  1:2305976/100‑1 (MQ=255)
                                                                                                 cTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGACCGACCAt     >  2:3148444/1‑99 (MQ=255)
                                                                                                  tGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCa   <  2:31096/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGC‑GGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAA  >  NC_000913/594850‑595047

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