Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 502,121 A→T 7.2% 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.2% 196.2 / 5.1 69Y137F (TAT→TTT) PP_0413polyamine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (32/32);  new base T (2/3);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

AAGGCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCAAGG  >  NC_002947/501973‑502267
                                                                                                                                                    |                                                                                                                                                  
aaGGCGTTGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTta                                                                                                                                                    >  3:335974/1‑149 (MQ=255)
 aGGCGTTGAGCTAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTtat                                                                                                                                                   <  4:335574/149‑1 (MQ=255)
      ttGAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAaga                                                                                                                                              <  7:324211/149‑1 (MQ=255)
        gAGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAgg                                                                                                                                            >  7:267292/1‑149 (MQ=255)
         aGCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGc                                                                                                                                           <  4:335974/149‑1 (MQ=255)
          gCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGca                                                                                                                                          <  7:158409/149‑2 (MQ=255)
          gCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTTATAAGAGGca                                                                                                                                          <  7:306729/149‑2 (MQ=255)
          gCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGc                                                                                                                                           <  5:219620/148‑1 (MQ=255)
          gCCAGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCg                                                                                                                                          >  3:173237/1‑149 (MQ=255)
             aGGATTTGGCCGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTg                                                                                                                                       >  4:74242/1‑149 (MQ=255)
                      ccGAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGtt                                                                                                                              >  7:317385/1‑149 (MQ=255)
                        gAGTCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGa                                                                                                                            <  7:26975/149‑1 (MQ=255)
                           tCGCGCAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGa                                                                                                                         <  4:177870/149‑1 (MQ=255)
                                cAAGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGtt                                                                                                                    >  5:175165/1‑149 (MQ=255)
                                 aaGAACCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTg                                                                                                                   <  6:4301/149‑1 (MQ=255)
                                     aCCAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGc                                                                                                               <  4:173237/149‑1 (MQ=255)
                                      ccAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCg                                                                                                              >  8:79672/1‑149 (MQ=255)
                                       cAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTgg                                                                                                                  <  7:334940/144‑1 (MQ=255)
                                       cAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTgg                                                                                                                  >  8:334940/1‑144 (MQ=255)
                                       cAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGa                                                                                                             >  4:135252/1‑149 (MQ=255)
                                       cAGACCCTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGa                                                                                                             >  1:238042/1‑149 (MQ=255)
                                            ccTTGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGCTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCgg                                                                                                        >  1:76865/1‑149 (MQ=255)
                                              ttGGTGACCTTTCCAAACGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCt                                                                                                      >  2:201370/1‑149 (MQ=255)
                                                               cGCTTGAACATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATAcg                                                                                     <  3:135252/149‑1 (MQ=255)
                                                                     aaCATGGAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCaa                                                                               <  8:317385/149‑1 (MQ=255)
                                                                           gAGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCa                                                                          <  7:79672/148‑1 (MQ=255)
                                                                            aGCTGGCCGGTTACCTCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCa                                                                          <  4:220488/147‑1 (MQ=255)
                                                                            aGCTGGCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAg                                                                         >  3:220488/1‑148 (MQ=255)
                                                                                ggCCGGTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCgg                                                                    >  1:22595/1‑149 (MQ=255)
                                                                                     gTTACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGAct                                                               >  8:139237/1‑149 (MQ=255)
                                                                                      ttACCGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGActt                                                              <  3:74242/149‑1 (MQ=255)
                                                                                         ccGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGActtct                                                            <  2:238042/148‑1 (MQ=255)
                                                                                          cGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACttcttt                                                          >  6:272401/1‑149 (MQ=255)
                                                                                          cGCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACttcttt                                                          >  5:343433/1‑149 (MQ=255)
                                                                                           gCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAg                                                                                                                                    <  7:334779/74‑1 (MQ=255)
                                                                                           gCAGCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAg                                                                                                                                    >  8:334779/1‑74 (MQ=255)
                                                                                              gCTTGTTGACCAAGACTGCCCGTGCGCTGCCCTTCACGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCTATACGCCGCAATACCAGTTCGGTGACTTCTTTCTAc                                                      <  5:188269/149‑1 (MQ=255)
                                                                                                      aCCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCa                                                                                              <  8:132441/101‑1 (MQ=255)
                                                                                                      aCCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCa                                                                                              >  7:127547/1‑101 (MQ=255)
                                                                                                      aCCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCa                                                                                              >  7:132441/1‑101 (MQ=255)
                                                                                                      aCCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCa                                                                                              <  8:127547/101‑1 (MQ=255)
                                                                                                      aCCAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGc                                              >  5:185551/1‑149 (MQ=255)
                                                                                                       ccAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGcc                                             >  5:146867/1‑149 (MQ=255)
                                                                                                        cAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATAGGCCGCAATACCa                                                                          >  7:285788/1‑119 (MQ=255)
                                                                                                        cAAGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCa                                                                          <  8:285788/119‑1 (MQ=255)
                                                                                                         aaGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTc                                           >  3:278567/1‑149 (MQ=255)
                                                                                                          aGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTCTTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCg                                          >  8:238695/1‑149 (MQ=255)
                                                                                                          aGACTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCg                                          <  6:175165/149‑1 (MQ=255)
                                                                                                            aCTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGt                                        <  8:267292/149‑1 (MQ=255)
                                                                                                            aCTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGt                                        <  5:226127/149‑1 (MQ=255)
                                                                                                             cTGCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTg                                       <  6:185551/149‑1 (MQ=255)
                                                                                                               gCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGt                                                                   >  2:259698/1‑119 (MQ=255)
                                                                                                               gCCCGTGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGt                                                                   <  1:259698/119‑1 (MQ=255)
                                                                                                                    tGCGCTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatc                                <  2:22595/149‑1 (MQ=255)
                                                                                                                      cgcTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatcgc                              <  1:201370/149‑1 (MQ=255)
                                                                                                                       gcTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCg                                          <  3:91825/136‑1 (MQ=255)
                                                                                                                       gcTGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCg                                          >  4:91825/1‑136 (MQ=255)
                                                                                                                         tGCCCTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACcatcgcatc                           >  4:332145/1‑149 (MQ=255)
                                                                                                                           cccTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGa                         <  7:139237/149‑1 (MQ=255)
                                                                                                                             cTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACctg                                                   <  6:160456/121‑1 (MQ=255)
                                                                                                                             cTTCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACctg                                                   >  5:160456/1‑121 (MQ=255)
                                                                                                                              ttCAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTCCTGGCCTCGGTGGACCATCGAATCGATGa                      >  6:217739/1‑149 (MQ=255)
                                                                                                                                cAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGa                      >  7:238607/1‑147 (MQ=255)
                                                                                                                                cAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGa                      <  8:238607/147‑1 (MQ=255)
                                                                                                                                cAAGGCTGAGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGAcc                    >  3:187613/1‑149 (MQ=255)
                                                                                                                                        aGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTg              <  4:278567/147‑1 (MQ=255)
                                                                                                                                        aGCCTGCCTCTTATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAg            <  4:187613/149‑1 (MQ=255)
                                                                                                                                            cgcctCTTTATAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGt                                                                   <  2:62701/89‑1 (MQ=255)
                                                                                                                                            cgcctCTTTATAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGt                                                                   >  1:62701/2‑90 (MQ=255)
                                                                                                                                            cgcctCTTTATAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAg        >  2:310317/2‑149 (MQ=255)
                                                                                                                                                  ttATAAAGAGGCGTTGCAGGCTGTTGACGAGCGTTGGGGCGACCCGGCTTATTGGCAGGCGATACGCCGCAATACCAGTTCGGTGACTTCTTTCTACCTGCTGGCCTCGGTGGACCATCGCATCGATGACCTTGGTGAGCTAGCCAAgg  >  8:202048/1‑149 (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: