Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 130,603 G→A 7.6% intergenic (+47/+38) PP_mr02 → / ← engB Pseudomon‑1 ncRNA/GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947130,6030GA 7.6% 195.9 / 4.2 68intergenic (+47/+38)PP_mr02/engBPseudomon‑1 ncRNA/GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base A (2/3);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

TTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTGGGTGCCGAGAACAGCTGGATGGTCACGCCATCC  >  NC_002947/130458‑130749
                                                                                                                                                 |                                                                                                                                                  
ttAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGTGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgcc                                                                                                                                                  <  4:281470/148‑1 (MQ=255)
ttAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgccg                                                                                                                                                 >  8:265759/1‑149 (MQ=255)
ttAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgccg                                                                                                                                                 >  1:75526/1‑149 (MQ=255)
  aaGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCa                                                                                                                                               >  2:138654/1‑149 (MQ=255)
      gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACTCCGCCGCATTcg                                                                                                                                           >  3:58336/1‑149 (MQ=255)
       ttGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTcgc                                                                                                                                          <  8:162846/149‑1 (MQ=255)
          gCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCggg                                                                                                                                       >  1:186770/1‑149 (MQ=255)
           cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGt                                                                                                                                      <  1:15925/149‑1 (MQ=255)
            aGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGt                                                                                                                                      >  6:2023/1‑148 (MQ=255)
            aGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCACCGCATTCGCGGGt                                                                                                                                      >  1:242306/1‑148 (MQ=255)
              atatCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCgg                                                                                                                                        >  3:221122/1‑144 (MQ=255)
              atatCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCgg                                                                                                                                        <  4:221122/144‑1 (MQ=255)
               tatCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTaaa                                                                                                                                   <  4:309037/148‑1 (MQ=255)
               tatCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                                                                                                                  <  1:283203/149‑1 (MQ=255)
                 tCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                                                                                                >  4:112046/1‑149 (MQ=255)
                   tGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                                                                                                              >  7:26035/1‑149 (MQ=255)
                    ggACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTATCGCGAATGCGATTTCGAAAACACCGCCGCATACCCGGGTAACCCCgg                                                                                                                              >  6:240275/1‑148 (MQ=255)
                     gACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                                                                                                              <  5:2023/147‑1 (MQ=255)
                     gACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGt                                                                                                                            >  7:159759/1‑149 (MQ=255)
                     gACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCACCGCATTCGCGGGTAAACCCggg                                                                                                                             <  2:242306/148‑1 (MQ=255)
                      aCCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACGTACCCCAACCCCTACGAAGTGCGGGGGTTTTTTTTCCCGGGAATTTCCGGGAGCCGGGTTTAAACCCAATTCCGTGTGGAATACACCCCCCCCTTCGCCGGGAAACCCGGGtg                                                                                                                           >  1:206962/1‑148 (MQ=255)
                         ccGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTACCCGCGAATGCGGCGGTGTATTCGACATCGCATTCGCGAGAAAACCCgctcccgg                                                                                                                        <  1:296926/149‑8 (MQ=18)
                          cGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                                                                                                       <  1:13900/149‑1 (MQ=255)
                           gCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcg                                                                                                                      <  2:206962/149‑1 (MQ=255)
                             ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTc                                                                                                                           >  3:88605/1‑142 (MQ=255)
                             ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTc                                                                                                                           <  4:88605/142‑1 (MQ=255)
                               tAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTACCCGCGAATGCGGCGGTGTATTCGACATCGCATTCGCGAGAAAACCCg                                                                                                                               <  8:286270/136‑1 (MQ=18)
                               tAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTACCCGCGAATGCGGCGGTGTATTCGACATCGCATTCGCGAGAAAACCCg                                                                                                                               >  7:286270/1‑136 (MQ=18)
                                 gCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAAc                                                                                                                 <  8:159759/148‑1 (MQ=255)
                                  cGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACTCCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCa                                                                                                               <  4:58336/149‑1 (MQ=255)
                                       cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAAc                                                                                                                 <  1:162454/142‑1 (MQ=255)
                                       cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAAc                                                                                                                 >  2:162454/1‑142 (MQ=255)
                                        ggACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAAc                                                                                                         <  3:185506/149‑1 (MQ=255)
                                        ggACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAAc                                                                                                         >  3:104058/1‑149 (MQ=255)
                                              ggACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCgg                                                                                                   >  5:265557/1‑149 (MQ=255)
                                               gACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCgg                                                                                                   <  2:186770/148‑1 (MQ=255)
                                               gACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGc                                                                                                  <  8:26035/149‑1 (MQ=255)
                                                    cGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc                                                                                             >  8:255566/1‑149 (MQ=255)
                                                         ttCCCCCTCCCCATACGAAGTCCTGGGTTTTTTTTTCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCa                                                                                        <  5:240275/149‑1 (MQ=255)
                                                         ttCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTc                                                                                          >  8:135843/1‑147 (MQ=255)
                                                           cccccTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAg                                                                                       <  8:61737/148‑1 (MQ=255)
                                                           cccccTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGt                                                                                      >  2:173985/1‑149 (MQ=255)
                                                           cccccTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCCCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAg                                                                                       <  6:265557/148‑1 (MQ=255)
                                                                    cATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                                                                                                       >  5:28139/1‑107 (MQ=255)
                                                                    cATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                                                                                                       <  6:28139/107‑1 (MQ=255)
                                                                             tCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAg                                                                     >  7:28560/1‑148 (MQ=255)
                                                                               cGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGca                                                                   <  8:28560/148‑1 (MQ=255)
                                                                               cGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGca                                                                   >  1:276562/1‑148 (MQ=255)
                                                                                 gggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCacac                                                                >  8:12745/1‑149 (MQ=255)
                                                                                            gCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTc                                                      <  2:276562/148‑1 (MQ=255)
                                                                                             cccGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc                                                                                <  5:89171/121‑1 (MQ=255)
                                                                                             cccGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc                                                                                >  6:89171/1‑121 (MQ=255)
                                                                                               cGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTc                                                   >  1:189626/1‑148 (MQ=255)
                                                                                                       gCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCAGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCa                                           <  7:265759/148‑1 (MQ=255)
                                                                                                         cGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGgc                                        <  5:202550/149‑1 (MQ=255)
                                                                                                           ggAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGc                                                                                                  <  5:321778/89‑1 (MQ=255)
                                                                                                           ggAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGc                                                                                                  >  6:321778/1‑89 (MQ=255)
                                                                                                              gCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTg                                   >  6:85265/1‑149 (MQ=255)
                                                                                                               cGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGc                                                                        <  7:24324/111‑1 (MQ=255)
                                                                                                               cGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGc                                                                        >  8:24324/1‑111 (MQ=255)
                                                                                                               cGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTgg                                  >  4:129632/1‑149 (MQ=255)
                                                                                                                gggTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTggg                                 >  3:243827/1‑149 (MQ=255)
                                                                                                                  gTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTGGGTg                               >  4:283132/1‑149 (MQ=255)
                                                                                                                      tctcGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTGGGTGCCg                            <  2:75526/148‑1 (MQ=255)
                                                                                                                       ctcGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTGGGTGCCga                           <  7:135843/148‑1 (MQ=255)
                                                                                                                          gcgAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTGGGTGCCGAGaa                        >  7:281780/1‑148 (MQ=255)
                                                                                                                              aTGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTGGGTGCCGAGAACAGCt                   <  3:112046/149‑1 (MQ=255)
                                                                                                                                 cGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTGGGTGCCGAGAACAGCTGGa                <  3:283132/149‑1 (MQ=255)
                                                                                                                                  gATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTGGGTGCCGAGAACAGCTGGAt               >  3:57819/1‑149 (MQ=255)
                                                                                                                                      tCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTGGGTGCCGAGAACAGCTGGATGGTc           <  4:57819/149‑1 (MQ=255)
                                                                                                                                          aTACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGTAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTGGGTGCCGAGAACAGCTGGATGGTCACGc       >  7:199552/1‑149 (MQ=255)
                                                                                                                                               ccgccgCATTCGCGGGTAAACCCGGGCCAGGCGTAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCACCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTGGGTGCCGAGAACAGCTGGATGGTCACGCCATcc  <  8:199552/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTCCGCCAGCACACGGTAGGCGTCTTCCACGCCCAGGCGCTTGGGTGCCGAGAACAGCTGGATGGTCACGCCATCC  >  NC_002947/130458‑130749

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: