Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 603,391 C→G 18.9% intergenic (+56/+25) 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,3910CG18.9% 491.6 / 109.3 196intergenic (+56/+25)pheP/ybdGphenylalanine transporter/mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (84/75);  new base G (23/14);  total (107/89)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

Only 100 of 205 total aligned reads displayed.

TTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACT  >  NC_000913/603292‑603489
                                                                                                   |                                                                                                  
ttGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATcc                                                                                                    <  2:2689922/100‑1 (MQ=255)
ttGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATcc                                                                                                    <  1:1708312/100‑1 (MQ=255)
ttGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATcc                                                                                                    >  2:2472164/1‑100 (MQ=255)
  gTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCACCCGAAACCGCTCTCATCCAt                                                                                                  >  2:4381056/1‑100 (MQ=255)
  gTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACAAACAAAAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAt                                                                                                  <  1:2409393/100‑1 (MQ=255)
   tATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAtt                                                                                                 <  1:1282383/100‑1 (MQ=255)
   tATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAtt                                                                                                 >  2:440404/1‑100 (MQ=255)
    aTTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAAGCGCTCTCATCCATTc                                                                                                >  2:4341722/1‑100 (MQ=255)
    aTTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTc                                                                                                >  1:1659272/1‑100 (MQ=255)
     ttCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg                                                                                               >  1:3144480/1‑100 (MQ=255)
     ttCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg                                                                                               >  2:800198/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)
          gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgag                                                                                          >  2:1578514/1‑100 (MQ=255)
          gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgag                                                                                          >  1:1878456/1‑100 (MQ=255)
           tttATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgaga                                                                                         >  2:140514/1‑100 (MQ=255)
            ttATGGCATTTGAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGagag                                                                                        <  1:860046/100‑1 (MQ=255)
             tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                       <  2:4250978/100‑1 (MQ=39)
             tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                       <  2:1410067/100‑1 (MQ=39)
              aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                       <  1:2920606/99‑1 (MQ=38)
              aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCg                                                                                      <  2:270196/100‑1 (MQ=38)
              aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCg                                                                                      <  2:849917/100‑1 (MQ=255)
               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)
                ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGt                                                                                    >  1:1739005/1‑100 (MQ=255)
                ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGt                                                                                    >  2:4218805/1‑100 (MQ=255)
                 gCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGtt                                                                                   <  1:1681152/100‑1 (MQ=38)
                  cATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTAATCCATTCGATGAGAGCGGttt                                                                                  >  2:126948/1‑100 (MQ=255)
                   aTTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGtttc                                                                                 <  2:1895299/100‑2 (MQ=38)
                      tacaaCGCTGCCCCGGATATAAGGCATTCACGCTACATCCGACAATACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttttt                                                                              <  1:745301/97‑1 (MQ=255)
                        aaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTaa                                                                            <  1:2365234/100‑1 (MQ=255)
                          aCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAtt                                                                          <  2:4057160/100‑1 (MQ=255)
                          aCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTAATCCATTCGATGAGAGCGGTCTTTTTAAtt                                                                          >  1:3896145/1‑100 (MQ=255)
                             cTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACt                                                                       >  2:2714858/1‑100 (MQ=255)
                                cGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCt                                                                    >  2:119638/1‑100 (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:4360244/1‑100 (MQ=255)
                                            ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAACCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                              >  1:416913/1‑94 (MQ=255)
                                              cATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGcc                                                      <  2:1497938/100‑1 (MQ=255)
                                              cATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACACGcc                                                      >  2:22173/1‑100 (MQ=255)
                                                  cACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagag                                                  >  1:3682667/1‑100 (MQ=255)
                                                       tACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGa                                             <  2:2837627/100‑1 (MQ=255)
                                                        aCATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGTGCGa                                             >  1:4293532/1‑99 (MQ=255)
                                                        aCATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGaa                                            <  2:454111/100‑1 (MQ=255)
                                                        aCATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGaa                                            <  2:3049350/100‑1 (MQ=255)
                                                          aTCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAata                                          <  1:1593130/100‑1 (MQ=255)
                                                          aTCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAata                                          <  2:4107195/100‑1 (MQ=255)
                                                           tCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTCTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAata                                          >  1:340140/1‑99 (MQ=255)
                                                            ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTACTTACTGCTTAAATGCACCCGCCAGAGAGCGAATATc                                        >  1:2180833/1‑100 (MQ=255)
                                                            ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTATATGCACCCGCCAGAGAGCGAATATc                                        >  1:4312369/1‑100 (MQ=255)
                                                              gACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAt                                      <  1:1411939/100‑1 (MQ=255)
                                                                 aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGGTTTTGTAATTACTGCTTAAATGCACCCGCCAGAGGGCGAATATCa                                       >  1:4344995/1‑96 (MQ=255)
                                                                  aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                  >  2:216361/1‑100 (MQ=255)
                                                                   aaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCg                                 >  2:1651647/1‑100 (MQ=255)
                                                                     cGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGt                               >  2:3995214/1‑100 (MQ=255)
                                                                      gATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                              >  2:3402478/1‑100 (MQ=255)
                                                                       aTGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTg                             <  2:1682752/100‑1 (MQ=255)
                                                                          tCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCg                          >  2:3065968/1‑100 (MQ=255)
                                                                           ggaattaaaaaaaCCGCTCTCATCGAATGGACGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGa                         >  2:2964011/11‑100 (MQ=25)
                                                                           cAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGTTTGGCGa                         >  1:1253157/1‑100 (MQ=255)
                                                                            aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCg                          <  1:911995/98‑1 (MQ=255)
                                                                               cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtga                     >  2:2613360/1‑100 (MQ=255)
                                                                               cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtga                     >  1:2044275/1‑100 (MQ=255)
                                                                                aTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgat                    >  2:2417650/1‑100 (MQ=255)
                                                                                aTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgat                    <  2:1935348/100‑1 (MQ=255)
                                                                                aTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTGTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgat                    <  2:3094167/100‑1 (MQ=255)
                                                                                 tCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatg                   <  2:3801185/100‑1 (MQ=255)
                                                                                 aaaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatg                   >  2:2064068/5‑100 (MQ=37)
                                                                                 aaaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCAAGAGAGCGAATATCATTGCCGGTTGGCGACtgatg                   >  1:27767/5‑100 (MQ=37)
                                                                                  ccGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatga                  <  2:3436259/100‑1 (MQ=255)
                                                                                   cGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACttcatca                 <  2:888395/100‑7 (MQ=255)
                                                                                   cGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGaa                 >  1:774297/1‑100 (MQ=255)
                                                                                    gAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAg                <  2:3156849/100‑1 (MQ=255)
                                                                                     aaaCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGt               <  2:3714724/100‑1 (MQ=255)
                                                                                      aaCCGCTCTCATCCATTCGATGAGGGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGtn              <  2:1790085/100‑2 (MQ=255)
                                                                                       aCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCg             >  2:1617218/1‑100 (MQ=255)
                                                                                        ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTATTTTCTTCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc            >  1:410058/1‑100 (MQ=37)
                                                                                        ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc            >  1:1553397/1‑100 (MQ=38)
                                                                                        ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc            >  1:483265/1‑100 (MQ=38)
                                                                                        ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc            >  1:1899259/1‑100 (MQ=38)
                                                                                        ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc            >  2:2810236/1‑100 (MQ=38)
                                                                                        ccGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc            <  2:3298516/100‑1 (MQ=255)
                                                                                        ccGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc            <  2:1141861/100‑1 (MQ=255)
                                                                                         cGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa           >  2:3724521/1‑100 (MQ=38)
                                                                                         cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa           >  2:2381690/1‑100 (MQ=255)
                                                                                         cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa           <  1:1008979/100‑1 (MQ=255)
                                                                                         cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa           >  2:485805/1‑100 (MQ=255)
                                                                                          gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg          >  1:2988059/1‑100 (MQ=39)
                                                                                          gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg          >  2:3587571/1‑100 (MQ=39)
                                                                                          gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg          >  2:1301267/1‑100 (MQ=39)
                                                                                           ctctCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGa         >  1:3571783/1‑100 (MQ=39)
                                                                                           ctctCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGa         >  2:3866089/1‑100 (MQ=255)
                                                                                               cATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaa     <  1:2710187/100‑1 (MQ=255)
                                                                                                 tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc   <  2:1784061/100‑1 (MQ=255)
                                                                                                  ccATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAct  <  2:550899/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACT  >  NC_000913/603292‑603489

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: