Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 985,932 G→C 5.9% E232Q (GAA→CAA)  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,9320GC 5.9% 216.3 / 3.9 68E232Q (GAA→CAA) pilFtype IV pili biogenesis protein PilF
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base C (2/2);  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 = 7.95e-01

TGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGT  >  NC_002947/985790‑986074
                                                                                                                                              |                                                                                                                                              
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                                                                                                                                   >  7:71337/1‑110 (MQ=255)
                                            actacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGTTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCCAATTCGGCCAgtg                                                                                                                                   <  8:71337/112‑3 (MQ=255)
                                             ctacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGTAGGGCACACTGGCCGAATTGGGCCAGCAAt                                                                                                                                 >  1:342449/1‑113 (MQ=255)
                                             ctacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGTAGGGCACACTGGCCGAATTGGGCCAGCAAt                                                                                                                                 <  2:342449/113‑1 (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                                                                                                                                 >  8:155310/1‑109 (MQ=255)
                                                 gATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAAt                                                                                                                                 <  7:155310/109‑1 (MQ=255)
                                                 gATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGAc                                                                                                                        >  3:276893/1‑118 (MQ=255)
                                                 gATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGAc                                                                                                                        <  4:276893/118‑1 (MQ=255)
                                                 gATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATa                                                                                         <  3:230323/149‑1 (MQ=255)
                                                  aTCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATAc                                                                                        <  1:99324/149‑1 (MQ=255)
                                                  aTCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATAc                                                                                        <  3:9944/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                                                                                        <  7:114193/127‑1 (MQ=255)
                                                                        ccACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATAc                                                                                        >  8:114193/1‑127 (MQ=255)
                                                                        ccACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGaa                                                                     >  4:139141/1‑146 (MQ=255)
                                                                         cACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAAcgcc                                                                 >  8:94589/1‑149 (MQ=255)
                                                                           cgacgCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGTTGAAcgccgc                                                               >  1:156430/1‑149 (MQ=255)
                                                                               gCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACg                                                                          >  7:195306/1‑134 (MQ=255)
                                                                               gCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACg                                                                          <  8:195306/134‑1 (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                                   >  6:192912/1‑135 (MQ=255)
                                                                                                                     gtTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGc                                   <  5:192912/135‑1 (MQ=255)
                                                                                                                      ttCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGa                     <  2:53173/148‑1 (MQ=255)
                                                                                                                      ttCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGa                     <  7:94589/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)
                                                                                                                                              |                                                                                                                                              
TGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCCGAATTGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGT  >  NC_002947/985790‑986074

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: