Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 502,121 A→T 7.6% Y137F (TAT→TTT)  PP_0413 → polyamine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947502,1210AT 7.6% 211.2 / 4.7 75Y137F (TAT→TTT) PP_0413polyamine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base T (2/4);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.76e-01

AAGGCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCAAGG  >  NC_002947/501973‑502267
                                                                                                                                                    |                                                                                                                                                   
aaGGCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTta                                                                                                                                                     >  3:51708/1‑149 (MQ=255)
aaGGCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTta                                                                                                                                                     >  6:156472/1‑149 (MQ=255)
 aGGCGTTGCGCCAGCATTTGGCCGAGTCGCGCAAGAACCAGACCCATCGTCACCTTTCCAAACGCTTCAACATGGAGCTCGCCGGTTACCGCAGCTTGTTGACCAACACTCCCCGTGCGCAGCCCTTCAAGGCTGAGCCTGCCTCTta                                                                                                                                                     >  5:236777/1‑148 (MQ=255)
   gCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATaa                                                                                                                                                  <  2:137019/149‑1 (MQ=255)
          gCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGca                                                                                                                                           <  2:117239/149‑2 (MQ=255)
           ccAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGcag                                                                                                                                          <  4:316771/149‑3 (MQ=255)
           ccAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGt                                                                                                                                          >  5:275124/1‑149 (MQ=255)
           ccAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGt                                                                                                                                          >  1:358056/1‑149 (MQ=255)
                       cGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGtt                                                                                                                               >  6:139460/1‑148 (MQ=255)
                         aGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGAc                                                                                                                            <  5:156472/149‑1 (MQ=255)
                           tCGCGCAAGAACCAGACCCTTTGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGa                                                                                                                          <  1:309098/149‑1 (MQ=255)
                           tCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGa                                                                                                                          <  4:355068/149‑1 (MQ=255)
                           tCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAAAGCTTGAACATGGAGCTGGCCGGTTCCCGCAGCTTGTTGACCAAGCATGCCCAGGCGCTGCCCTGCAAGGAGCACCATCCCGATTATAAATCGGCGTGGCAGGCTGTTgcgct                                                                                                                          >  8:69852/1‑145 (MQ=255)
                             gcgcAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGc                                                                                                                        <  1:41759/149‑1 (MQ=255)
                                cAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGtt                                                                                                                     >  5:20820/1‑149 (MQ=255)
                                   gAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTggg                                                                                                                  <  1:147494/149‑1 (MQ=255)
                                   gAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTggg                                                                                                                  <  6:275124/149‑1 (MQ=255)
                                   gAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTggg                                                                                                                  <  3:136566/149‑1 (MQ=255)
                                   gAACCAGAACCTTGGGGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTGACCGCAGCTTGGTGACCCAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCGTAGAAAGCGGCGTTGCAGGCGGTGGACGAGcttagcg                                                                                                                  >  4:276278/1‑143 (MQ=255)
                                     aCCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGc                                                                                                                <  6:225259/149‑1 (MQ=255)
                                            ccTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCgg                                                                                                        >  1:25869/1‑149 (MQ=255)
                                            ccTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCgg                                                                                                        >  7:25996/1‑149 (MQ=255)
                                             cTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGc                                                                                                       <  4:364833/149‑1 (MQ=255)
                                                 gTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTAt                                                                                                   >  5:65359/1‑149 (MQ=255)
                                                      cTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCATCTTGTTGACCTAGACTGCCCGTGCTCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCa                                                                                              <  5:19959/149‑1 (MQ=255)
                                                      cTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTgg                                                                                                >  4:94684/1‑147 (MQ=255)
                                                      cTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTAATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGAACCGGCTTATTgg                                                                                                >  4:99984/1‑147 (MQ=255)
                                                      cTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACAGTAGCTTCTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTgg                                                                                                >  4:107192/1‑147 (MQ=255)
                                                       tttCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAg                                                                                             >  4:60208/1‑149 (MQ=255)
                                                                 cTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGca                                                                                                                                           <  6:373840/94‑2 (MQ=255)
                                                                 cTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGca                                                                                                                                           >  5:373840/1‑93 (MQ=255)
                                                                    gAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGCCCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATAcgccgc                                                                                 >  6:226163/1‑148 (MQ=255)
                                                                    gAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATAcgccgc                                                                                 >  2:157937/1‑148 (MQ=255)
                                                                    gAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATAcgccg                                                                                  <  5:139460/147‑1 (MQ=255)
                                                                                 gCCGGTTTCCGCTGCTTGTTGACCAAGTCTGCCGGTGGGGTGCCCTTGTAGGCTGAGGCTGCGTCTTTTAAAGAGGCGTTGCAGGGTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCATTACCAGTTCGGt                                                                   <  6:236777/149‑1 (MQ=255)
                                                                                 gCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGt                                                                   <  2:358056/149‑1 (MQ=255)
                                                                                 gCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGt                                                                   <  4:51708/149‑1 (MQ=255)
                                                                                 gCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGATGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCGGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGt                                                                   <  6:65359/149‑1 (MQ=255)
                                                                                 gCCGGTTACCGCAGCTTGTTCACCAAGACTCCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCGGTTCg                                                                     <  8:135206/147‑1 (MQ=255)
                                                                                      ttACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGActt                                                              >  6:272895/1‑149 (MQ=255)
                                                                                           gCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACttctttc                                                         <  2:25869/149‑1 (MQ=255)
                                                                                           gCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACttctttc                                                         <  3:164631/149‑1 (MQ=255)
                                                                                              gCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTAc                                                      <  8:25996/149‑1 (MQ=255)
                                                                                                       ccAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGcc                                             >  8:218566/1‑149 (MQ=255)
                                                                                                        cAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCt                                            >  2:53578/1‑149 (MQ=255)
                                                                                                            aCTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGt                                        <  1:53578/149‑1 (MQ=255)
                                                                                                             cTGCCCGTGCGCTGCCCGTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGt                                        <  6:20820/148‑1 (MQ=255)
                                                                                                                cccGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTgg                                      <  3:60208/147‑1 (MQ=255)
                                                                                                                    tGCGCTGCCCTTCAAGGCTGAGCCTTCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGTCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatc                                <  7:69852/149‑1 (MQ=255)
                                                                                                                    tGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatc                                <  5:112104/149‑1 (MQ=255)
                                                                                                                       gcTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatcgca                             >  7:22485/1‑149 (MQ=255)
                                                                                                                            ccTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGAt                        <  3:15896/149‑1 (MQ=255)
                                                                                                                             cTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATg                       >  1:334899/1‑149 (MQ=255)
                                                                                                                              ttCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTGTACCTGCTGGCCTCGGTGGAGCATCGCATCGATGa                      >  3:29342/1‑149 (MQ=255)
                                                                                                                              ttCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatcgcatc                           <  5:43726/144‑1 (MQ=255)
                                                                                                                              ttCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatcgcatc                           >  6:43726/1‑144 (MQ=255)
                                                                                                                              ttCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGa                      >  4:146314/1‑149 (MQ=255)
                                                                                                                              ttCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGAACCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGa                      >  6:98477/1‑149 (MQ=255)
                                                                                                                                cAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCGGCAATACCAGTTCGGTGACTTCTTTCTAGCTGCTGGCGTCGGTGGACCATCGCATCGATGa                      <  3:107192/147‑1 (MQ=255)
                                                                                                                                 aaGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAATTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGAcc                    <  1:157937/148‑1 (MQ=255)
                                                                                                                                   ggCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTGTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTg                 >  6:84395/1‑149 (MQ=255)
                                                                                                                                   ggCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTg                 <  5:272895/149‑1 (MQ=255)
                                                                                                                                    gCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGCCGCGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACCCCGCAATACCAGTTCGGTGCCTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGCCCTTgg                >  5:288224/1‑149 (MQ=255)
                                                                                                                                     cTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTgg                                      >  1:226777/1‑126 (MQ=255)
                                                                                                                                     cTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTgg                                      <  2:226777/126‑1 (MQ=255)
                                                                                                                                     cTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGt               <  5:226163/149‑1 (MQ=255)
                                                                                                                                      tGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTg                                       <  8:214752/124‑1 (MQ=255)
                                                                                                                                      tGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGt               >  1:108121/1‑148 (MQ=255)
                                                                                                                                      ctgcaacgcctCTTTATAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGt               >  4:118188/8‑148 (MQ=255)
                                                                                                                                        aGCCTGCCTCTTATACAGAGGCGTTGGAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTg              <  3:99984/147‑1 (MQ=255)
                                                                                                                                        aGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTg              <  3:94684/147‑1 (MQ=255)
                                                                                                                                        aGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTg              <  3:146314/147‑1 (MQ=255)
                                                                                                                                        aGCCTGCCTCTTATAAAGAGGCGGTGCAGGCTGTTGACCAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTTCCTGCTGGCCTCGGTGGACCATCCCAGCGATGACCTTgggg              <  6:208075/147‑3 (MQ=255)
                                                                                                                                                tctTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGcc      >  3:350480/1‑147 (MQ=255)
                                                                                                                                                 ctTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCAGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTAcat                                                    >  1:269017/1‑99 (MQ=255)
                                                                                                                                                  ttATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTgg                                               >  7:192381/1‑104 (MQ=255)
                                                                                                                                                  ttATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTgg                                               <  8:192381/104‑1 (MQ=255)
                                                                                                                                                  ttATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCAAgg  >  4:135750/1‑149 (MQ=255)
                                                                                                                                                    ataAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCAATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatc                                <  4:98697/117‑1 (MQ=255)
                                                                                                                                                    ataAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCAATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatc                                >  3:98697/1‑117 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAGGCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCAAGG  >  NC_002947/501973‑502267

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