Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 502,121 A→T 9.1% 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 9.1% 176.9 / 7.5 66Y137F (TAT→TTT) PP_0413polyamine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (23/37);  new base T (3/3);  total (26/40)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

AAGGCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCAAG  >  NC_002947/501973‑502266
                                                                                                                                                    |                                                                                                                                                 
aaGGCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTta                                                                                                                                                   >  4:192723/1‑149 (MQ=255)
aaGGCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTta                                                                                                                                                   <  3:62052/149‑1 (MQ=255)
aaGGCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTta                                                                                                                                                   <  6:71274/149‑1 (MQ=255)
 aGGCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTtat                                                                                                                                                  <  1:15542/149‑1 (MQ=255)
  ggCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTtat                                                                                                                                                  >  7:111849/1‑148 (MQ=255)
      ttGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAaga                                                                                                                                             <  1:220280/149‑1 (MQ=255)
          gCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGca                                                                                                                                         <  2:26386/149‑2 (MQ=255)
          gCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGca                                                                                                                                         <  5:150576/149‑2 (MQ=255)
           ccAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGt                                                                                                                                        >  5:357699/1‑149 (MQ=255)
           ccAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGt                                                                                                                                        >  2:384203/1‑149 (MQ=255)
            cAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGtt                                                                                                                                       <  3:192723/149‑1 (MQ=255)
             aGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAACCGCTTGCACATGGAGCTGGCCGGTTACCGCAGATTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTg                                                                                                                                      >  8:274623/1‑149 (MQ=255)
                aTTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAg                                                                                                                                   <  1:158412/149‑1 (MQ=255)
                 tttGGCCGAGTCGGTCAAGAACCAGAACCTTTGTGACATTTCCAAACGCTTGAACATGTAGCTGGCCTGTTACTGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAATTCTTAGCCTTCCTCTTATAAGGAGGCGTTTCAgg                                                                                                                                  <  2:135602/149‑1 (MQ=255)
                 tttGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGcag                                                                                                                                        >  1:231945/1‑141 (MQ=255)
                 tttGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGcag                                                                                                                                        <  2:231945/143‑3 (MQ=255)
                                  aGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTgg                                                                                                                 <  8:209694/149‑1 (MQ=255)
                                   gAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTggg                                                                                                                >  5:388941/1‑149 (MQ=255)
                                   gAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTggg                                                                                                                <  2:215930/149‑1 (MQ=255)
                                   gAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTggg                                                                                                                <  4:38757/149‑1 (MQ=255)
                                    aaCCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTgggg                                                                                                               <  6:388941/149‑1 (MQ=255)
                                      ccAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGc                                                                                                              >  2:176716/1‑148 (MQ=255)
                                       cAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGa                                                                                                            >  3:289676/1‑149 (MQ=255)
                                         gACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGAcc                                                                                                          <  6:91490/149‑1 (MQ=255)
                                             cTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGc                                                                                                      <  5:77487/149‑1 (MQ=255)
                                                 gTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTAt                                                                                                  <  3:181387/149‑1 (MQ=255)
                                                  tGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTAtt                                                                                                 <  8:222316/149‑1 (MQ=255)
                                                   gACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTAtt                                                                                                 <  6:357699/148‑1 (MQ=255)
                                                     ccTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGc                                                                                              <  7:22782/149‑1 (MQ=255)
                                                      cTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTgg                                                                                               >  8:69700/1‑147 (MQ=255)
                                                      cTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCACTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTgg                                                                                               >  5:28846/1‑147 (MQ=255)
                                                          ccAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGc                                                                                          >  7:322727/1‑148 (MQ=255)
                                                          ccAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGc                                                                                          <  8:111849/148‑1 (MQ=255)
                                                          ccAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGc                                                                                          <  3:35547/148‑1 (MQ=255)
                                                           cAAACGCTTGTACTTGGAGCTGGCAGGTTACCGCAGCTTGGTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCTTTGGGGCGACCCGGCTTATTTGCAGGCGa                                                                                        <  2:130766/149‑1 (MQ=255)
                                                           cAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGa                                                                                        <  7:243183/149‑1 (MQ=255)
                                                                    gAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATAcgccgc                                                                                >  4:225440/1‑148 (MQ=255)
                                                                    gAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATAcgccg                                                                                 <  8:115959/147‑1 (MQ=255)
                                                                    gAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATAcgccg                                                                                 <  6:140470/147‑1 (MQ=255)
                                                                      acctGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAACGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAAt                                                                             <  6:150934/146‑1 (MQ=255)
                                                                      aCATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCaa                                                                              >  8:165601/1‑148 (MQ=255)
                                                                       cATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATa                                                                            >  5:227575/1‑149 (MQ=255)
                                                                        aTGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATAc                                                                           <  7:69700/149‑1 (MQ=255)
                                                                             gCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGtt                                                                      <  3:225440/149‑1 (MQ=255)
                                                                                 gCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGt                                                                  <  8:181741/149‑1 (MQ=255)
                                                                                   cGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGa                                                                <  6:227575/149‑1 (MQ=255)
                                                                                            cAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCTATACGCCGCAATACCAGTTCGGTGACttctttct                                                       <  7:274623/149‑1 (MQ=255)
                                                                                              gcgtgttgACGAAGACTGCCCGCTCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAACGTCGCGTTGCAGGCTGTTGACGAGCGGTGGGGCGACCCGGCTTATTGGCTGGCGATACGCCGCAATACCAGTGCGGTGACTTCTTTCAAc                                                     <  8:303719/146‑1 (MQ=255)
                                                                                              gCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTAc                                                     <  1:384203/149‑1 (MQ=255)
                                                                                                     gACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTggg                                                                                                                <  2:133791/83‑1 (MQ=255)
                                                                                                     gACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTggg                                                                                                                >  1:133791/1‑83 (MQ=255)
                                                                                                      aCCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGc                                             >  3:228442/1‑149 (MQ=255)
                                                                                                          aGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCg                                         >  4:382252/1‑149 (MQ=255)
                                                                                                           gACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACctgctg                                               >  5:65481/1‑142 (MQ=255)
                                                                                                           gACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACctgctg                                               <  6:65481/142‑1 (MQ=255)
                                                                                                           gACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCgg                                        >  6:44711/1‑149 (MQ=255)
                                                                                                             cTGCCCGTGCGCTGCCCTTCAAGGCTGATCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTg                                      <  4:289676/149‑1 (MQ=255)
                                                                                                             cTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTg                                      <  6:28846/149‑1 (MQ=255)
                                                                                                              tGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTgg                                     >  7:118632/1‑149 (MQ=255)
                                                                                                                cccGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTgg                                     <  4:228442/147‑1 (MQ=255)
                                                                                                                cccGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGAc                                   >  2:259825/1‑149 (MQ=255)
                                                                                                                      cgcTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatcg                              <  3:382252/148‑1 (MQ=255)
                                                                                                                            ccTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGAt                       <  1:176716/149‑1 (MQ=255)
                                                                                                                              ttCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGCCCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTgg                                     >  8:240038/1‑133 (MQ=255)
                                                                                                                              ttCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTgg                                     <  7:240038/133‑1 (MQ=255)
                                                                                                                                       gAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGt              >  5:8903/1‑147 (MQ=255)
                                                                                                                                           acgcctCTTTATAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGc          >  6:313338/3‑147 (MQ=255)
                                                                                                                                            tgcctCTTATAAAGATGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGa            <  2:71500/144‑1 (MQ=255)
                                                                                                                                            tgcctCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGa            >  1:71500/1‑144 (MQ=255)
                                                                                                                                            cgcctCTTTATAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAg       >  3:25984/2‑149 (MQ=255)
                                                                                                                                               ctctTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGcc     >  5:313547/1‑148 (MQ=255)
                                                                                                                                                    ataAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCAAg  <  8:322727/146‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AAGGCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCAAG  >  NC_002947/501973‑502266

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: