Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 603,393 T→A 19.0% intergenic (+58/+23) 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,3930TA19.0% 444.2 / 98.1 197intergenic (+58/+23)pheP/ybdGphenylalanine transporter/mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base T (86/73);  new base A (23/14);  total (110/87)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-02

Only 100 of 204 total aligned reads displayed.

GTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCT  >  NC_000913/603294‑603491
                                                                                                   |                                                                                                  
gTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAt                                                                                                    >  2:4207159/1‑100 (MQ=255)
gTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACAAACAAAAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAt                                                                                                    <  1:2409393/100‑1 (MQ=255)
 tATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAtt                                                                                                   <  1:1282383/100‑1 (MQ=255)
   ttCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg                                                                                                 >  1:3144480/1‑100 (MQ=255)
     ccTGTTTATGGCATTTAAAACGCTGCGTCGGACATAAGGCATTCACGCTACATCCGACAAAACGCTCTCAACCATCCGAAACCGCTCTCCTCCCTTCg                                                                                                 >  1:3814221/1‑98 (MQ=255)
      cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                              <  1:1122321/100‑1 (MQ=255)
      cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                              <  1:3740200/100‑1 (MQ=255)
       tGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATga                                                                                             <  2:1283439/100‑1 (MQ=255)
        gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgag                                                                                            >  1:1878456/1‑100 (MQ=255)
        gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgag                                                                                            >  2:1578514/1‑100 (MQ=255)
        gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgag                                                                                            >  2:179353/1‑100 (MQ=255)
           tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGc                                                                                         <  2:1109509/100‑1 (MQ=255)
            aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCg                                                                                        <  1:1625508/100‑1 (MQ=38)
            aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCg                                                                                        >  2:2768893/1‑100 (MQ=255)
             tGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                                                                                       >  1:1672467/1‑100 (MQ=255)
             tGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                                                                                       >  1:3914416/1‑100 (MQ=255)
              ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAATGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGt                                                                                      <  1:2441768/100‑1 (MQ=38)
              ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGt                                                                                      <  2:2692238/100‑1 (MQ=38)
              ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGt                                                                                      >  1:1739005/1‑100 (MQ=255)
                cATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttt                                                                                    >  2:1405457/1‑100 (MQ=255)
                    tAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttttt                                                                                >  2:2393226/1‑100 (MQ=255)
                     aaaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTa                                                                               <  1:585772/100‑1 (MQ=255)
                     aaaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTa                                                                               <  2:559612/100‑1 (MQ=255)
                      aaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTaa                                                                              <  1:2365234/100‑1 (MQ=255)
                       aaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAt                                                                             >  1:1682747/1‑100 (MQ=255)
                        aCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTAATCCATTCGATGAGAGCGGTCTTTTTAAtt                                                                            >  1:3896145/1‑100 (MQ=255)
                          gCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTAc                                                                          >  2:1859096/1‑100 (MQ=255)
                             gCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTg                                                                        <  2:4355410/99‑1 (MQ=255)
                               gTCGGAAATAAGGCATTCACGCTACATCCGACAAAATGATGTCACCCATCCGACACCGTTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCtt                                                                     >  1:2100566/1‑100 (MQ=255)
                               gTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCtt                                                                     <  1:2260126/100‑1 (MQ=255)
                                  ggACATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCCCTCATCCACTCGATGAGAGCGGCTTTTTTAATTACTGCTTcaa                                                                  >  1:3650033/1‑100 (MQ=255)
                                  ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                  <  2:1898215/100‑1 (MQ=255)
                                  ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                  >  1:2655741/1‑100 (MQ=255)
                                       tAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAc                                                             >  2:886861/1‑100 (MQ=255)
                                          ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAACCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                                >  1:416913/1‑94 (MQ=255)
                                          ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCg                                                          >  2:3913540/1‑100 (MQ=255)
                                          ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCg                                                          >  1:800088/1‑100 (MQ=255)
                                            cATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGcc                                                        <  2:1497938/100‑1 (MQ=255)
                                            cATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACACGcc                                                        >  2:22173/1‑100 (MQ=255)
                                                cACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagag                                                    >  1:888391/1‑100 (MQ=255)
                                                cACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagag                                                    >  1:3682667/1‑100 (MQ=255)
                                                  cGCTACAGCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagaga                                                   <  2:4293762/99‑1 (MQ=255)
                                                    cTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCg                                                <  2:4360776/100‑1 (MQ=255)
                                                      aCATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGTGCGa                                               >  1:4293532/1‑99 (MQ=255)
                                                         tCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTCTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAata                                            >  1:340140/1‑99 (MQ=255)
                                                          ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTACTTACTGCTTAAATGCACCCGCCAGAGAGCGAATATc                                          >  1:2180833/1‑100 (MQ=255)
                                                          ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTATATGCACCCGCCAGAGAGCGAATATc                                          >  1:4312369/1‑100 (MQ=255)
                                                           cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCa                                         <  2:4230710/100‑1 (MQ=255)
                                                            gACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAt                                        <  1:1411939/100‑1 (MQ=255)
                                                             aCAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAtt                                       >  2:2710202/1‑100 (MQ=255)
                                                              cAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTg                                      <  2:4381680/100‑1 (MQ=255)
                                                               aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGc                                     >  2:1187041/1‑100 (MQ=255)
                                                               aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGGTTTTGTAATTACTGCTTAAATGCACCCGCCAGAGGGCGAATATCa                                         >  1:4344995/1‑96 (MQ=255)
                                                                aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                    <  2:2465862/100‑1 (MQ=255)
                                                                aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                    >  2:216361/1‑100 (MQ=255)
                                                                 aaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCg                                   >  2:1651647/1‑100 (MQ=255)
                                                                  aCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCgg                                  <  2:3843259/100‑1 (MQ=255)
                                                                    gATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                                >  2:3402478/1‑100 (MQ=255)
                                                                    gATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                                >  2:282654/1‑100 (MQ=255)
                                                                     aTGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTg                               <  2:1682752/100‑1 (MQ=255)
                                                                      tGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                                >  1:757721/1‑98 (MQ=255)
                                                                      tGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTgg                              <  2:531715/100‑1 (MQ=255)
                                                                       gTCAACCATCCGAAACCGCTCTCATCCATGCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGc                             <  2:304986/100‑1 (MQ=255)
                                                                        tCAACCATCCGAAACCGCTCTCATGCATTCGATGAGAGCGGTTTCTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCg                            <  2:839589/100‑1 (MQ=255)
                                                                        tCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCg                            <  2:1978290/100‑1 (MQ=255)
                                                                        tCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCg                            >  2:484385/1‑100 (MQ=255)
                                                                         cAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGTTTGGCGa                           >  1:1253157/1‑100 (MQ=255)
                                                                             cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtga                       >  2:2613360/1‑100 (MQ=255)
                                                                             cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtga                       >  1:2044275/1‑100 (MQ=255)
                                                                              aTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgat                      <  2:126429/100‑1 (MQ=255)
                                                                              aTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTGTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgat                      <  2:3094167/100‑1 (MQ=255)
                                                                               aaaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatg                     >  2:2064068/5‑100 (MQ=37)
                                                                                ccGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatg                     <  1:292320/99‑1 (MQ=255)
                                                                                  gAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAg                  <  2:3156849/100‑1 (MQ=255)
                                                                                   ataCCGCTCTCATAGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGt                 >  2:2989613/3‑100 (MQ=38)
                                                                                   aaaCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGt                 <  2:3714724/100‑1 (MQ=255)
                                                                                     aCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCg               >  2:1617218/1‑100 (MQ=255)
                                                                                      ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTATTTTCTTCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc              >  1:410058/1‑100 (MQ=37)
                                                                                      ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAGCAGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc              >  1:2880758/1‑100 (MQ=38)
                                                                                      ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTGGc              >  1:3005526/1‑100 (MQ=38)
                                                                                      ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc              >  1:1899259/1‑100 (MQ=38)
                                                                                      ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc              >  1:483265/1‑100 (MQ=38)
                                                                                      ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATACCATTGCCGGTTGGCGACTGATGAAGTCGc              >  2:1305126/1‑100 (MQ=38)
                                                                                      ccGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc              <  2:3470602/100‑1 (MQ=255)
                                                                                      ccGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc              >  2:4002719/1‑100 (MQ=255)
                                                                                      ccGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc              <  2:3298516/100‑1 (MQ=255)
                                                                                      ccGCTCTCATCCATTCGATGAGAGCGGTTGTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc              <  2:3793100/100‑1 (MQ=255)
                                                                                       cGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa             >  2:4075530/1‑100 (MQ=38)
                                                                                       cGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa             >  1:1583264/1‑100 (MQ=39)
                                                                                       cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa             <  1:1008979/100‑1 (MQ=255)
                                                                                        gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTTTTTTTTAAATGCACCCGCCAGAGAGCGACTATCATTGCCGGTTGGCGACTGATGAAGTCTCAg            >  1:2742830/1‑100 (MQ=25)
                                                                                        gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg            >  2:3587571/1‑100 (MQ=39)
                                                                                         ctctCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGa           >  1:3571783/1‑100 (MQ=39)
                                                                                             cATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaa       >  2:2915632/1‑100 (MQ=255)
                                                                                               tCCATTCGATGAGAGCGGTTTTTTTATTTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc     <  2:441931/100‑1 (MQ=255)
                                                                                               tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc     <  2:1787887/100‑1 (MQ=255)
                                                                                                  aTTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAActct  >  2:2290954/1‑100 (MQ=255)
                                                                                                  aTTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAActct  >  2:1729503/1‑100 (MQ=255)
                                                                                                  aTTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAActct  >  2:928781/1‑100 (MQ=255)
                                                                                                   ttCGATGAGAGCGGTTTTTTTAATTACTTCTTAAATGCACCCGCCAGCGTGCGAATATCATTGCCGTTTGGCTACTGATGAAGTCTCAAACCAAActct  >  1:1494924/1‑99 (MQ=255)
                                                                                                   |                                                                                                  
GTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCT  >  NC_000913/603294‑603491

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: