Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 796,817 T→A 6.2% intergenic (+19/‑171) PP_0681 → / → PP_0682 hypothetical protein/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947796,8170TA 6.2% 188.4 / 2.7 65intergenic (+19/‑171)PP_0681/PP_0682hypothetical protein/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base A (2/2);  total (34/31)
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.06e-01

ACCATGCATGGTTATCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTA  >  NC_002947/796670‑796963
                                                                                                                                                   |                                                                                                                                                   
aCCATGCATGGTTATGAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGTTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCttt                                                                                                                                                     <  3:343916/148‑1 (MQ=255)
 ccATGCATGGTTATCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGc                                                                                                                                                   <  4:155857/149‑1 (MQ=255)
   atgcatgGTTATCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCgg                                                                                                                                                 >  8:345922/1‑149 (MQ=255)
    tgcatgCTTATCAAGCGGCCGATGACCCGCAGATCTAGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGATCGCTCTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGc                                                                                                                                                <  8:170389/149‑1 (MQ=255)
        tgGTTATCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCAAAGCGGccccg                                                                                                                                            <  3:339811/149‑3 (MQ=255)
           ttATCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATcgc                                                                                                                                         >  8:284756/1‑149 (MQ=255)
            tATCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATcgc                                                                                                                                         >  2:325112/1‑148 (MQ=255)
             aTCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGa                                                                                                                                       >  1:196995/1‑149 (MQ=255)
              tCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGa                                                                                                                                       >  7:194879/1‑148 (MQ=255)
                        gATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGc                                                                                                                           >  5:345754/1‑149 (MQ=255)
                        gATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGc                                                                                                                           >  2:269210/1‑149 (MQ=255)
                            aCCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTcc                                                                                                                        <  5:47651/148‑1 (MQ=255)
                            aCCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAAAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTcc                                                                                                                        <  1:269210/148‑1 (MQ=255)
                             cccGCAGAGCTGGTGTCTGCATTTGCAGGCGCAGTTGGGGCGGAGCATTGGTGAGCGGTGTCTCAACGTGGGGGATTACCAGTTTCATCCTTAGACCTGAATAATGCCGGGCCGGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCca                                                                                                                      <  6:41617/149‑1 (MQ=38)
                                  aGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGt                                                                                                                 <  7:110087/149‑1 (MQ=255)
                                    aTCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTa                                                                                                                >  6:344399/1‑148 (MQ=255)
                                         gtgtCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACa                                                                                                              <  8:194879/145‑1 (MQ=255)
                                          tgtCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGGCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATc                                                                                                         <  3:113638/149‑1 (MQ=255)
                                            tCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCgg                                                                                                       >  4:107068/1‑149 (MQ=255)
                                                cATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGc                                                                                                      <  2:196995/146‑1 (MQ=255)
                                                  tttGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCt                                                                                                 <  7:345922/149‑1 (MQ=255)
                                                      cAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTgg                                                                                             <  2:34884/149‑1 (MQ=255)
                                                         gcgcAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATc                                                                                                                                           <  7:346758/101‑1 (MQ=255)
                                                         gcgcAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATc                                                                                                                                           >  8:346758/1‑101 (MQ=255)
                                                         gcgcAGTTGCGGCGGATCAATGGTGAGCGGTGACGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAg                                                                                          >  8:108819/1‑149 (MQ=255)
                                                           gcAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGa                                                                                           <  2:116632/146‑1 (MQ=255)
                                                           gcAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGa                                                                                           >  1:116632/1‑146 (MQ=255)
                                                           gcAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAgtg                                                                                        <  7:284756/149‑1 (MQ=255)
                                                           gcAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAgtg                                                                                        >  5:272219/1‑149 (MQ=255)
                                                           gcAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGCGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGAAACAACGCCGCACCCACAAATACACCAAACGCTTGCTGTGGGAgtg                                                                                        >  1:176553/1‑149 (MQ=255)
                                                            cagaggcggcgGATCACTGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGCGACCTGACGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAgtgt                                                                                       <  6:283568/144‑1 (MQ=255)
                                                               ttGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGc                                                                                    >  5:365216/1‑149 (MQ=255)
                                                                  cggcggATCAATGGGGCTCGGTGGCGCAACCTGGGGGATTACCAGTTTCATCCAGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCtt                                                                                 <  3:273616/149‑1 (MQ=38)
                                                                    gcggATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTtgt                                                                               <  6:235858/149‑1 (MQ=255)
                                                                     cggATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTtgtg                                                                              <  2:310628/149‑1 (MQ=255)
                                                                       gATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTc                                                                            >  8:80529/1‑149 (MQ=255)
                                                                         tCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTcgc                                                                          <  5:344399/149‑1 (MQ=255)
                                                                         tCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTcgc                                                                          >  6:326780/1‑149 (MQ=255)
                                                                               gTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCa                                                                     <  7:80529/148‑1 (MQ=255)
                                                                               gTGAGCGGTGGCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGACGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCa                                                                     <  3:290727/148‑1 (MQ=255)
                                                                                           gCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGagcag                                                        >  5:303455/1‑149 (MQ=255)
                                                                                             aaCCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGcc                                                      >  1:47714/1‑149 (MQ=255)
                                                                                             aaCCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCCCGACACAAGGCCGCTACCACAGGTACAGCATCCGCGTGCTGTGGGAGTGTTGCCTTTTGGAGCGATAAGGAAGCGGAGCAGcc                                                      >  8:309340/1‑149 (MQ=37)
                                                                                                        ttACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCtttt                                           <  6:272219/149‑1 (MQ=255)
                                                                                                          aCCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCGTGTGc                                         >  8:197719/1‑149 (MQ=255)
                                                                                                           ccAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCt                                        >  7:271772/1‑149 (MQ=255)
                                                                                                            cAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGGGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTg                                       >  5:286107/1‑149 (MQ=255)
                                                                                                            cAGTTGCATCCTGAGACCTGACCAATGCCGGGGCCGCTTTGGGGCCCAGCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTg                                       <  3:208499/149‑1 (MQ=39)
                                                                                                               tttCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTg                                       >  3:328863/1‑146 (MQ=255)
                                                                                                               tttCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTg                                       <  4:328863/146‑1 (MQ=255)
                                                                                                                   aTCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGGGTCGCGCTCAGGCTGCGGAGCCTCCCCCCGGCTTTTGCAGTAACCgg                                >  8:241735/1‑149 (MQ=38)
                                                                                                                     ccTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGtg                              >  2:4956/1‑149 (MQ=255)
                                                                                                                      cTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGtgg                             <  1:325112/149‑1 (MQ=255)
                                                                                                                        gagaCCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGAGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCACCCGGGCTTTTGCGGTAACAGGtggct                           >  4:271763/1‑149 (MQ=255)
                                                                                                                         agaCCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGtggctg                          <  6:345754/149‑1 (MQ=255)
                                                                                                                               gaagaaTGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCtata                    >  5:128734/1‑149 (MQ=255)
                                                                                                                                aagaaTGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCtata                    <  6:365216/148‑1 (MQ=255)
                                                                                                                                aagaaTGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATAc                   <  8:271772/149‑1 (MQ=255)
                                                                                                                                   aaTGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCg                >  7:244627/1‑149 (MQ=255)
                                                                                                                                      gCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCg             >  3:246348/1‑149 (MQ=255)
                                                                                                                                      cgatgggCCGCAAAGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCg             >  4:304944/5‑149 (MQ=255)
                                                                                                                                       gatgggCCGCAAAGCGGCCCATCGCGACAGCAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCGGTGGGAGTGTGGcc                                                                                   >  1:29421/4‑79 (MQ=25)
                                                                                                                                       gatgggCCGCAAAGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGcc                                                                                   <  2:29421/75‑1 (MQ=25)
                                                                                                                                         gggggtgCTTTTCGGCCCATCGCGACACAAGGCCGGTCCCACAGGTAGAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGTTCAGGCTGCGTAGCAGGCCCCGGGCTTTTGCTGTAACAGGtgg                             <  7:250881/124‑1 (MQ=38)
                                                                                                                                         ggggCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGCTACAGCATCGGCTTGCTGTGGGAGTGTGTACTTGTGTCGCGAACAGGCTGCGGAGCAGCCCCCGGGCTATTGCTGTAAAATGtgg                             >  8:250881/1‑130 (MQ=37)
                                                                                                                                             ccGCTTTGCGGCCCATCGCGGCACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGAAACAGGTGGCTGGCTATACTCGCCGTCGCttt      >  5:355897/1‑149 (MQ=255)
                                                                                                                                                   tGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTa  <  8:244627/147‑1 (MQ=255)
                                                                                                                                                   tGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTa  <  2:47714/147‑1 (MQ=255)
                                                                                                                                                   tGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTa  <  5:326780/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCATGCATGGTTATCAAGTGGCCGATGACCCGCAGATCTGGTGTCTGCATTTGCAGGCGCAGTTGCGGCGGATCAATGGTGAGCGGTGTCGCAACCTGGGGGATTACCAGTTTCATCCTGAGACCTGAAGAATGCCGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGGTACAGCATCGGCTTGCTGTGGGAGTGTGGCCTTGTGTCGCGATCAGGCTGCGGAGCAGCCCCCGGGCTTTTGCTGTAACAGGTGGCTGGCTATACTCGCCGTCGCTTTTTTA  >  NC_002947/796670‑796963

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