Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 465,446 A→G 6.2% D122D (GAT→GAC) ‡ PP_0384 ← AsnC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947465,4460AG 6.2% 187.3 / 2.2 64D122D (GAT→GAC) ‡PP_0384AsnC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base G (2/2);  total (30/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 = 9.28e-01

ACCCGCTCCTACAACCGGCCTCATCGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCTTGCCGGAGCGTGGCTGCACCTTGATCCCGACATACG  >  NC_002947/465298‑465594
                                                                                                                                                    |                                                                                                                                                    
aCCCGCGCCTACAACCCGCCCCATCCAAAACCCCACCCGAGCCCCAGCCCGTTCAAACCGGCTTTCCACGGTCCACCTTGCGTCTGACAATGATCCACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGa                                                                                                                                                      <  3:301946/149‑1 (MQ=255)
     cTCCTACACCCGGCCTCATCGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAg                                                                                                                                                 <  8:288714/149‑1 (MQ=255)
        cTACAACCGGCCTCATCGATAACCCCACCCGCGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAgctg                                                                                                                                              >  5:96762/1‑149 (MQ=255)
        cTACAACCGGCCTCATCGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAgct                                                                                                                                               <  5:141079/148‑1 (MQ=255)
         tACAACCGGCCTCATCGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAgctgc                                                                                                                                             >  6:23525/1‑149 (MQ=255)
             aCCGGCCTCATCGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGATGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAGCTGGATCAGTTGCTCga                                                                                                                                         <  2:267108/149‑2 (MQ=255)
             aCCGGCCTCATCGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGATGTCAAGGTTCCACCTTGCTGCTGATAATGATCGACGTGGTGGTTTTCTCCACACCGTAGACACTGCCAATCTGGTCGAGCAGCTGGATCAGTTGCTCga                                                                                                                                         <  6:343918/149‑2 (MQ=255)
               cGGCCTCATCGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGc                                                                                                                                       <  8:265908/149‑1 (MQ=255)
               cGGCCTCATCGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGc                                                                                                                                       <  3:76159/149‑1 (MQ=255)
                 gCCTCATCGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGa                                                                                                                                     <  7:217634/149‑1 (MQ=255)
                  ccTCATCGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAg                                                                                                                                    <  7:67583/149‑1 (MQ=255)
                      aTCGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCg                                                                                                                                <  8:314166/149‑1 (MQ=255)
                      aTCGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGACCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTc                                                                                                                                  <  2:127836/147‑1 (MQ=255)
                        cGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTa                                                                                                                              >  5:272429/1‑149 (MQ=255)
                         gATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTAc                                                                                                                             >  8:153743/1‑149 (MQ=255)
                                        gcgcCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACAt                                                                                                              <  6:353598/149‑1 (MQ=255)
                                        gcgcCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACAt                                                                                                              <  6:349527/149‑1 (MQ=255)
                                            cAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGc                                                                                                                   <  6:233580/140‑1 (MQ=255)
                                            cAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGc                                                                                                                   >  5:233580/1‑140 (MQ=255)
                                                ccGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAAt                                                                                                      >  8:275934/1‑149 (MQ=255)
                                                   ttCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATcc                                                                                                                                                   <  6:359295/101‑1 (MQ=255)
                                                   ttCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATcc                                                                                                                                                   >  5:359295/1‑101 (MQ=255)
                                                    tCAAAGCGGCTGTCCACGGTCCTCCTTGCTGCAGCGAATGATCGACGTGTTGGTTTTCTCCACACCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTcgc                                                                                                  <  8:262206/149‑1 (MQ=255)
                                                         gCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCt                                                                                              >  7:333221/1‑148 (MQ=255)
                                                                  cACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGCCACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGcaca                                                                                     >  4:331909/1‑148 (MQ=255)
                                                                   aCGGTCCACCTTGCTGATGAGAATGATCGACGTGGTGGTTTTCTCCACCCAGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGCGCGAGTACGTACGCAGCCAGCCGACATAGTCGAATTAGCCGCTTACCGcacaca                                                                                   >  6:53192/1‑149 (MQ=255)
                                                                      gTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAgct                                                                                <  2:46625/149‑1 (MQ=255)
                                                                       tCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAgctg                                                                               <  7:153743/149‑1 (MQ=255)
                                                                        ccACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAgctg                                                                               <  8:333221/148‑1 (MQ=255)
                                                                         cACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAgctgct                                                                             >  7:155727/1‑149 (MQ=255)
                                                                              tgctgcTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCAcc                                                                        >  7:343840/1‑149 (MQ=255)
                                                                                ctgctgAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTg                                                                      <  1:32939/149‑1 (MQ=255)
                                                                                 tgctgAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTgg                                                                     <  8:239967/149‑1 (MQ=255)
                                                                                 tgctgAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTgg                                                                     >  4:37543/1‑149 (MQ=255)
                                                                                      agaATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCAt                                                                <  8:155727/149‑1 (MQ=255)
                                                                                       gaATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATg                                                               <  8:22246/149‑1 (MQ=255)
                                                                                        aaTGATCGTCGTGGTGTTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATgg                                                              <  6:380107/149‑1 (MQ=255)
                                                                                        aaTGATCGACGTGGTGGTTTTCTCCCCCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATgg                                                              >  3:14478/1‑149 (MQ=255)
                                                                                        aaTGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATgg                                                              >  7:168602/1‑149 (MQ=255)
                                                                                          ttaTCGACGTGGTGGTTTTCTCCACCCCGTGGATACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCTAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGc                                                             <  5:53192/146‑1 (MQ=255)
                                                                                          tGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCa                                                            <  6:96762/149‑1 (MQ=255)
                                                                                           gATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCAc                                                           >  4:86636/1‑149 (MQ=255)
                                                                                            aTCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACt                                                          <  5:23525/149‑1 (MQ=255)
                                                                                                       ggtTTTCTCCACCCCGTCGACACTGCCAATGTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGTAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGca                                               >  1:31791/1‑149 (MQ=255)
                                                                                                       ggtTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGcgccgcc                                               >  8:177627/1‑148 (MQ=255)
                                                                                                       ggtTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGca                                               >  2:163664/1‑149 (MQ=255)
                                                                                                       ggtTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGca                                               >  2:342337/1‑149 (MQ=255)
                                                                                                          tttCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCCCACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGcacca                                            >  2:75153/1‑149 (MQ=255)
                                                                                                          tttCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGcacca                                            >  3:217135/1‑149 (MQ=255)
                                                                                                          gaaaaccACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCAcc                                                                        >  8:239363/6‑121 (MQ=255)
                                                                                                          gaaaaccACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCAcc                                                                        <  7:239363/116‑1 (MQ=255)
                                                                                                             ctcCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACg                                          <  6:272429/148‑1 (MQ=255)
                                                                                                              tcCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTc                                        >  2:278920/1‑149 (MQ=255)
                                                                                                                             cTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCTTGCCGGAGCGTgg                         <  4:14478/149‑1 (MQ=255)
                                                                                                                             cTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCTTGCCGGAGCGTgg                         <  4:217135/149‑1 (MQ=255)
                                                                                                                             cTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCTTGCCGGAGCGTgg                         <  1:195375/149‑1 (MQ=255)
                                                                                                                              tGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCTTGCCGGAGCGTGGc                        >  2:100240/1‑149 (MQ=255)
                                                                                                                                ccAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCTTGCCGGAGCGTGGCTg                      >  6:57398/1‑149 (MQ=255)
                                                                                                                                        gTCGATCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCTTGCCGGAGCGTGGCTGCACCt                 <  8:343840/146‑1 (MQ=255)
                                                                                                                                         ccGAGCAGCTGGATCAGTTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCGATGGCACTCAGGCGCCGCACCACGTCCTTGCCGGAGCGTGGCTGCACCTTGAt             >  1:308960/2‑149 (MQ=255)
                                                                                                                                          cGAGCAGCTGGATCAGTTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCTTGCCGGAGCGTGGCTGCACCTTGATc            >  3:53655/1‑149 (MQ=255)
                                                                                                                                             gCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCTTGCCGGAGCGTGGCTGCACCTTGATCCCg         <  1:100240/149‑1 (MQ=255)
                                                                                                                                               aaCTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCtt                                     <  1:110916/119‑1 (MQ=255)
                                                                                                                                               aaCTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCtt                                     >  2:110916/1‑119 (MQ=255)
                                                                                                                                               aaCTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCTTGCCGGAGCGTGGCTGCACCTTGATCCCGa        >  8:49277/1‑148 (MQ=255)
                                                                                                                                                    aTCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCTTGCCGGAGCGTGGCTGCACCTTGATCCCGACATACg  >  2:235825/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCCGCTCCTACAACCGGCCTCATCGATAACCCCACCCGAGCGCCAGCCCGTTCAAAGCGGCTGTCCACGGTCCACCTTGCTGCTGAGAATGATCGACGTGGTGGTTTTCTCCACCCCGTCGACACTGCCAATCTGGTCGAGCAACTGATCCAGCTGCTCGGGCGAGTCCGTACGCAGCCAGGCGACATAGTCGAATTCGCCGCTTACCGCACACAGCTGCTGCACCTGGCCCATGGCACTCAGGCGCCGCACCACGTCCTTGCCGGAGCGTGGCTGCACCTTGATCCCGACATACG  >  NC_002947/465298‑465594

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: ‑