Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 502,123 A→T 9.1% K138* (AAA→TAA)  PP_0413 → polyamine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947502,1230AT 9.1% 165.8 / 7.5 66K138* (AAA→TAA) 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 = 4.60e-01

TTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCAAG  >  NC_002947/501979‑502266
                                                                                                                                                |                                                                                                                                               
ttGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAaga                                                                                                                                             <  1:220280/149‑1 (MQ=255)
    gCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGca                                                                                                                                         <  5:150576/149‑2 (MQ=255)
    gCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGca                                                                                                                                         <  2:26386/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                                                                                                                                        <  2:231945/143‑3 (MQ=255)
           tttGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGcag                                                                                                                                        >  1:231945/1‑141 (MQ=255)
                            aGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTgg                                                                                                                 <  8:209694/149‑1 (MQ=255)
                             gAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTggg                                                                                                                <  2:215930/149‑1 (MQ=255)
                             gAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTggg                                                                                                                >  5:388941/1‑149 (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                                                                                          <  3:35547/148‑1 (MQ=255)
                                                    ccAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGc                                                                                          <  8:111849/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                                                                                                                >  1:133791/1‑83 (MQ=255)
                                                                                               gACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTggg                                                                                                                <  2:133791/83‑1 (MQ=255)
                                                                                                aCCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGc                                             >  3:228442/1‑149 (MQ=255)
                                                                                                    aGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCg                                         >  4:382252/1‑149 (MQ=255)
                                                                                                     gACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACctgctg                                               <  6:65481/142‑1 (MQ=255)
                                                                                                     gACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACctgctg                                               >  5:65481/1‑142 (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)
                                                                                                                                                |                                                                                                                                               
TTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCAAG  >  NC_002947/501979‑502266

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