Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,031,761 T→C 47.9% intergenic (‑49/‑378) serT ← / → hyaA tRNA‑Ser/hydrogenase 1, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,031,7610TC47.9% ‑4.3 / 127.8 90intergenic (‑49/‑378)serT/hyaAtRNA‑Ser/hydrogenase 1, small subunit
Reads supporting (aligned to +/- strand):  ref base T (19/28);  new base C (19/24);  total (38/52)
Fisher's exact test for biased strand distribution p-value = 8.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GGGTCGCCGGTTTTCAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTTCTCTTTTAC  >  NC_000913/1031663‑1031856
                                                                                                  |                                                                                               
gggTCGCCGGTTTTCAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAAtt                                                                                                >  2:1517503/1‑100 (MQ=255)
    cGCCGGTTTTCAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAAtttgtt                                                                                            <  2:925840/100‑1 (MQ=255)
     gCCGGTTTTCAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAAtttgttt                                                                                           <  1:883908/100‑1 (MQ=255)
     gCCGGTTTTCAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAAtttgttt                                                                                           <  1:899288/100‑1 (MQ=255)
     gCCGGTTTTCAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAActtgtt                                                                                            <  2:1598693/99‑1 (MQ=255)
      ccGGTTTTCAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACTCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAAtttgtttg                                                                                          <  2:1321902/100‑1 (MQ=255)
         gTTTTCAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAAtttgtttgttt                                                                                       >  1:133884/1‑100 (MQ=255)
         gTTTTCAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAActtgtttgttt                                                                                       <  1:53636/100‑1 (MQ=255)
          ttttCAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAAtttgtttgttt                                                                                       <  2:1730828/99‑1 (MQ=255)
          ttttCAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAActtgtttgtttg                                                                                      >  1:1209533/1‑100 (MQ=255)
            ttCAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGGCGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGcc                                                                                    <  2:1779654/100‑1 (MQ=255)
              cAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCCATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTg                                                                                  <  2:1273972/100‑1 (MQ=255)
                aGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGaa                                                                                >  1:1899969/1‑100 (MQ=255)
                aGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGaa                                                                                >  2:548335/1‑100 (MQ=255)
                 gACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGaaa                                                                               >  1:398522/1‑100 (MQ=255)
                  aCCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGaaaa                                                                              >  1:1809136/1‑100 (MQ=255)
                  aCCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGaaaa                                                                              <  2:418489/100‑1 (MQ=255)
                  aCCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGaaaa                                                                              <  1:1324246/100‑1 (MQ=255)
                   ccGGTGCCTTCAAGCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGaaaaa                                                                             >  2:1742200/1‑100 (MQ=255)
                     ggTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACa                                                                           >  1:110379/1‑100 (MQ=255)
                       tGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGc                                                                         <  2:1157008/100‑1 (MQ=255)
                        gCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGcc                                                                        >  1:1003758/1‑100 (MQ=255)
                         ccTTCAACCGCTCGGCCACACTTCCGGAATGACTCGCACTACAAACATCCCTATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCa                                                                       <  2:983729/100‑1 (MQ=255)
                          cTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCaa                                                                      >  1:461550/1‑100 (MQ=255)
                          cTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCaa                                                                      >  1:845498/1‑100 (MQ=255)
                           ttCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCaaa                                                                     >  1:1832748/1‑100 (MQ=255)
                            tCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCaaaa                                                                    >  1:281268/1‑100 (MQ=255)
                             cAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAg                                                                   >  2:233555/1‑100 (MQ=255)
                              aaCCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGt                                                                  <  1:2016888/100‑1 (MQ=255)
                               aCCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAg                                                                   >  1:2030143/1‑98 (MQ=255)
                               aCCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTg                                                                 >  1:1970434/1‑100 (MQ=255)
                                  gCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCAt                                                              <  1:1335730/100‑1 (MQ=255)
                                  gCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCAt                                                              <  2:586920/100‑1 (MQ=255)
                                   cTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCAtt                                                             <  2:1003747/100‑1 (MQ=255)
                                   cTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCAtt                                                             <  2:853830/100‑1 (MQ=255)
                                   cTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCAtt                                                             >  1:856094/1‑100 (MQ=255)
                                    tCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTg                                                            <  1:191074/100‑1 (MQ=255)
                                         cacaCTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAg                                                        >  1:2186106/1‑99 (MQ=255)
                                            acTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTg                                                    <  2:133881/100‑1 (MQ=255)
                                            acTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTCTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTg                                                    <  2:1683117/100‑1 (MQ=255)
                                             cTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTg                                                    <  2:1512924/99‑1 (MQ=255)
                                             cTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGa                                                   >  2:482040/1‑100 (MQ=255)
                                             cTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGa                                                   <  2:1461329/100‑1 (MQ=255)
                                              ttCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGaa                                                  <  1:268133/100‑1 (MQ=255)
                                              ctCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGaa                                                  <  1:1954972/99‑1 (MQ=255)
                                               tCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGaaa                                                 <  1:2096341/100‑1 (MQ=255)
                                                ccGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGaaag                                                >  2:1954950/1‑99 (MQ=255)
                                                ccGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAAt                                                >  2:1878210/1‑100 (MQ=255)
                                                ccGGAATGACGCGCACTATAAACATCCCGATGCAGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAAt                                                <  2:1381801/100‑1 (MQ=255)
                                                 cGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATa                                               <  2:1725974/100‑1 (MQ=255)
                                                 cGGAATGACGCGCACTATAAACATCCCGATGCCGAGTGTAAACCTCTAATTTGTTTGTTTGCCTGAAATACAGCCAAAAGTGCATTGATAGCGTGAAAt                                                <  2:1393707/99‑1 (MQ=255)
                                                  ggAATGACGCGCACTATAAACATCCAGATGCGGCGTGTAAACACCTAATTTGTTTGTTTGCCTGAAAAAAAGACAAAAGTGCATTGATAGCGTGAAATaa                                              >  1:1182537/1‑100 (MQ=255)
                                                   gAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCGTTGATAGCGTGAAATAAc                                             <  2:839601/100‑1 (MQ=255)
                                                    aaTGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAAca                                            >  1:2118080/1‑100 (MQ=255)
                                                    aaTGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAAca                                            <  2:834735/100‑1 (MQ=255)
                                                    aaTGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAAca                                            <  2:1345693/100‑1 (MQ=255)
                                                    aaTGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAAca                                            <  1:233558/100‑1 (MQ=255)
                                                    aaTGACGCGCACTATAAACAGCCCGATGCGGCGTGTAAGCCCCTTACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAAca                                            <  2:466300/100‑1 (MQ=255)
                                                       gACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAAcagca                                         <  2:506938/100‑1 (MQ=255)
                                                       gACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAAcagca                                         <  1:2348750/100‑1 (MQ=255)
                                                         cgcgCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGa                                       >  1:2224306/1‑100 (MQ=255)
                                                         cgcgCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAACAACAGCCAAAAATGCATTGATAGCGTGAAATAACAGCAGa                                       >  1:2124309/1‑100 (MQ=255)
                                                         cgcgCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGa                                       >  2:162544/1‑100 (MQ=255)
                                                           cgcACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGAtt                                     <  2:2295036/100‑1 (MQ=255)
                                                           cgcACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGAtt                                     >  1:142554/1‑100 (MQ=255)
                                                             cACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGa                                   >  1:2217983/1‑100 (MQ=255)
                                                             cACTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGa                                   <  2:2106017/100‑1 (MQ=255)
                                                               cTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGAt                                  <  2:1626209/99‑1 (MQ=255)
                                                               cTATAAACATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGATTGAATGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATc                                 <  2:414291/100‑1 (MQ=255)
                                                                tataAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCa                                <  1:1130505/100‑1 (MQ=255)
                                                                 ataAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCAt                               >  2:2192384/1‑100 (MQ=255)
                                                                  taAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCAtt                              <  1:68917/100‑1 (MQ=255)
                                                                   aaaCATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGAGTGAAATAACAGCAGATTGATCAttt                             <  1:1453033/100‑1 (MQ=255)
                                                                    aaCATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTtc                            >  2:1314407/1‑100 (MQ=255)
                                                                    aaCATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCAttt                             >  1:952628/1‑99 (MQ=255)
                                                                    aaCATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTtc                            <  2:2093689/100‑1 (MQ=255)
                                                                     aCATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTtca                           <  1:1550634/100‑1 (MQ=255)
                                                                      cATCCCGATGCGGCGTGTAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTtcat                          <  1:767076/100‑1 (MQ=255)
                                                                          ccGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCAcc                      >  2:1093610/1‑100 (MQ=255)
                                                                             aTGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATg                   <  2:2118056/100‑1 (MQ=255)
                                                                                 ggcgTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAAtt               >  2:53155/1‑100 (MQ=255)
                                                                                    gtgtAAGCCCCTAACTTGTTTGTTTGCCTTAAAAACAGCCATAAGTGCATTGATAGCGTGAAATATCAGCAGATTGATCATTTCATCACCATGAAtt               <  2:1426005/97‑1 (MQ=255)
                                                                                    gtgtAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCt            <  2:2405235/100‑1 (MQ=255)
                                                                                      gtAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTtc          <  2:375409/100‑1 (MQ=255)
                                                                                      gtAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTtc          <  2:1768824/100‑1 (MQ=255)
                                                                                      gtAAACCCCTAATGTGTTTGTTTGTCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTGATCACCATGAATTCCTtc          <  2:2041394/100‑1 (MQ=255)
                                                                                      gtAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATCTCATCACCATGAATTCCTtc          >  2:1873856/1‑100 (MQ=255)
                                                                                       tAAACCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTtct         <  1:794861/100‑1 (MQ=255)
                                                                                        aaaCCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTtctc        >  2:2173346/1‑100 (MQ=255)
                                                                                        aaaCCCCTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTtctc        <  1:555939/100‑1 (MQ=255)
                                                                                           ccccTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTTCTCttt     >  2:1946825/1‑100 (MQ=255)
                                                                                           ccccTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTTCTCttt     >  2:658578/1‑100 (MQ=255)
                                                                                           ccccTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTTCTCttt     >  2:626074/1‑100 (MQ=255)
                                                                                            cccTAATTTGTTTGTTTGCATGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCATATTGATCATTTCATCACCATGAATTCCTTCTCtttt    <  2:1598138/100‑1 (MQ=255)
                                                                                             ccTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTTCTCTTTTa   >  2:470282/1‑100 (MQ=255)
                                                                                             ccTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTTCTCTTTTa   <  2:1444949/100‑1 (MQ=255)
                                                                                              cTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTTCTCTTTTAc  <  2:398519/100‑1 (MQ=255)
                                                                                              cTAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTTCTCTTTTAc  <  1:1812554/100‑1 (MQ=255)
                                                                                               tAACTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTTCTCTTTTa   >  1:145817/1‑98 (MQ=255)
                                                                                                 aTTTGTTTGTTTGCATGAATAACAGCCAAATGTGCATTGATAGCGTGAAATTACATCATATTGATCATTTGATTACCATGAATTCCTTCTCttt     <  2:1171463/94‑1 (MQ=255)
                                                                                                  |                                                                                               
GGGTCGCCGGTTTTCAAGACCGGTGCCTTCAACCGCTCGGCCACACTTCCGGAATGACGCGCACTATAAACATCCCGATGCGGCGTGTAAACCCCTAATTTGTTTGTTTGCCTGAAAAACAGCCAAAAGTGCATTGATAGCGTGAAATAACAGCAGATTGATCATTTCATCACCATGAATTCCTTCTCTTTTAC  >  NC_000913/1031663‑1031856

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: