Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 603,395 C→G 15.2% 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,3950CG15.2% 434.8 / 65.8 158intergenic (+60/+21)pheP/ybdGphenylalanine transporter/mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (74/60);  new base G (16/8);  total (90/68)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.36e-01

Only 100 of 163 total aligned reads displayed.

ATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCG  >  NC_000913/603296‑603494
                                                                                                   |                                                                                                   
 ttCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCg                                                                                                    <  1:3137749/100‑1 (MQ=255)
   ccTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                                  <  2:1315204/100‑1 (MQ=255)
   ccTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                                  >  2:1499887/1‑100 (MQ=255)
   ccTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGAt                                                                                                  <  2:3557069/100‑1 (MQ=255)
    cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                                 <  2:3274568/100‑1 (MQ=255)
    cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                                 >  2:3012261/1‑100 (MQ=255)
    cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                                 >  1:213910/1‑100 (MQ=255)
    cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATg                                                                                                 >  1:2091287/1‑100 (MQ=255)
    cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCAGACAAAACGATATCAACCATCCGAAACCGCTATCATCCATTCg                                                                                                    >  2:929761/1‑97 (MQ=255)
     tGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATga                                                                                                >  1:3283521/1‑100 (MQ=255)
        ttATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGagag                                                                                             <  2:1805503/100‑1 (MQ=255)
        ttATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACAATTCGAAACCGCTCTCATCCATTCGATGagag                                                                                             >  2:25947/1‑100 (MQ=255)
         tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGagag                                                                                             <  1:1435640/99‑1 (MQ=39)
          ggtgCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCg                                                                                           <  1:2828105/97‑1 (MQ=255)
          aTGGCATTTAAACCGCTGCGTCGGAAATAAGGCATTCACGCTACACCCGACAAAACGATGTCAGCCATCCGAAACCGCTCTCATCGAATGGATGAGAGCg                                                                                           <  2:2430723/100‑1 (MQ=37)
          aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCg                                                                                           <  1:2786486/100‑1 (MQ=38)
          aTGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCg                                                                                           >  2:1653075/1‑100 (MQ=255)
          aTGGCATTTAAAACGCTGCGTCGAAAATAAGGCATTCACGCTACATCCTACAAAACGATGTCAACCATCCGAAACCGCTCTCATCAAATGGATGAGAGCg                                                                                           <  2:1169556/100‑1 (MQ=38)
           tcgCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCAAAACCGCTCTCATCGAATGGATGAGAGCgg                                                                                          <  2:2213091/98‑1 (MQ=38)
           tGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                                                                                          >  1:2640160/1‑100 (MQ=255)
            ggCATTTAAGACGCTGCGTCGGAAAGAAGGCATTCGCGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGt                                                                                         <  1:1586948/100‑1 (MQ=37)
            ggCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCg                                                                                           >  2:2167266/1‑98 (MQ=255)
             gCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCACTCATCCATTCGATGagag                                                                                             >  2:1129752/1‑95 (MQ=255)
                tttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCGGtttcg                                                                                     <  1:1078283/100‑3 (MQ=37)
                tttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGATCggaagag                                                                                     >  2:2725160/1‑95 (MQ=255)
                tttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttt                                                                                     >  1:2740425/1‑100 (MQ=255)
                 ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttttt                                                                                    >  1:133728/1‑100 (MQ=255)
                  tAAAACGCTGCGTGGGAAATAAGGCATTCACACTACATCCGAGAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttttt                                                                                   >  1:2460420/1‑100 (MQ=255)
                  tAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttttttt                                                                                   >  1:934809/1‑100 (MQ=255)
                  tAAAACGCTGCGTCGCAAATAAGGCATTCACGCTACATCCCACAAAACGATGTCAACCATCCGTAACCGCTCTCATCCATTCGATGAGAGCGGtttttt                                                                                    <  2:3335150/99‑1 (MQ=255)
                   aaaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTa                                                                                  <  1:2441957/100‑1 (MQ=255)
                   aaaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTa                                                                                  <  1:1929175/100‑1 (MQ=255)
                     aaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAt                                                                                <  2:442322/100‑1 (MQ=255)
                     aaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAt                                                                                <  1:1243602/100‑1 (MQ=255)
                      aCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAtt                                                                               >  1:11167/1‑100 (MQ=255)
                       cGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTa                                                                              <  2:1564962/100‑1 (MQ=255)
                        gCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTAc                                                                             >  1:2412131/1‑100 (MQ=255)
                         cTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACt                                                                            >  2:273402/1‑100 (MQ=255)
                           gCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGc                                                                          <  2:201623/100‑1 (MQ=255)
                             gTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCt                                                                         >  2:2381999/1‑99 (MQ=255)
                               cGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaa                                                                      >  2:2106695/1‑100 (MQ=255)
                               cGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaa                                                                      >  1:3578931/1‑100 (MQ=255)
                                ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                     >  1:3047107/1‑100 (MQ=255)
                                ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTaaa                                                                     <  1:1896674/100‑1 (MQ=255)
                                 gAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAAt                                                                    <  1:257625/100‑1 (MQ=255)
                                  aaaTAAGGCAATCACGCTACATCCGACAAAACTATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                                   <  2:2444848/100‑1 (MQ=255)
                                   aaTAAGGCATTCACGCTACATCCGACAAAGCGATGTCAACCATCCTAAACCGCTTTCATCCATTCGATGAGAGCTGTTTTTTTAATTACTGCTTAAATGc                                                                  <  2:476268/100‑1 (MQ=255)
                                   aaTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGc                                                                  >  1:2693516/1‑100 (MQ=255)
                                    aTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                                   >  2:1063168/1‑98 (MQ=255)
                                     tAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAc                                                                >  1:3480843/1‑100 (MQ=255)
                                     tAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAc                                                                >  2:3002174/1‑100 (MQ=255)
                                       aGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCAccc                                                              <  1:1767813/100‑1 (MQ=255)
                                              cACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagag                                                       >  1:3141196/1‑100 (MQ=255)
                                                     cATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAat                                                >  1:3340025/1‑100 (MQ=255)
                                                     cATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAat                                                <  1:247916/100‑1 (MQ=255)
                                                       tCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAatat                                              <  1:2277376/100‑1 (MQ=255)
                                                       tCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAatat                                              >  1:1818994/1‑100 (MQ=255)
                                                       tCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAatat                                              >  1:848770/1‑100 (MQ=255)
                                                        ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATc                                             >  2:538411/1‑100 (MQ=255)
                                                         cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCa                                            >  1:2388068/1‑100 (MQ=255)
                                                         cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCa                                            >  2:30624/1‑100 (MQ=255)
                                                          gACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACTCGCCAGAGAGCGAATATCAt                                           >  2:3098379/1‑100 (MQ=255)
                                                          gACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAt                                           <  2:2743126/100‑1 (MQ=255)
                                                           aCAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAtt                                          <  1:96712/100‑1 (MQ=255)
                                                            cAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTg                                         <  1:748167/100‑1 (MQ=255)
                                                            cAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTg                                         >  2:644610/1‑100 (MQ=255)
                                                             aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGc                                        >  1:576857/1‑100 (MQ=255)
                                                                     gTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGc                                <  1:2155532/100‑1 (MQ=255)
                                                                        aaCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGa                              >  2:1886294/1‑99 (MQ=255)
                                                                         aCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACt                            >  1:2266414/1‑100 (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:996338/100‑1 (MQ=255)
                                                                              ccGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatg                        <  1:549057/99‑1 (MQ=255)
                                                                              aaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatga                       >  2:1356944/4‑100 (MQ=37)
                                                                              aaaaaaCCGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatga                       >  1:2222244/4‑100 (MQ=37)
                                                                               cGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGGGCGAATATCATGGCCGGTTGGCGACTGATGaa                      >  2:1439690/1‑100 (MQ=255)
                                                                                gAATCCGCTCTCATCCATTCGATGAGAGCAGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAg                     <  1:428489/100‑1 (MQ=255)
                                                                                gAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAg                     >  1:30338/1‑100 (MQ=255)
                                                                                     cGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                >  2:618318/1‑100 (MQ=38)
                                                                                     cGCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCAGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                >  1:2635596/1‑100 (MQ=38)
                                                                                     cGCTCTCAGGCATTCGATGAGAGCGGTGTTTTTAAGTACGGCTGAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCg                  <  2:3382167/98‑1 (MQ=255)
                                                                                      gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               >  2:3211228/1‑100 (MQ=39)
                                                                                      gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               >  1:2939848/1‑100 (MQ=39)
                                                                                      gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               >  1:1664429/1‑100 (MQ=255)
                                                                                      gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               <  1:987804/100‑1 (MQ=255)
                                                                                      gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               <  2:2816771/100‑1 (MQ=255)
                                                                                      gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               <  1:2843541/100‑1 (MQ=255)
                                                                                       ctctCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAg               >  1:186/1‑99 (MQ=39)
                                                                                       ctctCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGa              >  1:1450978/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)
                                                                                            aTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaaa         >  2:1978440/1‑100 (MQ=255)
                                                                                               cATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTTGCGACTGATGAAGTCGCAGACCAAActc      >  2:1407272/1‑100 (MQ=255)
                                                                                                  tCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTc   >  2:1617576/1‑100 (MQ=255)
                                                                                                   cGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCg  >  2:3285969/1‑100 (MQ=255)
                                                                                                   cGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCg  <  2:1878997/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: