Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 281,430 T→A 6.5% E185V (GAG→GTG)  yagA ← CP4‑6 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,4300TA 6.5% 193.3 / 3.4 61E185V (GAG→GTG) yagACP4‑6 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (31/26);  new base A (2/2);  total (33/28)
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.64e-01

CCCACCCGGATACCCAGGCGCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTC  >  NC_000913/281283‑281575
                                                                                                                                                   |                                                                                                                                                 
cccACCCGGATACCCAGGCGCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGtctc                                                                                                                                                  >  6:188079/1‑149 (MQ=255)
   aCCCGGATACCCAGGCGCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTcgcg                                                                                                                                               >  6:223631/1‑149 (MQ=255)
      cGGATACCCAGGCGCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTACCGCTCAATCACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCg                                                                                                                                            >  8:94986/1‑149 (MQ=255)
        gATACCCAGGCGCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCCCGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGtt                                                                                                                                          >  1:178514/1‑149 (MQ=255)
           aCCCAGGCGCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCAt                                                                                                                                       >  8:36560/1‑149 (MQ=255)
             ccAGGCGCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCg                                                                                                                                     >  5:86949/1‑149 (MQ=255)
                ggCGCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTa                                                                                                                                  <  5:98486/149‑1 (MQ=255)
                 gcgcATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTa                                                                                                                                  >  4:82352/1‑148 (MQ=255)
                    cATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCAGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGt                                                                                                                              >  5:22839/1‑149 (MQ=255)
                    cATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGGGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGt                                                                                                                              <  2:40854/149‑1 (MQ=255)
                      tCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACGGTgc                                                                                                                            >  5:155230/1‑149 (MQ=255)
                         gCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCTGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGcc                                                                                                                         <  6:22839/149‑1 (MQ=255)
                          ccACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTgcgc                                                                                                                          >  5:41066/1‑147 (MQ=255)
                            acaGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAgg                                                                                                                      <  5:188079/149‑1 (MQ=255)
                              aGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCACGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGca                                                                                                                    >  6:206067/1‑149 (MQ=255)
                              aGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCACGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGca                                                                                                                    <  5:206067/149‑1 (MQ=255)
                                 tCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTc                                                                                                                                         <  6:234212/125‑1 (MQ=255)
                                 tCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTc                                                                                                                                         >  5:234212/1‑125 (MQ=255)
                                 tCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAg                                                                                                                 >  1:132590/1‑149 (MQ=255)
                                  ccAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTcgc                                                                                                                                                >  6:15889/1‑117 (MQ=255)
                                  ccAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTcgc                                                                                                                                                <  5:15889/117‑1 (MQ=255)
                                  ccAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAgg                                                                                                                >  2:47195/1‑149 (MQ=255)
                                  ccAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAgg                                                                                                                >  3:128772/1‑149 (MQ=255)
                                   cAGCGCCGTCCAGGTGCCGGTGGTGGCGCCCCACGGGGAGCCGTTATCCATGGTCATCCGGGCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAggg                                                                                                               <  5:7179/149‑1 (MQ=255)
                                     gcgcCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGa                                                                                                              >  1:248893/1‑148 (MQ=255)
                                       gcCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGaaa                                                                                                            >  2:199347/1‑148 (MQ=255)
                                        ccGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGaaaaa                                                                                                          <  4:218133/149‑1 (MQ=255)
                                           tCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACgg                                                                                                       <  6:41066/149‑1 (MQ=255)
                                                 tGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTgg                                                                                                 <  5:223631/149‑1 (MQ=255)
                                                  gCCGGTGTTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCAGTTCCGGCAGGCCTTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGgt                                                                                                >  4:72597/1‑149 (MQ=255)
                                                     ggtggtGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGgtcgt                                                                                             <  1:29554/149‑1 (MQ=255)
                                                      gtggtgTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGgtcgtc                                                                                            <  1:199347/149‑1 (MQ=255)
                                                            tCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCcagcag                                                                                      >  5:180676/1‑149 (MQ=255)
                                                              gCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAgg                                                                                     <  6:95977/148‑1 (MQ=255)
                                                               ccccACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGa                                                                                                              <  3:5564/122‑1 (MQ=255)
                                                               ccccACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGa                                                                                                              >  4:5564/1‑122 (MQ=255)
                                                               ccccACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGt                                                                                   <  2:248893/149‑1 (MQ=255)
                                                                    cGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCg                                                                              <  3:82352/149‑1 (MQ=255)
                                                                     ggTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCACGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCg                                                                              >  1:104116/1‑148 (MQ=255)
                                                                     ggTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCACGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCg                                                                              >  1:100825/1‑148 (MQ=255)
                                                                           ccGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGATGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCCCAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCa                                                                       >  1:48946/1‑149 (MQ=255)
                                                                             gTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTa                                                                                                                                  <  6:150932/88‑1 (MQ=255)
                                                                              ttATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCACGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCg                                                                    <  2:104116/149‑1 (MQ=255)
                                                                              ttATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCACGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCg                                                                    <  2:100825/149‑1 (MQ=255)
                                                                                   cATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACCGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCAc                                                               <  6:155230/149‑1 (MQ=255)
                                                                                    aTGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACa                                                                                                                                >  7:172512/1‑83 (MQ=255)
                                                                                    aTGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACa                                                                                                                                <  8:172512/83‑1 (MQ=255)
                                                                                    aTGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCAcc                                                              >  2:98781/1‑149 (MQ=255)
                                                                                    aTGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCAcc                                                              >  6:103412/1‑149 (MQ=255)
                                                                                    aTGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGGTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCAcc                                                              >  7:205792/1‑149 (MQ=255)
                                                                                        tCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGcc                                                           >  1:214953/1‑148 (MQ=255)
                                                                                           tccggtccggCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCaaaa                                                       <  2:132590/149‑1 (MQ=255)
                                                                                                              aCGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAAt                                    <  1:47195/149‑1 (MQ=255)
                                                                                                                gCTCAAACACGCTGACCAGCTGCTGCTGCCCGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATcc                                  <  8:222196/149‑1 (MQ=255)
                                                                                                                  tCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCAt                                >  5:67630/1‑149 (MQ=255)
                                                                                                                    aaaCACGCTGACCAGCTGCTGATGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCt                              <  2:48946/149‑1 (MQ=255)
                                                                                                                      acacGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGc                            >  5:80071/1‑149 (MQ=255)
                                                                                                                       cacGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGcc                           <  3:72597/149‑1 (MQ=255)
                                                                                                                          gCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCa                          >  8:232267/1‑147 (MQ=255)
                                                                                                                            tGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAgg                      >  5:177326/1‑149 (MQ=255)
                                                                                                                              aCCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACTGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCg                    >  8:217545/1‑149 (MQ=255)
                                                                                                                                   ctgctgctgCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGTGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCg               <  7:36560/149‑1 (MQ=255)
                                                                                                                                                 tctcGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTc  >  7:194850/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCCACCCGGATACCCAGGCGCATCAGCCACAGCTCCAGCGCCGTCCAGGTGCCGGTGGTGTCGCCCCACGGTGAGCCGTTATCCATGGTCATCCGGTCCGGCAGGCCGTAACGCTCAAACACGCTGACCAGCTGCTGCTGCACGGTCTCGCGCCGTTCATCGGTACAGTGCGCCAGGCACAGGGAAAAACGGGAGTGGTCGTCCAGCAGGGTGAGCGGATGGCAGCGTCCACCGCCAAAAGGAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTC  >  NC_000913/281283‑281575

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