Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 603,393 T→A 15.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,3930TA15.0% 407.6 / 53.1 158intergenic (+58/+23)pheP/ybdGphenylalanine transporter/mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base T (73/61);  new base A (16/8);  total (89/69)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

TATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTC  >  NC_000913/603295‑603490
                                                                                                  |                                                                                                 
 aTTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTc                                                                                                 <  1:1483674/100‑1 (MQ=255)
    ccTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                              >  2:1499887/1‑100 (MQ=255)
    ccTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                              <  2:1315204/100‑1 (MQ=255)
    ccTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                              >  1:2021982/1‑100 (MQ=255)
     cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                             <  1:945077/100‑1 (MQ=255)
     cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                             <  1:2941250/100‑1 (MQ=255)
     cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                             >  1:213910/1‑100 (MQ=255)
     cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                             >  1:2091287/1‑100 (MQ=255)
     cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCAGACAAAACGATATCAACCATCCGAAACCGCTATCATCCATTCg                                                                                                >  2:929761/1‑97 (MQ=255)
      tGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                              >  1:2168718/1‑98 (MQ=255)
      tGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATga                                                                                            >  1:3283521/1‑100 (MQ=255)
        tttATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgaga                                                                                          <  2:1367719/100‑1 (MQ=255)
         ttATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGagag                                                                                         <  2:1805503/100‑1 (MQ=255)
          tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAATCATCCGAAACCGCTCTCATCCATTCGATgagaag                                                                                        >  1:1729521/1‑98 (MQ=255)
          tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGagag                                                                                         <  1:1435640/99‑1 (MQ=39)
           aTGGCATTTAAACCGCTGCGTCGGAAATAAGGCATTCACGCTACACCCGACAAAACGATGTCAGCCATCCGAAACCGCTCTCATCGAATGGATGAGAGCg                                                                                       <  2:2430723/100‑1 (MQ=37)
           aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCg                                                                                       <  1:2786486/100‑1 (MQ=38)
           aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCg                                                                                       >  2:1653075/1‑100 (MQ=255)
           aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCg                                                                                       >  2:2506991/1‑100 (MQ=255)
            ttgCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCAAAACCGCTCTCATCGAATGGATGAGAGCgg                                                                                      <  2:2212878/98‑1 (MQ=37)
            tGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                                                                                      >  1:2640160/1‑100 (MQ=255)
             ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGt                                                                                     >  2:1601060/1‑100 (MQ=255)
              gCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCACTCATCCATTCGATGagag                                                                                         >  2:1129752/1‑95 (MQ=255)
                 tttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGtttcg                                                                                 <  1:1078283/100‑3 (MQ=37)
                 tttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttt                                                                                 >  1:2740425/1‑100 (MQ=255)
                  ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttttt                                                                                >  1:983250/1‑100 (MQ=255)
                  ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttttt                                                                                >  1:133728/1‑100 (MQ=255)
                   tAAAACGCTGCGTGGGAAATAAGGCATTCACACTACATCCGAGAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttttt                                                                               >  1:2460420/1‑100 (MQ=255)
                   tAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttttt                                                                               >  2:1896665/1‑100 (MQ=255)
                    aacgcGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGATACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTa                                                                              <  2:492898/96‑1 (MQ=255)
                    aaaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTa                                                                              <  1:2441957/100‑1 (MQ=255)
                     aaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTaa                                                                             <  1:786403/100‑1 (MQ=255)
                      aaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAt                                                                            <  1:1243602/100‑1 (MQ=255)
                      aaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAt                                                                            <  2:442322/100‑1 (MQ=255)
                       aCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAtt                                                                           >  1:11167/1‑100 (MQ=255)
                         gCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAGCCATCCGAATCCGCTCTCATCCATTCGATTAGAGCGGTTTTTTTa                                                                              <  2:1560217/95‑1 (MQ=255)
                         gCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTAc                                                                         >  1:2412131/1‑100 (MQ=255)
                          cTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACt                                                                        >  2:273402/1‑100 (MQ=255)
                            gCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGc                                                                      <  2:201623/100‑1 (MQ=255)
                             cGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCt                                                                     <  2:3283496/100‑1 (MQ=255)
                             cGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCt                                                                     >  2:2427064/1‑100 (MQ=255)
                                cGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaa                                                                  >  1:3578931/1‑100 (MQ=255)
                                 ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                 <  1:1896674/100‑1 (MQ=255)
                                 ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                 >  1:3070059/1‑100 (MQ=255)
                                 ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                 >  1:3047107/1‑100 (MQ=255)
                                  gAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAAt                                                                <  1:257625/100‑1 (MQ=255)
                                    aaTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGc                                                              >  1:2693516/1‑100 (MQ=255)
                                     aTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                               >  2:1063168/1‑98 (MQ=255)
                                      tAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAc                                                            >  2:3002174/1‑100 (MQ=255)
                                      tAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAc                                                            >  2:2318563/1‑100 (MQ=255)
                                        aGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAccc                                                          <  1:1767813/100‑1 (MQ=255)
                                         ggCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCg                                                         <  1:2633495/100‑1 (MQ=255)
                                           cATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGcc                                                       <  2:1879129/100‑1 (MQ=255)
                                              tCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCaga                                                    <  2:2728066/100‑1 (MQ=255)
                                               cACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagag                                                   >  1:3141196/1‑100 (MQ=255)
                                               cACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCaga                                                    >  1:2804991/1‑99 (MQ=255)
                                                      cATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAat                                            <  1:247916/100‑1 (MQ=255)
                                                        tCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAatat                                          >  1:848770/1‑100 (MQ=255)
                                                        tCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAatat                                          <  1:2277376/100‑1 (MQ=255)
                                                         ccGACAAAACGATGTCGAGCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATc                                         >  1:1590788/1‑100 (MQ=255)
                                                         ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATc                                         >  2:538411/1‑100 (MQ=255)
                                                          cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATc                                         >  1:3511895/1‑99 (MQ=255)
                                                          cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCa                                        >  2:30624/1‑100 (MQ=255)
                                                          cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCa                                        >  1:2388068/1‑100 (MQ=255)
                                                          cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCa                                        >  2:1059448/1‑100 (MQ=255)
                                                           gACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAt                                       <  2:2743126/100‑1 (MQ=255)
                                                           gACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAt                                       <  2:213910/100‑1 (MQ=255)
                                                             cAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTg                                     >  2:644610/1‑100 (MQ=255)
                                                             cAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTg                                     <  2:3029278/100‑1 (MQ=255)
                                                              aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGc                                    >  1:576857/1‑100 (MQ=255)
                                                                         aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGAc                         <  2:3340000/100‑1 (MQ=255)
                                                                         aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGAc                         <  1:1822396/100‑1 (MQ=255)
                                                                           ccATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCCTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtg                       <  2:2680618/100‑1 (MQ=255)
                                                                           attaaaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtg                       >  2:2142121/8‑100 (MQ=25)
                                                                             aTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgat                     <  1:3392066/100‑1 (MQ=255)
                                                                               ccGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatga                   <  1:333039/100‑1 (MQ=255)
                                                                               aaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatga                   >  1:2222244/4‑100 (MQ=37)
                                                                                cGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGGGCGAATATCATGGCCGGTTGGCGACTGATGaa                  >  2:1439690/1‑100 (MQ=255)
                                                                                 gAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAg                 <  1:3441958/100‑1 (MQ=255)
                                                                                   aaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTc               >  1:1876029/1‑100 (MQ=38)
                                                                                     ccGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGc             >  2:1270047/1‑100 (MQ=38)
                                                                                      cGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa            >  2:618318/1‑100 (MQ=38)
                                                                                      cGCTCTCAGGCATTCGATGAGAGCGGTGTTTTTAAGTACGGCTGAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCg              <  2:3382167/98‑1 (MQ=255)
                                                                                       gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGATAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa            >  2:2656434/1‑99 (MQ=38)
                                                                                       gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg           >  1:2939848/1‑100 (MQ=39)
                                                                                       gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg           >  2:3211228/1‑100 (MQ=39)
                                                                                       gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTCACGGTTGGCGACTGATGAAg                 >  2:3629191/1‑94 (MQ=37)
                                                                                       gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg           <  1:526979/100‑1 (MQ=255)
                                                                                       gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg           >  1:647847/1‑100 (MQ=255)
                                                                                       gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg           <  2:2816771/100‑1 (MQ=255)
                                                                                       gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg           >  2:2045338/1‑100 (MQ=255)
                                                                                       gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg           <  1:2843541/100‑1 (MQ=255)
                                                                                       gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg           <  1:987804/100‑1 (MQ=255)
                                                                                        ctctCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGa          >  1:2268883/1‑100 (MQ=39)
                                                                                         tctcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGAc         >  1:1613201/1‑100 (MQ=255)
                                                                                          ctcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGAcc        >  1:1174694/1‑100 (MQ=255)
                                                                                            cATCCATTCGATGAGAGCGGTTTTTTTTTTTTTTTTTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGTTAGGCGACTGATGAAGt                >  2:2309234/1‑90 (MQ=255)
                                                                                            cATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaa      >  2:464857/1‑100 (MQ=255)
                                                                                                cATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTTGCGACTGATGAAGTCGCAGACCAAActc  >  2:1407272/1‑100 (MQ=255)
                                                                                                cATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAActc  >  2:3230060/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
TATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTC  >  NC_000913/603295‑603490

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: