Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 603,395 C→G 17.3% 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,3950CG17.3% 491.4 / 89.4 185intergenic (+60/+21)pheP/ybdGphenylalanine transporter/mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (77/75);  new base G (21/11);  total (98/87)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 194 total aligned reads displayed.

ATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCG  >  NC_000913/603296‑603494
                                                                                                   |                                                                                                   
aTTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGACAACCATCCGAAACCGCTCTCATCCATTc                                                                                                     <  1:3835298/100‑1 (MQ=255)
 ttCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg                                                                                                    <  2:2250501/100‑1 (MQ=255)
   ccTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                                  >  2:3912642/1‑100 (MQ=255)
    cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                                 <  1:2754862/100‑1 (MQ=255)
     tGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATga                                                                                                >  2:2447503/1‑100 (MQ=255)
      gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgag                                                                                               >  2:3321042/1‑100 (MQ=255)
      gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgag                                                                                               >  1:2591616/1‑100 (MQ=255)
      gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATga                                                                                                <  1:3994105/99‑1 (MQ=255)
         tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                            <  2:88513/100‑1 (MQ=39)
         tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                            <  2:4213098/100‑1 (MQ=39)
         tATGGCATTTAAAACGCTGCGTCGGAAATAACGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGc                                                                                            <  1:1626578/100‑1 (MQ=38)
          aaggCATTTAAAAAGCTGCGTCGGAAATAAGCCATTCACGCTCCATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCg                                                                                           >  2:2975159/3‑100 (MQ=255)
          aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCg                                                                                           <  1:2471349/100‑1 (MQ=38)
          aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCg                                                                                           <  2:4138241/100‑1 (MQ=255)
           tGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                                                                                          <  1:1320218/100‑1 (MQ=255)
            ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCgga                                                                                         >  2:677567/1‑99 (MQ=38)
            ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGt                                                                                         >  1:3229748/1‑100 (MQ=255)
             gCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtt                                                                                        >  1:3811639/1‑100 (MQ=255)
             gCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtt                                                                                        >  2:302047/1‑100 (MQ=255)
                 ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttttt                                                                                    <  2:562748/100‑1 (MQ=255)
                    aaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTaa                                                                                 >  2:1812388/1‑100 (MQ=255)
                     aaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAt                                                                                <  1:3973776/100‑1 (MQ=255)
                     aaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAt                                                                                <  1:3351847/100‑1 (MQ=255)
                       cGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTa                                                                              <  1:2379992/100‑1 (MQ=255)
                          tGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTg                                                                           <  2:3195898/100‑1 (MQ=255)
                            cGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCt                                                                         >  2:3136668/1‑100 (MQ=255)
                              tCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTa                                                                       <  2:3079651/100‑1 (MQ=255)
                              tCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTa                                                                       <  2:1033586/100‑1 (MQ=255)
                                ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCtt                                                                        <  1:4291282/97‑1 (MQ=255)
                                   aaTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGc                                                                  <  2:3153089/100‑1 (MQ=255)
                                    aTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATACa                                                                 >  2:370851/1‑100 (MQ=255)
                                     tAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAc                                                                <  2:240955/100‑1 (MQ=255)
                                     tAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAc                                                                <  2:2139994/100‑1 (MQ=255)
                                      aaGGCATTCACGCTACATCCGACTAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAcc                                                               >  2:3737575/1‑100 (MQ=255)
                                        ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCg                                                             >  2:3939428/1‑100 (MQ=255)
                                        ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCg                                                             <  2:1505144/100‑1 (MQ=255)
                                         gCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGc                                                            >  1:2988779/1‑100 (MQ=255)
                                                cGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagagag                                                     >  1:284379/1‑100 (MQ=255)
                                                cGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagagag                                                     <  2:4164889/100‑1 (MQ=255)
                                                cGCTACATCCGACAAAACGATGTCAACCAGCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGTACCCGCCagagag                                                     <  1:1180412/100‑1 (MQ=255)
                                                 gCTACATCCGACAAAGCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGc                                                    >  2:1857319/1‑100 (MQ=255)
                                                   tACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGa                                                  >  1:2343449/1‑100 (MQ=255)
                                                    aCATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGaa                                                 <  2:2668113/100‑1 (MQ=255)
                                                    aCATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGaa                                                 >  2:146055/1‑100 (MQ=255)
                                                     cATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAat                                                <  2:1892211/100‑1 (MQ=255)
                                                      aTCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAata                                               >  2:178926/1‑100 (MQ=255)
                                                         cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCa                                            >  2:638302/1‑100 (MQ=255)
                                                          gACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAt                                           <  1:2675081/100‑1 (MQ=255)
                                                            cAACACGGTGTAAACCATCCGAAACTGCGATCATCCAGTAGATGAGATCGACATCTTTAATTACTGCTTAAATGCACCCGCCAGACAGCAAGTATCATTg                                         <  1:770450/100‑1 (MQ=255)
                                                            cAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTg                                         <  1:436836/100‑1 (MQ=255)
                                                             aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGc                                        >  1:3317750/1‑100 (MQ=255)
                                                              aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                       <  1:4229482/100‑1 (MQ=255)
                                                              aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                       >  2:3493458/1‑100 (MQ=255)
                                                              aaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGc                                        >  2:3741281/1‑99 (MQ=255)
                                                                ccGATCTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCgg                                     <  1:656297/99‑1 (MQ=255)
                                                                     gTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGc                                >  1:3475822/1‑100 (MQ=255)
                                                                     gTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGc                                >  2:1724146/1‑100 (MQ=255)
                                                                     gTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGc                                <  1:1864995/100‑1 (MQ=255)
                                                                       cAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGa                              >  2:1310682/1‑100 (MQ=255)
                                                                        aaCCATCCGAAAGCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGAc                             >  2:2903603/1‑100 (MQ=255)
                                                                        aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGAc                             <  2:3469128/100‑1 (MQ=255)
                                                                         aCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACt                            >  1:2688819/1‑100 (MQ=255)
                                                                         aCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACt                            >  1:567496/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)
                                                                           cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtga                          >  2:1219887/1‑100 (MQ=255)
                                                                           cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtga                          >  2:203056/1‑100 (MQ=255)
                                                                             aaaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACACGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatg                        >  2:1248024/5‑100 (MQ=37)
                                                                              ccGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACttcatc                       <  1:1036431/100‑6 (MQ=255)
                                                                              ccGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatga                       <  1:3073468/100‑1 (MQ=255)
                                                                               cGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGaa                      >  2:1270049/1‑100 (MQ=255)
                                                                                  aaCCGCTCTCCTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGATTATCATTGCCGGTTGGCGACTGATGAAGTc                   <  1:2903586/100‑1 (MQ=255)
                                                                                  aaCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTCGCGACTGATGAAGTc                   >  1:586041/1‑100 (MQ=255)
                                                                                   aCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTACTTACTACTTAAATGCACGAGCCAGAGAGCGAATATCATTACCAGTTGGCGCCtaaaaacgt                    >  1:3222335/1‑90 (MQ=16)
                                                                                   aCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTATTTTTTTCTTAAATGCACCAGCCAGAGAGCGAATATCATTTCCGATTGACGACtgttta                       >  1:4030115/1‑91 (MQ=255)
                                                                                    ttccgatcttCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 <  1:3572637/91‑1 (MQ=255)
                                                                                    ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 >  2:4242102/1‑100 (MQ=38)
                                                                                    ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 >  1:1981677/1‑100 (MQ=38)
                                                                                    ccGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc                 >  1:801487/1‑100 (MQ=255)
                                                                                     cGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                <  2:2213950/100‑1 (MQ=255)
                                                                                      gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               >  2:1531673/1‑100 (MQ=39)
                                                                                      gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               >  2:2361998/1‑100 (MQ=39)
                                                                                        tctcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGAc             <  1:4170635/100‑1 (MQ=255)
                                                                                        tctcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGAc             <  1:3581167/100‑1 (MQ=255)
                                                                                        tctcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGAc             <  1:1428686/100‑1 (MQ=255)
                                                                                          tcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCa           <  2:765114/100‑1 (MQ=255)
                                                                                           cATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaa          <  2:4350657/100‑1 (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:2859644/1‑100 (MQ=255)
                                                                                             tCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAAc        <  2:1173614/100‑1 (MQ=255)
                                                                                               cATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAActc      >  2:168071/1‑100 (MQ=255)
                                                                                                aTTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAActct     >  2:1311799/1‑100 (MQ=255)
                                                                                                 ttCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCtt    >  2:3847641/1‑100 (MQ=255)
                                                                                                  tCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTc   <  1:3297698/100‑1 (MQ=255)
                                                                                                   cGATGAGAGCGGTTTTTTTTTTTTTTTTTTAAATAGACCAGCAACAGCGCGATTTTTATTGCCGGTTGGCGAATGATGAAGTCGCTAACCaacatattcg  >  1:966507/1‑92 (MQ=255)
                                                                                                   cGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCg  <  2:1340408/100‑1 (MQ=255)
                                                                                                   cGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCg  <  2:3668483/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: