Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 360,391 G→A 5.5% G134S (GGC→AGC)  cynX → putative cyanate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913360,3910GA 5.5% 230.0 / 2.5 72G134S (GGC→AGC) cynXputative cyanate transporter
Reads supporting (aligned to +/- strand):  ref base G (37/31);  new base A (2/2);  total (39/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 = 9.01e-01

TACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCA  >  NC_000913/360244‑360539
                                                                                                                                                   |                                                                                                                                                    
tACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATggg                                                                                                                                                     >  2:485028/1‑149 (MQ=255)
       aaaGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGAGCGGCGGt                                                                                                                                              >  1:1370970/1‑149 (MQ=255)
       aaaGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGAGCGGCGGt                                                                                                                                              >  1:1364885/1‑149 (MQ=255)
         cggccgcTGCGGCTTCGCCGGGCACGGCTTGGGGGGGTGGGGAGCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTgg                                                                                                                                            <  1:1370087/145‑1 (MQ=255)
           tGCGCTGCTAAGCAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGcggcggcgcgt                                                                                                                                          >  2:1094217/1‑144 (MQ=255)
               ctgctgctTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTgg                                                                                                                                      <  2:271951/149‑1 (MQ=255)
                         gcagcGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATa                                                                                                                            >  1:1349725/1‑149 (MQ=255)
                           agcGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCAtcgc                                                                                                                          >  2:1026283/1‑146 (MQ=255)
                            gcgcACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACg                                                                                                                         >  2:525111/1‑149 (MQ=255)
                             cgcACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGc                                                                                                                        >  2:1343460/1‑149 (MQ=255)
                              gcACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGcc                                                                                                                       >  1:1245373/1‑149 (MQ=255)
                              gcACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTCCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGcc                                                                                                                       >  1:480819/1‑149 (MQ=255)
                                   gCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGt                                                                                                                  >  2:174293/1‑149 (MQ=255)
                                   gCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGt                                                                                                                  <  2:571963/149‑1 (MQ=255)
                                      ggggggggtggggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAg                                                                                                               <  2:1548235/145‑1 (MQ=255)
                                        gtggggtgggGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGtt                                                                                                             <  2:1296678/149‑1 (MQ=255)
                                         tggggtggggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTc                                                                                                            >  2:973934/1‑149 (MQ=255)
                                           gggtggggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGt                                                                                                              <  2:1498231/145‑1 (MQ=255)
                                           gggtggggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGt                                                                                                              >  1:1498231/1‑145 (MQ=255)
                                           gggtggggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCaa                                                                                                          >  2:1580378/1‑149 (MQ=255)
                                           gggtggggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCaa                                                                                                          <  1:1071802/149‑1 (MQ=255)
                                           gggtggggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGAGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCaa                                                                                                          <  2:1370970/149‑1 (MQ=255)
                                           gggtggggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGAGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCaa                                                                                                          <  2:1364885/149‑1 (MQ=255)
                                            ggtggggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAAc                                                                                                         <  2:1122248/149‑1 (MQ=255)
                                              tggggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACAt                                                                                                       <  1:640151/149‑1 (MQ=255)
                                               ggggAGCGGCATCATTCCGGCGGGGATGCCTTCGGGTATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGGGATGGGGCGGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATa                                                                                                      <  2:1198674/149‑1 (MQ=255)
                                                gggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAg                                                                                                     >  2:478394/1‑149 (MQ=255)
                                                 ggATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACAGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGc                                                                                                    >  1:1398376/1‑149 (MQ=255)
                                                       gCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAAcc                                                                                              >  1:490427/1‑149 (MQ=255)
                                                         atcttTCAGGCGGTGTTGCCTTCGGTGATTAACCGGCGGGTTCTGCAGCGCACGCCAGTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTg                                                                                            <  1:1326086/149‑1 (MQ=255)
                                                         atcatTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTg                                                                                            <  1:478394/149‑1 (MQ=255)
                                                         atcatTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTg                                                                                            <  1:1104514/149‑1 (MQ=255)
                                                           catTCAGGCGGTGATGCCTGCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGt                                                                                          <  2:251855/149‑1 (MQ=255)
                                                            atTCAGGCGGGGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTa                                                                                         <  2:236184/149‑1 (MQ=255)
                                                             ttCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTAt                                                                                        >  2:638353/1‑149 (MQ=255)
                                                             ttCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTAt                                                                                        >  1:232357/1‑149 (MQ=255)
                                                             ttCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGAGAAACCTGGTAt                                                                                        >  1:816031/1‑149 (MQ=255)
                                                                 ggCGGGGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCaaa                                                                                    <  1:1238200/149‑1 (MQ=255)
                                                                     gTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACt                                                                                >  1:409452/1‑149 (MQ=255)
                                                                            cTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCtgg                                                                         >  2:355373/1‑149 (MQ=255)
                                                                             tgcGGTGATTAAACGGCGGTTTCAGCAGCGCACGCGACTGGTGAGGGGGCCGGGGTCCGCGGCTCTGATGGGCGGCGGTGGGCGTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCtggt                                                                        <  1:1075621/147‑1 (MQ=255)
                                                                                    aTTAAACGGGGGTTTCAGCAGCGCACGCCACTGGTGATGTGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGGCATAACGCCGTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGTCGct                                                                 >  1:248047/1‑149 (MQ=255)
                                                                                    aTTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGct                                                                 <  1:913504/149‑1 (MQ=255)
                                                                                    aTTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGct                                                                 >  2:894989/1‑149 (MQ=255)
                                                                                     ttAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGctg                                                                >  2:421702/1‑149 (MQ=255)
                                                                                     ttAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGctg                                                                >  1:757308/1‑149 (MQ=255)
                                                                                     atAAACGGCGGATTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTATGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTCGTTAGATCAACATAGCGAAACCTGGTATCAAACACTCGCCAGGTGGGCGctg                                                                >  1:757314/2‑149 (MQ=255)
                                                                                      tAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGctgc                                                               <  1:734993/149‑1 (MQ=255)
                                                                                        aaCGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGctgcct                                                             >  2:102100/1‑149 (MQ=255)
                                                                                        aaCGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGctgcct                                                             >  2:696385/1‑149 (MQ=255)
                                                                                           ggcggTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCtgcctgcc                                                          <  2:257808/149‑1 (MQ=255)
                                                                                            gcggTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCg                                                         <  1:485028/149‑1 (MQ=255)
                                                                                                      cagcGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGct                                               >  1:827476/1‑149 (MQ=255)
                                                                                                        gcgcACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGctct                                             >  1:383542/1‑149 (MQ=255)
                                                                                                           cACGCCACTGGTGATGGGGCGGTGGTCCGCGGCTCTGATGGGCGGCGGGGGGCGTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTg                                          <  1:625948/149‑1 (MQ=255)
                                                                                                              gCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCt                                       <  1:1343460/149‑1 (MQ=255)
                                                                                                                cACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCtgg                                     <  2:480819/149‑1 (MQ=255)
                                                                                                                cACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCtgg                                     <  2:1245373/149‑1 (MQ=255)
                                                                                                                   tGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCtggtgg                                  >  2:1597794/1‑149 (MQ=255)
                                                                                                                     gTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCtggtggtg                                >  2:875694/1‑149 (MQ=255)
                                                                                                                      ggaaggggCTGGGGGCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCtggtggtgg                               <  2:86805/145‑1 (MQ=255)
                                                                                                                                   gTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCcgcg                  >  1:1322501/1‑149 (MQ=255)
                                                                                                                                    tCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGa                 >  1:128164/1‑149 (MQ=255)
                                                                                                                                     gcgcgGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAAGGCCCGGGTTAGTTCAACATAGCGAAACCTGGTCTCAAACACTCGGCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCTATGCGCCCGCGAg                <  2:757314/148‑1 (MQ=255)
                                                                                                                                     ccGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAg                <  2:757308/149‑1 (MQ=255)
                                                                                                                                        gggCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTc             <  2:1349725/148‑1 (MQ=255)
                                                                                                                                          gCTCTGATCGCCGGCGGTGGGCTTGGTTCCGCCATAACGCCCTGGTTAGATCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGTCAAAGCGCCCGCGAGGTCGc           <  1:777296/149‑1 (MQ=255)
                                                                                                                                           ctctGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGcc          >  2:1087661/1‑149 (MQ=255)
                                                                                                                                            tctGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCct         >  2:801653/1‑149 (MQ=255)
                                                                                                                                              tGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCctct       >  2:1353904/1‑149 (MQ=255)
                                                                                                                                               gATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCtt      >  1:1027397/1‑149 (MQ=255)
                                                                                                                                               gATGGGCGGCGGTGGGCTTGGGGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCtt      <  2:794259/149‑1 (MQ=255)
                                                                                                                                                 tGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAACGCGCCCGCGAGGTCGCCTCTTcc    >  1:832421/1‑149 (MQ=255)
                                                                                                                                                   ggcggcGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCca  >  1:623264/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TACCCGCAAAGTGCGCTGCTGCTTAGCAGCGCACTGCTTGGTGGGGTGGGGATCGGCATCATTCAGGCGGTGATGCCTTCGGTGATTAAACGGCGGTTTCAGCAGCGCACGCCACTGGTGATGGGGCTGTGGTCCGCGGCTCTGATGGGCGGCGGTGGGCTTGGTGCCGCCATAACGCCCTGGTTAGTTCAACATAGCGAAACCTGGTATCAAACACTCGCCTGGTGGGCGCTGCCTGCCGTTGTTGCGCTCTTTGCCTGGTGGTGGCAAAGCGCCCGCGAGGTCGCCTCTTCCCA  >  NC_000913/360244‑360539

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