Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,777,379 A→G 5.7% T89A (ACT→GCT)  yfjZ → CP4‑57 prophage; antitoxin of the YpjF‑YfjZ toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,777,3790AG 5.7% 215.5 / 2.9 71T89A (ACT→GCT) yfjZCP4‑57 prophage; antitoxin of the YpjF‑YfjZ toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (36/30);  new base G (2/2);  total (38/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 = 8.69e-01

GTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTA  >  NC_000913/2777231‑2777527
                                                                                                                                                    |                                                                                                                                                    
gTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATa                                                                                                                                                      <  1:10347/149‑1 (MQ=255)
     gTGACGCCGAATGTCCTAAGCTGGATGTGGTTTTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCtt                                                                                                                                                 <  7:96722/149‑1 (MQ=255)
     gTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCtt                                                                                                                                                 <  8:30933/149‑1 (MQ=255)
        aCGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGt                                                                                                                                              <  5:109029/149‑1 (MQ=255)
        aCGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGt                                                                                                                                              <  2:99956/149‑1 (MQ=255)
         cGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTa                                                                                                                                             <  7:273653/149‑1 (MQ=255)
                 gTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGGCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCgg                                                                                                                                      <  2:224842/148‑1 (MQ=255)
                           ggATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACa                                                                                                                           >  4:120718/1‑149 (MQ=255)
                            gATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTCCCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACAt                                                                                                                          >  2:137360/1‑149 (MQ=255)
                                     tttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGttt                                                                                                                  <  1:254967/148‑1 (MQ=255)
                                     tttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTa                                                                                                                 >  6:91989/1‑149 (MQ=255)
                                     tttCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGTAGTTGCGGCTACGTATACATCGCTGTTTAc                                                                                                                >  2:118317/1‑149 (MQ=255)
                                       tCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTAc                                                                                                                <  4:15778/148‑1 (MQ=255)
                                             ttttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACt                                                                                                          >  8:217459/1‑148 (MQ=255)
                                               ttATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCaa                                                                                                       >  7:69142/1‑149 (MQ=255)
                                                 aTCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCACCg                                                                                                     >  5:153780/1‑149 (MQ=255)
                                                  tCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATGCTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGc                                                                                                    >  5:300598/1‑149 (MQ=255)
                                                  tCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATGCTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGc                                                                                                    >  7:173447/1‑149 (MQ=255)
                                                    aGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATCCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGct                                                                                                   <  2:97548/148‑1 (MQ=255)
                                                     gCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaa                                                                                                 <  7:117816/149‑1 (MQ=255)
                                                     gCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaa                                                                                                 >  6:59248/1‑149 (MQ=255)
                                                        aGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGctaact                                                                                               >  2:252592/1‑148 (MQ=255)
                                                          aTAGAGTCGATGCTGCCCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAAt                                                                                            >  5:220008/1‑149 (MQ=255)
                                                             gagTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTc                                                                                         <  1:252592/149‑1 (MQ=255)
                                                             gagTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTGTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACCTCGCAGTTTACCCCACTCCACGCTAACGAATTTc                                                                                         >  2:41500/1‑149 (MQ=255)
                                                               gTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCAc                                                                                       <  7:158088/149‑1 (MQ=255)
                                                               gTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCAc                                                                                       >  3:206718/1‑149 (MQ=255)
                                                                tCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGtt                                                                                                                   >  3:278797/1‑120 (MQ=255)
                                                                tCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGtt                                                                                                                   <  4:278797/120‑1 (MQ=255)
                                                                  gaagCTCACCACTGGTGAACTGAATCCCCGCCATGCCCAAAGCGTCACCCTGTACCACAACGGTTTTACCTGCCAAGCCGATACTCTTGCTAGTTGCGGCTACCTATACATCGCTGATTACACCACTCAACCCTAAATAATTTCAc                                                                                       >  8:229351/4‑146 (MQ=255)
                                                                       tGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGagc                                                                                 >  6:24991/1‑147 (MQ=255)
                                                                              tGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAAc                                                                                                      <  2:210299/119‑1 (MQ=255)
                                                                              tGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAAc                                                                                                      >  1:210299/1‑119 (MQ=255)
                                                                              tGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACCTATACATCGCTGTTTACCCCACTCACCGCTAACTAATTTCACGAGAGCAAGCATGaa                                                                        >  7:276967/1‑149 (MQ=255)
                                                                                gTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACATCGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCTCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAaca                                                                      <  5:215947/149‑1 (MQ=255)
                                                                                gTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAaca                                                                      <  8:184480/149‑1 (MQ=255)
                                                                                 tGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAacac                                                                     >  3:296742/1‑149 (MQ=255)
                                                                                 tGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAacac                                                                     <  5:91989/149‑1 (MQ=255)
                                                                                    aCTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACActc                                                                   >  3:3372/1‑148 (MQ=255)
                                                                                     cTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTa                                                                 <  7:217459/149‑1 (MQ=255)
                                                                                      tGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTATCTAATTTCACGAGCGCAAGCATGAACACTCTAc                                                                >  4:284831/1‑149 (MQ=255)
                                                                                          tCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGc                                                            >  5:233852/1‑149 (MQ=255)
                                                                                            cccGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTa                                                          >  5:165493/1‑149 (MQ=255)
                                                                                                    gCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCg                                                  >  8:210646/1‑149 (MQ=255)
                                                                                                    gCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCg                                                  >  8:204234/1‑149 (MQ=255)
                                                                                                      ccAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATGCTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCa                                                <  8:173447/149‑1 (MQ=255)
                                                                                                      ccAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATGCTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCa                                                <  6:300598/149‑1 (MQ=255)
                                                                                                         aTGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTG‑‑‑CAATTTCGCAggcggc                                          <  5:24991/149‑1 (MQ=255)
                                                                                                          tGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACTCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAggcg                                            >  7:56232/1‑149 (MQ=255)
                                                                                                            cGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAggcgg                                           <  7:229351/148‑1 (MQ=255)
                                                                                                             gTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAggcggc                                          >  1:122668/1‑148 (MQ=255)
                                                                                                             gTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGgcggcg                                         <  4:296742/149‑1 (MQ=255)
                                                                                                             gTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGCTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGgcggcg                                         <  2:15113/149‑1 (MQ=255)
                                                                                                               cACCATGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTTGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGaa                                       <  1:128728/149‑1 (MQ=255)
                                                                                                                          acaACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGt                            <  7:204234/149‑1 (MQ=255)
                                                                                                                             aCGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACActc                                                                   >  1:160023/1‑107 (MQ=255)
                                                                                                                             aCGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACActc                                                                   <  2:160023/107‑1 (MQ=255)
                                                                                                                                gTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAg                      >  2:250598/1‑149 (MQ=255)
                                                                                                                                 gcccaCCTGCCACCCCGATCCTCTCGGGAGGTGCGGCTCCGTATACCTCGCTGCTTACCCCACTCAACCCTACCTAATTTCCCGAGAGCAAGCATGAACCCTCTACCTGCTACAATATCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGt                     <  4:206718/145‑1 (MQ=255)
                                                                                                                                     aCCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCt                 >  7:175462/1‑149 (MQ=255)
                                                                                                                                     aCCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCt                 >  3:262624/1‑149 (MQ=255)
                                                                                                                                      ccTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCACTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTg                >  8:295177/1‑149 (MQ=255)
                                                                                                                                       cTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGt               >  3:32794/1‑149 (MQ=255)
                                                                                                                                          cGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTc              <  4:101757/147‑1 (MQ=255)
                                                                                                                                           gAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCACTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTgg           <  7:295177/149‑1 (MQ=255)
                                                                                                                                            aaGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGc          <  5:59248/149‑1 (MQ=255)
                                                                                                                                             aGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTg                    <  8:1281/138‑1 (MQ=255)
                                                                                                                                             aGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTAACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTg                    >  7:1281/1‑138 (MQ=255)
                                                                                                                                                 gATACTCTTGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGc    <  1:118317/149‑1 (MQ=255)
                                                                                                                                                  aTACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGc    >  1:214856/1‑149 (MQ=255)
                                                                                                                                                   tACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGc    >  1:142614/1‑148 (MQ=255)
                                                                                                                                                   tACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGc    >  1:244222/1‑148 (MQ=255)
                                                                                                                                                    aCTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTa  <  3:127300/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTTTAGTGACGCCGAATGTCCTAAGCTGGATGTGGTATTTCCACATTTTATCAGCCAGATAGAGTCGATGCTGACCACTGGTGAACTGAATCCCCGCCATGCCCAATGCGTCACCCTGTACCACAACGGTTTTACCTGCGAAGCCGATACTCTTGGTAGTTGCGGCTACGTATACATCGCTGTTTACCCCACTCAACGCTAACTAATTTCACGAGAGCAAGCATGAACACTCTACCTGCTACAATTTCGCAGGCGGCGAAGCCCTGCCTGTCGCCAGTGGCTGTCTGGCAAATGCTA  >  NC_000913/2777231‑2777527

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: