Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 603,395 C→G 19.0% intergenic (+60/+21) pheP → / ← ybdG phenylalanine transporter/mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913603,3950CG19.0% 469.6 / 110.3 196intergenic (+60/+21)pheP/ybdGphenylalanine transporter/mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (86/72);  new base G (23/14);  total (110/86)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-03

Only 100 of 208 total aligned reads displayed.

ATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCG  >  NC_000913/603296‑603494
                                                                                                   |                                                                                                   
aTTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTc                                                                                                     >  1:1659272/1‑100 (MQ=255)
 ttCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg                                                                                                    <  2:2268023/100‑1 (MQ=255)
   ccTGTTTATGGCATTTAAAACGCTGCGTCGGACATAAGGCATTCACGCTACATCCGACAAAACGCTCTCAACCATCCGAAACCGCTCTCCTCCCTTCg                                                                                                    >  1:3814221/1‑98 (MQ=255)
   ccTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                                  >  2:1620802/1‑100 (MQ=255)
    cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                                 <  1:1122321/100‑1 (MQ=255)
    cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                                 <  1:3740200/100‑1 (MQ=255)
    cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                                 <  2:3076535/100‑1 (MQ=255)
     tGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATga                                                                                                <  2:1283439/100‑1 (MQ=255)
         tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                            <  2:4250978/100‑1 (MQ=39)
         tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGc                                                                                            <  2:1109509/100‑1 (MQ=255)
          aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                            <  1:2920606/99‑1 (MQ=38)
          aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCg                                                                                           <  2:270196/100‑1 (MQ=38)
           tGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCgg                                                                                          <  2:1799345/100‑1 (MQ=38)
           tGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                                                                                          >  1:3914416/1‑100 (MQ=255)
           tGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCg                                                                                           <  1:3032118/99‑1 (MQ=255)
            ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAATGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGt                                                                                         <  1:2441768/100‑1 (MQ=38)
            ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGt                                                                                         <  1:1880180/100‑1 (MQ=38)
            ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGt                                                                                         >  2:4218805/1‑100 (MQ=255)
             gCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGtt                                                                                        <  1:1681152/100‑1 (MQ=38)
              cATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttt                                                                                       >  2:1405457/1‑100 (MQ=255)
              cATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGt                                                                                         >  1:1755768/1‑98 (MQ=255)
               aTTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGtttc                                                                                      <  2:1895299/100‑2 (MQ=38)
                  tAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttttt                                                                                   >  2:2062030/1‑100 (MQ=255)
                  tAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttttt                                                                                   <  1:3394210/100‑1 (MQ=255)
                   aaaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTa                                                                                  <  1:585772/100‑1 (MQ=255)
                   aaaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTa                                                                                  <  2:559612/100‑1 (MQ=255)
                      aCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTAATCCATTCGATGAGAGCGGTCTTTTTAAtt                                                                               >  1:3896145/1‑100 (MQ=255)
                        gCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACATAACGATGTCAACCATCCGAGACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTAc                                                                             <  1:2393613/100‑1 (MQ=255)
                         cTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACt                                                                            >  2:2714858/1‑100 (MQ=255)
                           gCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTg                                                                           <  2:4355410/99‑1 (MQ=255)
                             gTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCtt                                                                        <  1:2260126/100‑1 (MQ=255)
                                ggACATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCCCTCATCCACTCGATGAGAGCGGCTTTTTTAATTACTGCTTcaa                                                                     >  1:3650033/1‑100 (MQ=255)
                                ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                     >  1:2655741/1‑100 (MQ=255)
                                ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                     <  1:143542/100‑1 (MQ=255)
                                     tAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAc                                                                >  2:886861/1‑100 (MQ=255)
                                      aaGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAcc                                                               <  1:1318848/100‑1 (MQ=255)
                                        ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAACCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                                   >  1:416913/1‑94 (MQ=255)
                                        ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCg                                                             >  1:800088/1‑100 (MQ=255)
                                        ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCg                                                             >  2:3913540/1‑100 (MQ=255)
                                          cATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGcc                                                           <  2:2173539/100‑1 (MQ=255)
                                          cATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGcc                                                           <  2:1497938/100‑1 (MQ=255)
                                              cACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagag                                                       >  1:888391/1‑100 (MQ=255)
                                              cACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagag                                                       >  1:3682667/1‑100 (MQ=255)
                                                cGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagagag                                                     >  2:3212131/1‑100 (MQ=255)
                                                cGCTACAGCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagaga                                                      <  2:4293762/99‑1 (MQ=255)
                                                  cTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCg                                                   <  2:4360776/100‑1 (MQ=255)
                                                   tACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGa                                                  <  2:2837627/100‑1 (MQ=255)
                                                    aCATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGTGCGa                                                  >  1:4293532/1‑99 (MQ=255)
                                                    aCATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGaa                                                 <  2:3049350/100‑1 (MQ=255)
                                                      aTCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAata                                               <  1:1593130/100‑1 (MQ=255)
                                                        ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTACTTACTGCTTAAATGCACCCGCCAGAGAGCGAATATc                                             >  1:2180833/1‑100 (MQ=255)
                                                        ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTATATGCACCCGCCAGAGAGCGAATATc                                             >  1:4312369/1‑100 (MQ=255)
                                                         cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCa                                            >  2:3118082/1‑100 (MQ=255)
                                                            cAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTg                                         >  2:1716434/1‑100 (MQ=255)
                                                             aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGc                                        >  2:1187041/1‑100 (MQ=255)
                                                             aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGGTTTTGTAATTACTGCTTAAATGCACCCGCCAGAGGGCGAATATCa                                            >  1:4344995/1‑96 (MQ=255)
                                                                aCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCgg                                     <  2:3843259/100‑1 (MQ=255)
                                                                  gATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                                   >  2:3402478/1‑100 (MQ=255)
                                                                    tGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                                   >  1:757721/1‑98 (MQ=255)
                                                                    tGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTgg                                 <  2:531715/100‑1 (MQ=255)
                                                                      tCAACCATCCGAAACCGCTCTCATGCATTCGATGAGAGCGGTTTCTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCg                               <  2:839589/100‑1 (MQ=255)
                                                                      tCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCg                               >  2:484385/1‑100 (MQ=255)
                                                                       ggaattaaaaaaaCCGCTCTCATCGAATGGACGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGa                              >  2:2964011/11‑100 (MQ=25)
                                                                       cAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGTTTGGCGa                              >  1:1253157/1‑100 (MQ=255)
                                                                       cAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTCTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGa                              >  2:2637891/1‑100 (MQ=255)
                                                                        aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGTTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGAc                             <  2:1672472/100‑1 (MQ=255)
                                                                           cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtga                          >  1:2044275/1‑100 (MQ=255)
                                                                            aTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgat                         >  2:2417650/1‑100 (MQ=255)
                                                                            aTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgat                         <  2:126429/100‑1 (MQ=255)
                                                                            aTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgat                         <  2:1935348/100‑1 (MQ=255)
                                                                             tCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatg                        >  2:1579500/1‑100 (MQ=255)
                                                                                 aaaCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGt                    <  2:3714724/100‑1 (MQ=255)
                                                                                  aaCCGCTCTCATCCATTCGATGAGGGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGtn                   <  2:1790085/100‑2 (MQ=255)
                                                                                   aCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCg                  >  2:2541394/1‑100 (MQ=255)
                                                                                    ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTATTTTCTTCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 >  1:410058/1‑100 (MQ=37)
                                                                                    ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAGCAGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 >  1:2880758/1‑100 (MQ=38)
                                                                                    ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 >  1:483265/1‑100 (MQ=38)
                                                                                    ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 >  1:1553397/1‑100 (MQ=38)
                                                                                    ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 >  2:2810236/1‑100 (MQ=38)
                                                                                    ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATACCATTGCCGGTTGGCGACTGATGAAGTCGc                 >  2:1305126/1‑100 (MQ=38)
                                                                                    ccGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 <  2:3298516/100‑1 (MQ=255)
                                                                                    ccGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 >  2:4002719/1‑100 (MQ=255)
                                                                                    ccGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 <  2:3470602/100‑1 (MQ=255)
                                                                                    ccGCTCTCATCCATTCGATGAGAGCGGTTGTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 <  2:3793100/100‑1 (MQ=255)
                                                                                     cGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                >  2:4075530/1‑100 (MQ=38)
                                                                                     cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                <  1:1008979/100‑1 (MQ=255)
                                                                                     cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                >  2:2381690/1‑100 (MQ=255)
                                                                                      gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTTTTTTTTAAATGCACCCGCCAGAGAGCGACTATCATTGCCGGTTGGCGACTGATGAAGTCTCAg               >  1:2742830/1‑100 (MQ=25)
                                                                                      gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               >  1:2988059/1‑100 (MQ=39)
                                                                                      gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               >  2:3587571/1‑100 (MQ=39)
                                                                                      gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               <  1:1736419/100‑1 (MQ=255)
                                                                                       ctctCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGa              >  1:3571783/1‑100 (MQ=39)
                                                                                             tCCATTCGATGAGAGCGGTTTTTTTATTTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc        <  2:441931/100‑1 (MQ=255)
                                                                                             tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc        <  1:2161441/100‑1 (MQ=255)
                                                                                             tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc        <  2:1784061/100‑1 (MQ=255)
                                                                                                aTTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAActct     >  2:928781/1‑100 (MQ=255)
                                                                                                  tCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTc   <  1:2322774/100‑1 (MQ=255)
                                                                                                   cGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCg  >  2:2125505/1‑100 (MQ=255)
                                                                                                   cGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCg  <  1:3413480/100‑1 (MQ=255)
                                                                                                   cGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCg  <  2:1567555/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCG  >  NC_000913/603296‑603494

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