Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 603,393 T→A 16.8% 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,3930TA16.8% 464.1 / 81.5 185intergenic (+58/+23)pheP/ybdGphenylalanine transporter/mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base T (78/75);  new base A (21/10);  total (99/86)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

Only 100 of 189 total aligned reads displayed.

GTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTT  >  NC_000913/603294‑603492
                                                                                                   |                                                                                                   
gTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAt                                                                                                     <  1:1102789/100‑1 (MQ=255)
 aaacctgcagtATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAtt                                                                                                    <  2:2702923/90‑1 (MQ=255)
  aTTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTc                                                                                                   >  2:2932038/1‑100 (MQ=255)
  aTTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTc                                                                                                   >  1:3148821/1‑100 (MQ=255)
   ttCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg                                                                                                  <  2:2250501/100‑1 (MQ=255)
    tCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGa                                                                                                 <  2:1811207/100‑1 (MQ=255)
     ccTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                                >  2:3912738/1‑100 (MQ=255)
     ccTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                                <  2:3866826/100‑1 (MQ=255)
        gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgag                                                                                             <  2:2416105/100‑1 (MQ=255)
        gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATga                                                                                              <  1:3994105/99‑1 (MQ=255)
           tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                          <  2:88513/100‑1 (MQ=39)
           tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                          <  2:4213098/100‑1 (MQ=39)
           tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                          <  2:1710709/100‑1 (MQ=39)
           tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                          <  1:583011/100‑1 (MQ=39)
           tATGGCATTTAAAACGCTGCGTCGGAAATAACGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                          <  1:1626578/100‑1 (MQ=38)
            aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCg                                                                                         <  2:3014568/100‑1 (MQ=255)
            aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCg                                                                                         <  2:1545782/100‑1 (MQ=38)
             tGGCATTTAAAGTGCTGCGTCGGAAATAAGGCATTCACGTTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                                                                                        <  1:3573783/100‑1 (MQ=255)
             tGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCgg                                                                                        <  1:677566/100‑1 (MQ=38)
              ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCgga                                                                                       >  2:677567/1‑99 (MQ=38)
              ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                                                                                        <  1:3616368/99‑1 (MQ=255)
              ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGt                                                                                       >  1:2753470/1‑100 (MQ=255)
                 aTTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttt                                                                                    >  1:4350625/1‑100 (MQ=255)
                  tttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGtttcg                                                                                   <  1:2577283/100‑3 (MQ=37)
                   ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttttt                                                                                  <  2:562748/100‑1 (MQ=255)
                   ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttttt                                                                                  >  2:1138453/1‑100 (MQ=255)
                       aaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAt                                                                              <  1:3351847/100‑1 (MQ=255)
                       aaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAt                                                                              <  2:309581/100‑1 (MQ=255)
                        aCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAtt                                                                             <  2:1844255/100‑1 (MQ=255)
                        aCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAtt                                                                             <  2:3147234/100‑1 (MQ=255)
                              cGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCt                                                                       >  2:3189249/1‑100 (MQ=255)
                              cGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCt                                                                       <  1:2833793/100‑1 (MQ=255)
                               gTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCtt                                                                      >  2:3483653/1‑100 (MQ=255)
                                  ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                   >  2:2010241/1‑100 (MQ=255)
                                   gAAATAAGGCATTCACGCTACATCCGACATAACGATGTCAACCATCCTAATCCGCTCTCATCCATTGGATGAGAGCGGTTTTTTTAATTACTGCTTAAAt                                                                  <  1:2650681/100‑1 (MQ=255)
                                    aaaTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                                 <  1:3815962/100‑1 (MQ=255)
                                     aaTAAGTCATTCACGCTACATCCGACAAAACTATGTCAACCATCCTAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                                 <  1:298254/99‑1 (MQ=255)
                                      aTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCa                                                               >  1:803499/1‑100 (MQ=255)
                                      aTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATACa                                                               >  2:370851/1‑100 (MQ=255)
                                        aaGGCATTCACGCTACATCCGACTAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAcc                                                             >  2:3737575/1‑100 (MQ=255)
                                        aaGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAcc                                                             >  2:2202887/1‑100 (MQ=255)
                                          ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCg                                                           <  2:1505144/100‑1 (MQ=255)
                                                cACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagag                                                     >  1:3981873/1‑100 (MQ=255)
                                                cACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagag                                                     >  2:3295663/1‑100 (MQ=255)
                                                  cGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagagag                                                   >  1:612965/1‑100 (MQ=255)
                                                  cGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagagag                                                   <  2:4164889/100‑1 (MQ=255)
                                                   gCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGc                                                  <  2:3061479/100‑1 (MQ=255)
                                                     tACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGa                                                >  1:2343449/1‑100 (MQ=255)
                                                      aCATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGaa                                               >  2:146055/1‑100 (MQ=255)
                                                        aTCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAata                                             >  2:263131/1‑100 (MQ=255)
                                                          ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATc                                           >  2:1327549/1‑100 (MQ=255)
                                                           cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCa                                          >  2:638302/1‑100 (MQ=255)
                                                           cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCa                                          >  1:3014438/1‑100 (MQ=255)
                                                              cAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTg                                       <  1:3904706/100‑1 (MQ=255)
                                                               aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGc                                      >  1:3317750/1‑100 (MQ=255)
                                                                aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                     <  1:540636/100‑1 (MQ=255)
                                                                aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                     >  2:3493458/1‑100 (MQ=255)
                                                                aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                     >  2:2130545/1‑100 (MQ=255)
                                                                aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                     <  1:4229482/100‑1 (MQ=255)
                                                                aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGc                                      >  2:3741281/1‑99 (MQ=255)
                                                                      tGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTgg                               >  1:773096/1‑100 (MQ=255)
                                                                       gTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGc                              >  1:3475822/1‑100 (MQ=255)
                                                                       gTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGc                              <  1:1864995/100‑1 (MQ=255)
                                                                        tCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCg                             >  1:132656/1‑100 (MQ=255)
                                                                         cAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGa                            >  1:4289927/1‑100 (MQ=255)
                                                                         cAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGa                            >  2:1310682/1‑100 (MQ=255)
                                                                          aaCCATCCGAAAGCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGAc                           >  2:2903603/1‑100 (MQ=255)
                                                                            ccATCCGAAACCGTTCTCATCCATTCGATGAGCGCTGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtg                         <  2:284380/100‑1 (MQ=255)
                                                                            ccATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTTTTTTTTTCTTAAATGCACACGCCAGATAGAGAATATTATTGCCGGttttagactg                         >  1:2967812/1‑92 (MQ=255)
                                                                            ccAGCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACt                          <  1:1712895/99‑1 (MQ=255)
                                                                            attaaaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGAcgg                         >  2:1727842/8‑98 (MQ=25)
                                                                             cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtga                        <  2:1882414/100‑1 (MQ=255)
                                                                             cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtga                        >  2:1219887/1‑100 (MQ=255)
                                                                               aaaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACACGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatg                      >  2:1248024/5‑100 (MQ=37)
                                                                                ccGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACttcatc                     <  1:1036431/100‑6 (MQ=255)
                                                                                 cGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGaa                    >  2:1270049/1‑100 (MQ=255)
                                                                                    aaCCGCTCTCCTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGATTATCATTGCCGGTTGGCGACTGATGAAGTc                 <  1:2903586/100‑1 (MQ=255)
                                                                                    aaCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTc                 >  2:4156551/1‑100 (MQ=255)
                                                                                     aCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTGAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAACTCg                >  1:2448815/1‑100 (MQ=37)
                                                                                     aCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCg                >  1:3223304/1‑100 (MQ=38)
                                                                                      ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc               >  1:1981677/1‑100 (MQ=38)
                                                                                      ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc               >  2:4242102/1‑100 (MQ=38)
                                                                                       cGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa              >  2:3467200/1‑100 (MQ=38)
                                                                                       cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa              <  2:2213950/100‑1 (MQ=255)
                                                                                        gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa              >  1:1216913/1‑99 (MQ=38)
                                                                                        gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa              >  1:3283141/1‑99 (MQ=38)
                                                                                        gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg             >  2:1531673/1‑100 (MQ=39)
                                                                                        gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg             >  2:2361998/1‑100 (MQ=39)
                                                                                        gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg             <  1:2686151/100‑1 (MQ=255)
                                                                                         ctctCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGa            >  1:3527686/1‑100 (MQ=39)
                                                                                          tctcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGAc           <  1:4170635/100‑1 (MQ=255)
                                                                                          tctcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGAc           <  2:3981903/100‑1 (MQ=255)
                                                                                            tcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCa         <  2:765114/100‑1 (MQ=255)
                                                                                             cATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaa        <  2:4350657/100‑1 (MQ=255)
                                                                                             cATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTAAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaa        >  1:2360337/1‑100 (MQ=255)
                                                                                              aTCCATTCGATGAGAGCGGTTTTTTTTTTTTTTTTTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaaa       >  1:3334526/1‑100 (MQ=255)
                                                                                               tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATTATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc      <  1:2350149/100‑1 (MQ=255)
                                                                                               tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGatcggaa      >  2:3572662/1‑94 (MQ=255)
                                                                                               tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc      <  2:1173614/100‑1 (MQ=255)
                                                                                               tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc      >  2:2859644/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTT  >  NC_000913/603294‑603492

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: