Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 985,937 G→C 5.9% L233F (TTG→TTC pilF → type IV pili biogenesis protein PilF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947985,9370GC 5.9% 213.4 / 4.0 68L233F (TTG→TTCpilFtype IV pili biogenesis protein PilF
Reads supporting (aligned to +/- strand):  ref base G (30/34);  new base C (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

TGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAACGG  >  NC_002947/985790‑986083
                                                                                                                                                   |                                                                                                                                                  
tGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTgg                                                                                                                                                   >  8:316496/1‑149 (MQ=255)
       tCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCCAATTCGGCCAgtg                                                                                                                                            <  6:37949/149‑3 (MQ=255)
          taagAAAACAGGCATTATGTGCCGGCTCGTGACTACTACTATCGTTTGAGCCAGTTGAGGGAACACGACGGCCGTAGCCTGCTGCTGGTCAGCCGCCTTGCCAGTGTGTTCGACGAGCAGGGCACACTGGCCTAATTGTGGCAGCAAtt                                                                                                                                         <  6:200168/146‑1 (MQ=255)
            cGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTa                                                                                                                                        <  6:235604/148‑1 (MQ=255)
             gAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACCCCCCTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGCGAGTTCGACGAGCAGGGCACACTGGCCGAATTGCGCCAGCAATTAc                                                                                                                                       >  7:332375/1‑148 (MQ=255)
                aaCAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAAc                                                                                                                                    >  5:200270/1‑148 (MQ=255)
                      cATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTa                                                                                                                             >  7:121109/1‑149 (MQ=255)
                             tgCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCgg                                                                                                                       >  2:99324/1‑148 (MQ=255)
                              gCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTa                                                                                                                     >  4:227468/1‑149 (MQ=255)
                                cGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACg                                                                                                                   >  4:9944/1‑149 (MQ=255)
                                       ggACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGaa                                                                                                             >  2:121752/1‑148 (MQ=255)
                                        gACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGc                                                                                                                                                 >  7:156828/1‑111 (MQ=255)
                                        gACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGc                                                                                                                                                 <  8:156828/111‑1 (MQ=255)
                                           tactacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCg                                                                                                                >  1:99779/1‑141 (MQ=255)
                                           tactacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAatat                                                                                                          >  4:230323/1‑147 (MQ=255)
                                           tactacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCa                                                                                                        <  6:200270/149‑1 (MQ=255)
                                            actacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGTTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCCAATTCGGCCAgtg                                                                                                                                            <  8:71337/112‑3 (MQ=255)
                                            actacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGTTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCCAATTCGGCCAgtg                                                                                                                                            >  7:71337/1‑110 (MQ=255)
                                             ctacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGTAGGGCACACTGGCCGAATTGGGCCAGCAAt                                                                                                                                          <  2:342449/113‑1 (MQ=255)
                                             ctacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGTAGGGCACACTGGCCGAATTGGGCCAGCAAt                                                                                                                                          >  1:342449/1‑113 (MQ=255)
                                             ctacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAg                                                                                                       <  6:170913/148‑1 (MQ=255)
                                             ctacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAg                                                                                                       <  7:316496/148‑1 (MQ=255)
                                              tacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCa                                                                                                        <  3:48816/146‑1 (MQ=255)
                                              tacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCa                                                                                                        >  4:48816/1‑146 (MQ=255)
                                                 gATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTg                                                                                                                                                    >  6:234603/1‑99 (MQ=255)
                                                 gATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAAt                                                                                                                                          <  7:155310/109‑1 (MQ=255)
                                                 gATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAAt                                                                                                                                          >  8:155310/1‑109 (MQ=255)
                                                 gATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGAc                                                                                                                                 <  4:276893/118‑1 (MQ=255)
                                                 gATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGAc                                                                                                                                 >  3:276893/1‑118 (MQ=255)
                                                 gATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATa                                                                                                  <  3:230323/149‑1 (MQ=255)
                                                  aTCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATAc                                                                                                 <  3:9944/149‑1 (MQ=255)
                                                  aTCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATAc                                                                                                 <  1:99324/149‑1 (MQ=255)
                                                     gTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCt                                                                                               <  8:121109/148‑1 (MQ=255)
                                                      tttCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATc                                                                                                                           <  7:91079/119‑1 (MQ=255)
                                                      tttCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATc                                                                                                                           >  8:91079/1‑119 (MQ=255)
                                                          aGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCg                                                                                           >  4:87249/1‑147 (MQ=255)
                                                           gCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCt                                                                                               <  5:219907/142‑1 (MQ=255)
                                                           gCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCt                                                                                               >  6:219907/1‑142 (MQ=255)
                                                           gCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAg                                                                                        >  1:53173/1‑149 (MQ=255)
                                                                   aGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGTTGTTCGAGGAGCAGGGCACACTGGCCGAATTGGGCGAGCAATTACAACGACTTTATCCCGGTACGCCGTAATATCAGCAATACCCGTCGGAGCAACGAt                                                                                 <  8:332375/148‑1 (MQ=255)
                                                                    gCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATg                                                                                >  6:35842/1‑148 (MQ=255)
                                                                        ccACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCTGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGaa                                                                              <  3:139141/146‑1 (MQ=255)
                                                                        ccACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATAc                                                                                                 >  8:114193/1‑127 (MQ=255)
                                                                        ccACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATAc                                                                                                 <  7:114193/127‑1 (MQ=255)
                                                                        ccACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGaa                                                                              >  4:139141/1‑146 (MQ=255)
                                                                         cACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAAcgcc                                                                          >  8:94589/1‑149 (MQ=255)
                                                                           cgacgCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGTTGAAcgccgc                                                                        >  1:156430/1‑149 (MQ=255)
                                                                               gCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACg                                                                                   <  8:195306/134‑1 (MQ=255)
                                                                               gCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACg                                                                                   >  7:195306/1‑134 (MQ=255)
                                                                                       ccTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagt                                                            <  5:35842/149‑1 (MQ=255)
                                                                                       ccTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagt                                                            <  3:227468/149‑1 (MQ=255)
                                                                                           ctgctgGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagtag                                                          <  3:87249/147‑1 (MQ=255)
                                                                                                         ccTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGc                                            >  6:10161/1‑147 (MQ=255)
                                                                                                                aGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGaa                                                              >  3:20137/1‑122 (MQ=255)
                                                                                                                aGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGaa                                                              <  4:20137/122‑1 (MQ=255)
                                                                                                                  ggTGTTCGACGAGCAGGTCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagtagca                                                        <  5:280099/126‑1 (MQ=255)
                                                                                                                  ggTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagtagca                                                        >  6:280099/1‑126 (MQ=255)
                                                                                                                   gtgtTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCgg                                 <  1:121752/148‑1 (MQ=255)
                                                                                                                   gtggTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGt                                <  6:230589/149‑1 (MQ=255)
                                                                                                                     gtTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGc                                            <  5:192912/135‑1 (MQ=255)
                                                                                                                     gtTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGc                                            >  6:192912/1‑135 (MQ=255)
                                                                                                                      ttCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGa                              <  7:94589/148‑1 (MQ=255)
                                                                                                                      ttCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGa                              <  2:53173/148‑1 (MQ=255)
                                                                                                                        cgacgaGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCt                           >  2:289460/1‑149 (MQ=255)
                                                                                                                                 gggCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCt                  >  7:134944/1‑149 (MQ=255)
                                                                                                                                      cacTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGccc             <  8:134944/149‑1 (MQ=255)
                                                                                                                                      cacTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGccc             <  2:102261/149‑1 (MQ=255)
                                                                                                                                           gCCCAATTCGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGt           >  4:211737/1‑146 (MQ=255)
                                                                                                                                                 ttGGGCCAGCAATTACAACTACTTTATGCCGGTACTCCGTAATATCAGCAATACCTCTCGGATCAACGATGAACGCCGCGCATCCAGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGATTTTGCGCCTGTCCCTTGAGAAACgg  <  6:184113/149‑1 (MQ=255)
                                                                                                                                                 ttGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAACgg  <  1:289460/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAACGG  >  NC_002947/985790‑986083

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: