Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,962,561 G→A 15.7% N68N (AAC→AAT) ‡ ppdC ← putative prepilin peptidase‑dependent protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,962,5610GA15.7% 284.1 / 39.0 121N68N (AAC→AATppdCputative prepilin peptidase‑dependent protein
Reads supporting (aligned to +/- strand):  ref base G (49/53);  new base A (8/11);  total (57/64)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

Only 100 of 123 total aligned reads displayed.

CAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTT  >  NC_000913/2962462‑2962660
                                                                                                   |                                                                                                   
cAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAg                                                                                                     >  2:2142575/1‑100 (MQ=255)
  ggCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGtt                                                                                                   <  1:1306368/100‑1 (MQ=255)
      ggTCATCTCGCCTTCTCTGCCCCCGGGTGAATCTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCa                                                                                               <  1:1536138/100‑1 (MQ=255)
      ggTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGCCACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCa                                                                                               >  1:1555459/1‑100 (MQ=255)
        tCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCTAGTTGGCAgg                                                                                             >  1:1695360/1‑100 (MQ=255)
         cATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAgg                                                                                             >  1:2393506/1‑99 (MQ=255)
          aTCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTg                                                                                           <  2:2047712/100‑1 (MQ=255)
          aTCTCGCCTTCCCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCAGTCGGTTGACCTGCCAGTTGGCAGGTg                                                                                           <  2:1580080/100‑1 (MQ=255)
           tctcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTgg                                                                                          <  1:837787/100‑1 (MQ=255)
           tctcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTgg                                                                                          >  2:406711/1‑100 (MQ=255)
           tctcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTgg                                                                                          <  2:618757/100‑1 (MQ=255)
            ctcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGc                                                                                         >  2:1252523/1‑100 (MQ=255)
            ctcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGc                                                                                         >  2:2512089/1‑100 (MQ=255)
             tcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCg                                                                                        >  2:2004117/1‑100 (MQ=255)
             tcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtcaac                                                                                        <  2:570243/100‑5 (MQ=255)
             tcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtcaac                                                                                        <  1:2213389/100‑5 (MQ=255)
             tcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtcaac                                                                                        <  1:1140100/100‑5 (MQ=255)
             tcGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGCTGACCTGCCAACTGGCAGGtcaac                                                                                        <  2:1693088/100‑5 (MQ=255)
              cGCCTTCTCTGCCCCCGGGTGACACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTTACCTGCCAACTGGCAGGtcaacc                                                                                       <  1:1772883/100‑6 (MQ=255)
              cGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGgaga                                                                                       >  1:2425753/1‑100 (MQ=255)
              cGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGa                                                                                       <  1:646016/100‑1 (MQ=255)
              cGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtcaac                                                                                        <  2:1445255/99‑5 (MQ=255)
              cGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtcaac                                                                                        <  2:1110510/99‑5 (MQ=255)
               gCCTTCTCTGCCCCCGGGTGAAACTATCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGaa                                                                                      >  1:1158994/1‑100 (MQ=255)
               gCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGaa                                                                                      <  1:1443112/100‑1 (MQ=255)
               gCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGaa                                                                                      >  2:359039/1‑100 (MQ=255)
                 cTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAAt                                                                                    >  2:2412863/1‑100 (MQ=255)
                 cTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAAt                                                                                    <  2:305836/100‑1 (MQ=255)
                 cTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtcaaccaga                                                                                    >  1:1177926/1‑92 (MQ=255)
                  ttCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAACTGGCAGGtcaaccgaat                                                                                   <  1:1154881/100‑10 (MQ=38)
                      ctGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgcgc                                                                               >  2:1420648/1‑100 (MQ=255)
                      ctGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgcgc                                                                               <  1:1252519/100‑1 (MQ=255)
                        gcccCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGCTGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATcgcgcg                                                                              >  1:487732‑M2/5‑99 (MQ=255)
                         cccCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCa                                                                            >  1:525582‑M2/4‑100 (MQ=255)
                          ccCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAg                                                                           >  1:1641878‑M2/3‑100 (MQ=255)
                               ggTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCg                                                                      >  2:907894/1‑100 (MQ=255)
                                   aaaCTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGNAATNGCGCGCAGTTGCGTTtg                                                                  >  1:1531011/1‑100 (MQ=255)
                                   aaaCTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtg                                                                  >  2:2228254/1‑100 (MQ=255)
                                      cTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgctg                                                               <  2:1249480/100‑1 (MQ=255)
                                      cTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTtgctg                                                               <  2:2225774/100‑1 (MQ=255)
                                        aGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgcc                                                             <  2:783650/100‑1 (MQ=255)
                                          cGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccag                                                           <  1:2282729/100‑1 (MQ=255)
                                          cGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccag                                                           >  1:1372676/1‑100 (MQ=255)
                                           gTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTgccagc                                                          >  2:294643/1‑100 (MQ=255)
                                                gCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATg                                                      >  1:1947725/1‑99 (MQ=255)
                                                   gATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATgccgc                                                  <  2:844295/100‑1 (MQ=255)
                                                    aTGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATgccgcc                                                 <  2:931827/100‑1 (MQ=255)
                                                     tGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCa                                                <  1:1472145/100‑1 (MQ=255)
                                                     tGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCa                                                >  1:1962891/1‑100 (MQ=255)
                                                      gCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCag                                               >  2:372291/1‑100 (MQ=255)
                                                       cTGACACATCCCGCCTGCGATGTCTGCATTCAGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCaga                                              >  2:1257588/1‑100 (MQ=255)
                                                        tGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCagag                                             >  2:1160541/1‑100 (MQ=255)
                                                        tGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCagag                                             <  1:1923558/100‑1 (MQ=255)
                                                        tGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCagag                                             >  1:1644431/1‑100 (MQ=255)
                                                        tGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTCGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCaga                                              >  1:1033668/1‑99 (MQ=255)
                                                         gACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGc                                            <  2:1589225/100‑1 (MQ=255)
                                                          acacATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCt                                           >  1:1086362/1‑100 (MQ=255)
                                                           cacaTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCCGTTGCGTTTGCTGCCAGCCATGCCGCCag                                               >  1:630719/1‑95 (MQ=255)
                                                           cacaTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTg                                          >  2:1050803/1‑100 (MQ=255)
                                                             caTCCCGCCTGCTATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGtt                                        >  2:1703013/1‑100 (MQ=255)
                                                             caTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTGTGCTGCCAGCCATGCGGCCAGAGCTGt                                         >  1:155487/1‑99 (MQ=255)
                                                             caTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGATGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGtt                                        >  2:1298127/1‑100 (MQ=255)
                                                              aTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGTCAGTTGGCAGGTGGCGAAATCGCGTGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTg                                       <  1:1933923/100‑1 (MQ=255)
                                                              aTCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTg                                       <  2:2441634/100‑1 (MQ=255)
                                                               tCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTgg                                      >  1:71171/1‑100 (MQ=255)
                                                                     cTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTgg                                <  2:2325586/100‑1 (MQ=255)
                                                                      tGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGt                               <  1:1759286/100‑1 (MQ=255)
                                                                       gCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGtt                              <  2:545032/100‑1 (MQ=255)
                                                                         gATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTc                            <  2:525584/100‑1 (MQ=255)
                                                                         gATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTc                            <  1:670142/100‑1 (MQ=255)
                                                                              cTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTctgct                        >  1:230858/1‑99 (MQ=255)
                                                                               tGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTgg                      <  2:1279889/100‑1 (MQ=255)
                                                                                gCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGc                     <  2:1695369/100‑1 (MQ=255)
                                                                                   tgcGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAg                  <  1:62359/98‑1 (MQ=255)
                                                                                     ttcgccaCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTCGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGAc                >  2:2531142/7‑100 (MQ=255)
                                                                                     cGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGAc                <  1:372289/100‑1 (MQ=255)
                                                                                     cGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGAc                <  2:1213804/100‑1 (MQ=255)
                                                                                      ggTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACt               <  1:2348408/100‑1 (MQ=255)
                                                                                       cgccaCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTg              >  2:1499309/5‑100 (MQ=255)
                                                                                       cgccaCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTACCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTg              >  1:176075/5‑100 (MQ=255)
                                                                                        ttGACCTGCCAGTTGGCAGGTAGAGAAATCGAGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGt             >  1:182120/1‑100 (MQ=255)
                                                                                        gccaCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGt             >  2:1708726/4‑100 (MQ=255)
                                                                                         ccaCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGtt            >  1:404499/3‑100 (MQ=255)
                                                                                         ccaCCTGCCAACTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGtt            >  1:604438/3‑100 (MQ=255)
                                                                                          gACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTc           >  2:1604667/1‑100 (MQ=255)
                                                                                          gACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTc           <  2:1995015/100‑1 (MQ=255)
                                                                                           aCCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCa          >  1:386712/1‑100 (MQ=255)
                                                                                           aCCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCa          <  1:1257584/100‑1 (MQ=255)
                                                                                           aCCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCa          <  1:1257604/100‑1 (MQ=255)
                                                                                           aCCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCa          <  1:2135465/100‑1 (MQ=255)
                                                                                           aCCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCa          <  2:2071653/100‑1 (MQ=255)
                                                                                            ccTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAt         >  1:2350398/1‑100 (MQ=255)
                                                                                            ccTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAt         >  1:160814/1‑100 (MQ=255)
                                                                                            ccTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTCGTTTCTGCTGGCAAGACTGTTCAt         >  1:949272/1‑100 (MQ=255)
                                                                                             cTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAtt        <  1:1863142/100‑1 (MQ=255)
                                                                                             cTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAtt        <  2:162907/100‑1 (MQ=255)
                                                                                             cTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCAtt        <  1:2508124/100‑1 (MQ=255)
                                                                                              tGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTGGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTCCTGGTTTCTGCTGGCAAGACCGTTCATTa       >  1:1552994/1‑100 (MQ=255)
                                                                                                 cAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTGATg    <  2:2370243/100‑1 (MQ=255)
                                                                                                   gTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCGGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGtt  <  1:880826/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CAGGCGGGTCATCTCGCCTTCTCTGCCCCCGGGTGAAACTAGCGTAACGCTGATGCTGACACATCCCGCCTGCGATGTCTGCATTCGGTTGACCTGCCAGTTGGCAGGTGGCGAAATCGCGCGCAGTTGCGTTTGCTGCCAGCCATGCCGCCAGAGCTGTTGGTACTGGTTTCTGCTGGCAAGACTGTTCATTAATGTT  >  NC_000913/2962462‑2962660

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: