Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,087,771 C→T 5.7% R163Q (CGG→CAG)  yeeY ← LysR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,087,7710CT 5.7% 223.7 / 2.4 70R163Q (CGG→CAG) yeeYLysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (34/31);  new base T (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

TGATGCCGTAGAAGCGCCTGCCTGCGCAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACGCTATGATGTTGATAAAAAGC  >  NC_000913/2087626‑2087919
                                                                                                                                                 |                                                                                                                                                    
tgATGCCGTAGAAGCGCCTGCCTGCGCAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGt                                                                                                                                                    <  1:64352/148‑1 (MQ=255)
   tGCCGTAGAAGCGCCTGCCTGCGCAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCAAGATCGCCCAACCGTGaa                                                                                                                                                 >  6:80074/1‑148 (MQ=255)
      cGTAGAAGCGCCTGCCTGCGCAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAAc                                                                                                                                             <  3:254731/149‑1 (MQ=255)
        tAGAAGCGCCTGCCTGCGCAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAAcc                                                                                                                                            >  3:327320/1‑148 (MQ=255)
         aGAAGCGCCTGCCTGCGCAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCa                                                                                                                                           >  4:320919/1‑148 (MQ=255)
          gAAGCGCCTGCCTGCGCAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCaaa                                                                                                                                         <  2:6621/149‑1 (MQ=255)
          gAAGCGCCTGCCTGCGCAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCaaa                                                                                                                                         >  3:341114/1‑149 (MQ=255)
           aaGCGCCTGCTTGCGCAGTATCCACCACCACAATCGCCCTATATCGCTTAATAATCCGACGATTTAACGGCTGTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCTACCTTGAAAAACCATAc                                                                                                                                        <  6:60202/149‑1 (MQ=255)
               gcctgcctGCGCAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCAc                                                                                                                                    <  4:341114/149‑1 (MQ=255)
                cctgcctgCGCAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACt                                                                                                                                   <  7:168864/149‑1 (MQ=255)
                      tgcGCAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTgg                                                                                                                                 >  8:392669/1‑145 (MQ=255)
                        cgcAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGAtg                                                                                                                           <  7:302919/149‑1 (MQ=255)
                         gcAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGAtgg                                                                                                                          <  7:392669/149‑1 (MQ=255)
                           aGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGAtggtg                                                                                                                        <  8:33986/149‑1 (MQ=255)
                                   ccaccaCAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGc                                                                                                                <  3:146989/149‑1 (MQ=255)
                                    caccacAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCAc                                                                                                                                    >  5:334342/1‑128 (MQ=255)
                                    caccacAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCAc                                                                                                                                    <  6:334342/128‑1 (MQ=255)
                                      ccacAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCAt                                                                                                              >  3:52433/1‑148 (MQ=255)
                                      ccacAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCAt                                                                                                              <  4:52433/148‑1 (MQ=255)
                                       cacaATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATcg                                                                                                            >  4:167869/1‑149 (MQ=255)
                                       cacaATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATcg                                                                                                            >  2:198618/1‑149 (MQ=255)
                                        acaATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATcgc                                                                                                           <  3:167869/149‑1 (MQ=255)
                                                      tCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAataata                                                                                              <  3:211739/148‑1 (MQ=255)
                                                                 taCGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGcccc                                                                                  >  1:127702/1‑149 (MQ=255)
                                                                 taCGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGcccc                                                                                  >  7:351044/1‑149 (MQ=255)
                                                                 taCGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACAGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGcccc                                                                                  >  5:107414/1‑149 (MQ=255)
                                                                    gacgaTTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACTGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCt                                                                                 >  1:239488/1‑147 (MQ=255)
                                                                    gacgaTTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACTGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCt                                                                                 >  1:245612/1‑147 (MQ=255)
                                                                    gacgaTTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACaata                                                                                                 >  6:341576/1‑131 (MQ=255)
                                                                    gacgaTTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACaata                                                                                                 <  5:341576/131‑1 (MQ=255)
                                                                    gacgaTTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCt                                                                                 >  1:182466/1‑147 (MQ=255)
                                                                    gacgaTTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGa                                                                               <  4:327320/149‑1 (MQ=255)
                                                                      cgaTTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGCCCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGt                                                                             >  7:336645/1‑149 (MQ=255)
                                                                         tttAACGGCTCTTCTTCCCGGGCTAAAGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCATACTCACTGGACGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAg                                                                          <  5:223517/149‑1 (MQ=255)
                                                                             aCGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCa                                                                      >  6:193074/1‑149 (MQ=255)
                                                                                gCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATcc                                                                   <  4:142005/149‑1 (MQ=255)
                                                                                  tcttctTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCa                                                                 >  6:351255/1‑149 (MQ=255)
                                                                                  tcttctTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCa                                                                 >  4:355287/1‑149 (MQ=255)
                                                                                  tcttctTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCa                                                                 >  1:363171/1‑149 (MQ=255)
                                                                                    ttcttcCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCAt                                                                >  7:140756/1‑148 (MQ=255)
                                                                                       ttcCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCa                                                                 <  7:199366/144‑1 (MQ=255)
                                                                                       ttcCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCa                                                                 >  8:199366/1‑144 (MQ=255)
                                                                                           aGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACTGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAaccacc                                                        <  2:239488/149‑1 (MQ=255)
                                                                                           aGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACTGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAaccacc                                                        <  2:245612/149‑1 (MQ=255)
                                                                                           aGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAaccacc                                                        >  1:111162/1‑149 (MQ=255)
                                                                                           aGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAaccacc                                                        <  2:127702/149‑1 (MQ=255)
                                                                                           aGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAaccacc                                                        <  2:182466/149‑1 (MQ=255)
                                                                                                tAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTa                                                   >  2:68324/1‑149 (MQ=255)
                                                                                                       tgatgaGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCAtt                                            <  8:140756/149‑1 (MQ=255)
                                                                                                              gcgcGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGGCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAAc                                     <  1:198618/149‑1 (MQ=255)
                                                                                                                gcgACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTAc                                                  <  1:323263/134‑1 (MQ=255)
                                                                                                                  gACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCa                                 <  1:68324/149‑1 (MQ=255)
                                                                                                                    ccGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCcaacaa                                                                                                   >  6:100513/1‑81 (MQ=255)
                                                                                                                      gCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGgcg                             >  6:6522/1‑149 (MQ=255)
                                                                                                                        aaaGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGcgcg                           >  3:59101/1‑149 (MQ=255)
                                                                                                                           gACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCa                        >  8:165649/1‑149 (MQ=255)
                                                                                                                             cTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACg                      >  4:254362/1‑149 (MQ=255)
                                                                                                                             cTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACg                      >  2:307739/1‑149 (MQ=255)
                                                                                                                              ttGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGa                                          >  4:366357/1‑128 (MQ=255)
                                                                                                                              ttGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGa                                          <  3:366357/128‑1 (MQ=255)
                                                                                                                                gTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACGCTa                   >  8:294289/1‑149 (MQ=255)
                                                                                                                                 ttgcAGATCGCCCAACCGTTAAAAACCAAACTCAGTGGAAGATGGTGGGTCATGCATCGCGCCAACAATAATATCCGCTCGCGCCTGAGTCAGGGCATCCGATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACGCTat                  <  3:339963/146‑1 (MQ=255)
                                                                                                                                    cAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACGCTatga                <  2:111162/148‑1 (MQ=255)
                                                                                                                                     aGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGCGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACGCTatgatg              >  5:360760/1‑149 (MQ=255)
                                                                                                                                          gCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACGCTATGATGTTGa          <  8:351044/148‑1 (MQ=255)
                                                                                                                                           cccAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCCTCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCAGTGATAAACATCAGGCGCGTCACGCTATGATGTTGAt         >  2:385152/1‑148 (MQ=255)
                                                                                                                                           cccAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACGCTATGATGTTGAt         >  6:213374/1‑148 (MQ=255)
                                                                                                                                              aaCCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACGCTATGATGTTGATaaaa     <  3:320919/149‑1 (MQ=255)
                                                                                                                                               aCCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAg                                >  4:185530/1‑121 (MQ=255)
                                                                                                                                               aCCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAg                                <  3:185530/121‑1 (MQ=255)
                                                                                                                                                ccGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACGCTATGATGTTGATAAAAAg   <  6:107414/149‑1 (MQ=255)
                                                                                                                                                 cGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACGCTATGATGTTGATAAAAAGc  >  1:288444/1‑149 (MQ=255)
                                                                                                                                                 cGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACGCTATGATGTTGATAAAAAGc  >  8:412529/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TGATGCCGTAGAAGCGCCTGCCTGCGCAGTATCCCCCACCACAATCGCCCGATATCGCTTAATAATACGACGATTTAACGGCTCTTCTTCCAGGGCTAACGGATGATGAGGCGCGACCGCAAAGACTTGTTCCAGATCGCCCAACCGTGAAAAACCAAACTCACTGGAAGATGGTGGCTCATGCATCGCGCCAACAATAATATCCGCTCGCCCCTGAGTCAGGGCATCCCATGAACCACCGAGTACGCCATTGATAAACTTCAGGCGCGTCACGCTATGATGTTGATAAAAAGC  >  NC_000913/2087626‑2087919

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: