Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 603,393 T→A 5.2% 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,3930TA 5.2% 319.8 / 2.8 128intergenic (+58/+23)pheP/ybdGphenylalanine transporter/mechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base T (59/58);  new base A (3/5);  total (65/63)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 132 total aligned reads displayed.

CACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCATA  >  NC_000913/603245‑603541
                                                                                                                                                    |                                                                                                                                                    
cACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAt                                                                                                                                                      >  2:778625/1‑149 (MQ=255)
  agatggatgATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGATATAAGGCATTCTAGCTACATCCGAGAAAACGTTGTCAACCATCCGAAACCGCTCTCATCCATTc                                                                                                                                                    <  2:1337665/148‑1 (MQ=255)
     ttgatgatTTGCGATTGTCAGAGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGTCATTTAAAACGCTTCTTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCCTCCATTCGAt                                                                                                                                                 <  1:1008951/147‑1 (MQ=255)
        atgatATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATgag                                                                                                                                              >  1:1275708/1‑149 (MQ=255)
            tatGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATTTCAACCATACGAAACCGCTCTCATCCATTCGATgagtgag                                                                                                                                          >  2:1040043/1‑145 (MQ=255)
            tatGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCg                                                                                                                                          <  1:290025/149‑1 (MQ=255)
            tatGCGCTTGTCACCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCAAGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCg                                                                                                                                          >  1:1014307/1‑149 (MQ=255)
             atGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCgg                                                                                                                                         >  2:521239/1‑149 (MQ=255)
             atGCGCTTGTCAGCGATCCTGCTGCCGGGGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCGGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCGAATGGATGAGAGCgg                                                                                                                                         <  2:170244/149‑1 (MQ=255)
               gcgcTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtt                                                                                                                                       >  1:1535538/1‑149 (MQ=255)
                 gcTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGtttt                                                                                                                                     >  1:1247350/1‑149 (MQ=255)
                 gcTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTAGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGttat                                                                                                                                     >  1:229132/1‑147 (MQ=255)
                      tCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTaa                                                                                                                                <  1:824838/149‑1 (MQ=255)
                        aGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACTCTACATCCGACAAAACGATGTCACCCCTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAtt                                                                                                                              >  2:857873/1‑149 (MQ=255)
                        aGCGATCCTGCTGCCGGTGTCGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCCTCCATTCGATGATAGCCGTTTTTTTAAtt                                                                                                                              >  1:1286219/1‑149 (MQ=255)
                         gCGATCCTGCTGCCGGTGTGGATTGTATTCCTGATTATGGAATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTa                                                                                                                             <  2:696872/149‑1 (MQ=255)
                           gATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTTAAACGCAGCGTCGTATATAAGGCATTCGCGCTACATCCGACAAAACGATGGAAACCATCCGAAACCGCGCTGCACCCGTCGCTGAGAGCGGGTTATTTAATTCCt                                                                                                                           >  2:1011201/1‑149 (MQ=255)
                           gATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACt                                                                                                                           >  1:1568992/1‑149 (MQ=255)
                               ctgctgCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAATACTGCtt                                                                                                                       >  1:467941/1‑149 (MQ=255)
                                   tgcCGGTGTGGATTGTATTCCTTTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAAt                                                                                                                   <  2:697189/149‑1 (MQ=255)
                                   tgcCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACTCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAAt                                                                                                                   <  1:229628/149‑1 (MQ=255)
                                    gcCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATg                                                                                                                  >  1:1344050/1‑149 (MQ=255)
                                     ccGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGc                                                                                                                 <  1:198454/149‑1 (MQ=255)
                                      gggTGTGGATTGTTTTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCa                                                                                                                <  2:1527175/148‑1 (MQ=255)
                                          gtgGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCTCTCTCATACATTAGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCg                                                                                                            >  1:1039782/1‑149 (MQ=255)
                                            ggATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGcc                                                                                                          >  2:590770/1‑149 (MQ=255)
                                                 gTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCAGTCGATGAGCGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagagc                                                                                                     >  1:1323650/1‑148 (MQ=255)
                                                  tATTCCTGTTTATTGCATTTAAAACGCTGCGTCGTAAATAAGGCATTCACGCTACATCCGACAAAACTATGTCAACCATCCGAAACCTCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagagag                                                                                                    <  1:240209/149‑1 (MQ=255)
                                                  tATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCagagag                                                                                                    <  2:502960/149‑1 (MQ=255)
                                                   aTTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGc                                                                                                   <  1:297994/149‑1 (MQ=255)
                                                     tCCTGTTTATGTCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCTATTAGAGCGGTTTTTTTAATTACTGCTAAAATGCCCCCGCCAgagagaga                                                                                                 >  2:1039529/1‑149 (MQ=255)
                                                      ccTGTTTATGGCATTTAACACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACCAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGGTTTTTTAATTACTGCTTAAATGCCCCCGCCAGAGAGCGaa                                                                                                >  1:140659/1‑149 (MQ=255)
                                                       cTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAat                                                                                               >  1:116671/1‑149 (MQ=255)
                                                        tGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAata                                                                                              >  2:1453564/1‑149 (MQ=255)
                                                         gTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAatat                                                                                             >  2:780465/1‑149 (MQ=255)
                                                           ttATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCa                                                                                           <  2:269738/149‑1 (MQ=255)
                                                            tATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCTACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCAt                                                                                          >  1:1305017/1‑149 (MQ=255)
                                                                gCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGcc                                                                                      >  2:1334681/1‑149 (MQ=255)
                                                                 cATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCg                                                                                     <  1:1012843/149‑1 (MQ=255)
                                                                    ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                                                                                  <  2:658371/149‑1 (MQ=255)
                                                                    ttAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGtt                                                                                  <  2:1530038/149‑1 (MQ=255)
                                                                      aaacCGCTGCGTCGGAAATCCGGCATGCCCGCTACATCCGACAAAACGATGTCATCCATCCGAAACCGCGCTCATCCATTCGATGAGAGCGGTTTTTTTAATAACTGCTTAAATGCACCCGCCAGAGAGCGACTATCATTGCCGGTTgg                                                                                <  2:1454801/149‑1 (MQ=255)
                                                                      aaaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTTCTTAAATGCACACGCCAGAGAGCGAATATCATTGCCGGTTgg                                                                                >  2:860105/1‑149 (MQ=255)
                                                                       aaaCGCTTCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGc                                                                               <  2:123646/149‑1 (MQ=255)
                                                                       aaaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGAGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCAGGTTGGc                                                                               >  1:144221/1‑149 (MQ=255)
                                                                        aaCGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATAGCCGTTTGGcc                                                                              >  1:1080025/1‑148 (MQ=255)
                                                                          cGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCCTCCGAAACCTCTCACATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGGTGGCGAc                                                                            <  2:1020429/149‑1 (MQ=255)
                                                                            cTGCGTCGGAAATAAGGCATTCACTCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCACAGCCATTCGATGTGAGCGGTTTTTTTAAATACTGCTTAAATGCACCCGCCAGAGAGCGAATTTCATTTCCGGTTGGCGACtg                                                                          >  1:144074/1‑149 (MQ=255)
                                                                            cTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtg                                                                          <  1:778625/149‑1 (MQ=255)
                                                                             tGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtga                                                                         <  1:165241/149‑1 (MQ=255)
                                                                               ggTCGGAAATTACGCAGTCCCCCTTCCTCAGACAAACCGATGTCAAACATCCGAAACCGCTCTCATCCATTCGATGAGAGCGTTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACtgatg                                                                       <  1:1011201/148‑1 (MQ=255)
                                                                                 tCGGAAATAAGGCATTAACGCTACATCCGACAAAACGATGTCAACCAACCGAAACCGCAATCATCCATTCCAAGAGAGCGGTTTATTTAATTACTGCTAATATGCACCCGCCAGAGCCCGAATATCATTGCCGGTTGGCGAGTGATGaa                                                                     >  1:931900/1‑149 (MQ=255)
                                                                                  cGGAAATAAGGAATTCACGCTACATCCGACAAACCGATGTCACCCATCCGAAACCGATCTCATCGATTCGATGAGAGCGGTTATTTTACTTTCGGCTTAAATGCACCAGCCATAGAGCGAATATCAATGACGGTTGGCGACTGAGGAAg                                                                    >  2:1037067/1‑149 (MQ=255)
                                                                                   ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGt                                                                   <  1:375195/149‑1 (MQ=255)
                                                                                   ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGt                                                                   <  1:375197/149‑1 (MQ=255)
                                                                                   ggAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGt                                                                   >  2:129888/1‑149 (MQ=255)
                                                                                       aTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                                                               >  2:130113/1‑149 (MQ=255)
                                                                                       aTAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACAGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCa                                                               <  2:939980/149‑1 (MQ=255)
                                                                                         atggCATTCACGCTTCATACTAGAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGa                                                             <  2:674757/147‑1 (MQ=255)
                                                                                         aaGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGa                                                             <  1:780465/149‑1 (MQ=255)
                                                                                              aTTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaaa                                                        >  1:85603/1‑149 (MQ=255)
                                                                                              aTTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaaa                                                        >  2:797540/1‑149 (MQ=255)
                                                                                              aTTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCaaa                                                        >  2:1291096/1‑149 (MQ=255)
                                                                                                  aCGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAActct                                                    >  2:600969/1‑149 (MQ=255)
                                                                                                      tACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGa                                                <  1:431260/149‑1 (MQ=255)
                                                                                                      tACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGa                                                <  1:1334681/149‑1 (MQ=255)
                                                                                                      tACATACGACATAACGTTGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGa                                                <  2:1039782/149‑1 (MQ=255)
                                                                                                           ccGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATc                                           >  2:373938/1‑149 (MQ=255)
                                                                                                            cGACAAAACGATGTCACCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGCCCAAACTCTTCGACAATCg                                          >  2:693375/1‑149 (MQ=255)
                                                                                                            cGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCACAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCg                                          >  2:655377/1‑149 (MQ=255)
                                                                                                                aaaaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTGTTTTTAATTACTGCTTAAAAGCCCCCGCCAGAGCGCGAAGATCATTGCCGGTTGGCGACTGATGATGTCGCAGACCAAACTCTTCGACAATCGCaaa                                      >  2:1019118/1‑149 (MQ=255)
                                                                                                                  aaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGCCCAAACTCTTCGACAATCGCAAata                                    >  1:1039733/1‑149 (MQ=255)
                                                                                                                  aaCGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAATACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACACTACGACAATCGCAAata                                    >  2:993197/1‑149 (MQ=255)
                                                                                                                    cGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATAtg                                  >  1:988059/1‑149 (MQ=255)
                                                                                                                        gTCAACCATCCGAAACCGCTCTCATCCATTCGATGATAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTTATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGAt                              <  1:1039529/149‑1 (MQ=255)
                                                                                                                          cAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCg                            >  1:812345/1‑149 (MQ=255)
                                                                                                                            aCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGaa                          <  1:1147802/149‑1 (MQ=255)
                                                                                                                            aCCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGaa                          <  1:1453564/149‑1 (MQ=255)
                                                                                                                             ccATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAat                         <  2:1535538/149‑1 (MQ=255)
                                                                                                                             ccATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAAAGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAat                         <  1:797359/149‑1 (MQ=255)
                                                                                                                             caaTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAat                         <  2:1336121/147‑1 (MQ=255)
                                                                                                                              cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAata                        >  2:110780/1‑149 (MQ=255)
                                                                                                                              cATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAata                        >  2:959839/1‑149 (MQ=255)
                                                                                                                               aTCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAatat                       >  2:501731/1‑149 (MQ=255)
                                                                                                                                 ccGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAatatat                     >  1:1417896/1‑149 (MQ=255)
                                                                                                                                 ccGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAAATACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAatatat                     <  2:467941/149‑1 (MQ=255)
                                                                                                                                    aaaCCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAg                  >  1:47699/1‑149 (MQ=255)
                                                                                                                                         gCTCTCATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCtta             >  2:1181009/1‑147 (MQ=255)
                                                                                                                                         gCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGa             <  1:521239/149‑1 (MQ=255)
                                                                                                                                          ctctCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGCCCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGaa            >  1:900594/1‑149 (MQ=255)
                                                                                                                                          ctctCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGaa            <  1:567828/149‑1 (MQ=255)
                                                                                                                                           tctcATCGAATGGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAAt           >  1:1087704/1‑149 (MQ=255)
                                                                                                                                            ctcATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGGAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATg          <  1:115011/149‑1 (MQ=255)
                                                                                                                                              ccaccATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCt        <  2:286575/146‑1 (MQ=255)
                                                                                                                                               aTCCATTCGATGAGTGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGAGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCtt       <  1:400275/149‑1 (MQ=255)
                                                                                                                                               aTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCtt       >  1:322667/1‑149 (MQ=255)
                                                                                                                                               aTCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGAAAATCGCAAATATGTGATCTAATATATCAGCCTGAATGCtt       >  2:914470/1‑149 (MQ=255)
                                                                                                                                                   aTTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCat   >  1:1061974/1‑149 (MQ=255)
                                                                                                                                                    ttCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCata  <  1:600969/149‑1 (MQ=255)
                                                                                                                                                    ttCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCata  <  2:1417896/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACGATGGATGATATGCGCTTGTCAGCGATCCTGCTGCCGGTGTGGATTGTATTCCTGTTTATGGCATTTAAAACGCTGCGTCGGAAATAAGGCATTCACGCTACATCCGACAAAACGATGTCAACCATCCGAAACCGCTCTCATCCATTCGATGAGAGCGGTTTTTTTAATTACTGCTTAAATGCACCCGCCAGAGAGCGAATATCATTGCCGGTTGGCGACTGATGAAGTCGCAGACCAAACTCTTCGACAATCGCAAATATGTGATCGAATATATCAGCCTGAATGCTTTCATA  >  NC_000913/603245‑603541

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