Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 74,159 A→T 9.3% K212M (AAG→ATG)  PP_0063 → lipid A biosynthesis lauroyl acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294774,1590AT 9.3% 276.7 / 14.8 97K212M (AAG→ATG) PP_0063lipid A biosynthesis lauroyl acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (44/44);  new base T (3/6);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

GCGTGGCACCGTCGACCAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGC  >  NC_002947/74012‑74306
                                                                                                                                                   |                                                                                                                                                    
gcgTGGCACCGTCGACCAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCaa                                                                                                                                                      <  4:184718/148‑1 (MQ=255)
 cgTGGCACCGTCGACCAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGt                                                                                                                                                    <  4:71807/149‑1 (MQ=255)
       aCCGTCGACCAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGc                                                                                                                                              >  5:4569/1‑149 (MQ=255)
          gTCGACCAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGa                                                                                                                                           <  8:294869/149‑1 (MQ=255)
          gTCGACCAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGa                                                                                                                                           <  2:195721/149‑1 (MQ=255)
               ccAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATg                                                                                                                                      >  5:292889/1‑149 (MQ=255)
               ccAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGAGCCGAACATg                                                                                                                                      >  4:129740/1‑149 (MQ=255)
                cAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCCGGCGCTGACCAGCAAGTTCGTGCCGAACATGc                                                                                                                                     >  7:264355/1‑149 (MQ=255)
                cAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGc                                                                                                                                     >  8:273654/1‑149 (MQ=255)
                 aaggaaggCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCt                                                                                                                                    <  3:278689/149‑1 (MQ=255)
                  aggaaggCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCATGTTCGGCACGAACTTGCTg                                                                                                                                   <  8:27270/149‑1 (MQ=255)
                  aggaaggCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTg                                                                                                                                   <  3:196621/149‑1 (MQ=255)
                   ggaaggCATCCTCAGCGTGATCAAGGAAGTGTGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCATGTTCGGCACGAACTTGCTgg                                                                                                                                  <  7:31884/149‑1 (MQ=255)
                   ggaaggCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTgg                                                                                                                                  <  6:169163/149‑1 (MQ=255)
                   ggaaggCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTgg                                                                                                                                  >  6:125224/1‑149 (MQ=255)
                    gaaggCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGc                                                                                                                                 <  7:171175/149‑1 (MQ=255)
                      aggCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCg                                                                                                                               <  5:296196/149‑1 (MQ=255)
                       ggCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCgg                                                                                                                              >  7:143056/1‑149 (MQ=255)
                          aTCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTgg                                                                                                                           >  4:107464/1‑149 (MQ=255)
                          aTCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTgg                                                                                                                           >  1:158711/1‑149 (MQ=255)
                             cTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGc                                                                                                                          >  5:324277/1‑147 (MQ=255)
                             cTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCaa                                                                                                                        <  6:324277/149‑1 (MQ=255)
                              tCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGACGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAg                                                                                                                       >  6:214968/1‑149 (MQ=255)
                              tCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTACTGCCGACCAGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAg                                                                                                                       >  5:205771/1‑149 (MQ=255)
                                 gCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCg                                                                                                                   <  8:100245/149‑1 (MQ=255)
                                 gCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCg                                                                                                                   <  1:14087/149‑1 (MQ=255)
                                   gTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGt                                                                                                                 <  2:158711/149‑1 (MQ=255)
                                    tGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCgg                                                                                                                  >  8:305450/1‑147 (MQ=255)
                                       tCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgt                                                                                                             <  6:292889/149‑1 (MQ=255)
                                       tCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGCACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgt                                                                                                             <  3:94048/149‑1 (MQ=255)
                                        cAAGGAAGTGCGGCGTGGTGGCCAGGTGGGTTTGCCTGCCGACGCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCGAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgt                                                                                                             <  8:264355/148‑1 (MQ=255)
                                         aaGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgtgt                                                                                                           >  7:275478/1‑149 (MQ=255)
                                         aaGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCAGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgtgt                                                                                                           <  6:205771/149‑1 (MQ=255)
                                           ggAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCt                                                                                                                                    <  8:188618/123‑1 (MQ=255)
                                           ggAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCt                                                                                                                                    >  7:188618/1‑123 (MQ=255)
                                            gAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgt                                                                                                             >  8:144508/1‑144 (MQ=255)
                                            gAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgt                                                                                                             <  7:144508/144‑1 (MQ=255)
                                               gTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCt                                                                                                     <  6:4569/149‑1 (MQ=255)
                                               gTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCt                                                                                                     <  8:143056/149‑1 (MQ=255)
                                                  cgcCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTg                                                                                                                                   <  5:9649/117‑1 (MQ=255)
                                                  cgcCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTg                                                                                                                                   >  6:9649/1‑117 (MQ=255)
                                                  cgcCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGc                                                                                                   <  8:275478/148‑1 (MQ=255)
                                                         gtggCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTg                                                                                           <  8:75249/149‑1 (MQ=255)
                                                           ggCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCg                                                                                         <  7:305450/149‑1 (MQ=255)
                                                            gCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGATCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCgg                                                                                        >  2:77368/1‑149 (MQ=255)
                                                                  tGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATCCGCTGCGGTTGCCg                                                                                  <  7:273654/149‑1 (MQ=255)
                                                                      tATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATg                                                                              >  6:232431/1‑149 (MQ=255)
                                                                        ttCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGc                                                                            >  7:258641/1‑149 (MQ=255)
                                                                         tCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCt                                                                           >  3:115368/1‑149 (MQ=255)
                                                                             gCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGGGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCggg                                                                       <  7:95099/149‑1 (MQ=255)
                                                                                gACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCtt                                                                    >  2:21454/1‑149 (MQ=255)
                                                                                 aCCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCtt                                                                    >  2:18291/1‑148 (MQ=255)
                                                                                 aCCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCAGGCTTc                                                                   <  6:111158/149‑1 (MQ=255)
                                                                                   ccGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGc                                                                      >  2:323495/1‑144 (MQ=255)
                                                                                   ccGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGc                                                                      <  1:323495/144‑1 (MQ=255)
                                                                                   ccGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCaa                                                                 >  4:23343/1‑149 (MQ=255)
                                                                                    cGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAg                                                                <  5:214968/149‑1 (MQ=255)
                                                                                      gAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGgt                                                              >  1:291260/1‑149 (MQ=255)
                                                                                        gCCGGCGGAATCGGCGGGCGTGTTCGTACCGGTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGATCATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGGTGCCGGATGGCTCGGGCTTCAAGgtgt                                                            <  7:54209/149‑1 (MQ=255)
                                                                                         ccGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGtt                                                           >  1:13960/1‑149 (MQ=255)
                                                                                           ggcggAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTcc                                                         >  5:129605/1‑149 (MQ=255)
                                                                                           ggcggAATCGGCGGGCGTGTTCGCACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTcc                                                         >  3:195064/1‑149 (MQ=255)
                                                                                            gcggAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCt                                                        <  4:115368/149‑1 (MQ=255)
                                                                                             cggAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTc                                                       >  8:181468/1‑149 (MQ=255)
                                                                                                aaTCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGaa                                                    <  5:29282/149‑1 (MQ=255)
                                                                                                  tCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGc                                                  >  1:323702/1‑149 (MQ=255)
                                                                                                    ggcgggcgTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCt                                                 <  1:21454/148‑1 (MQ=255)
                                                                                                     gcgggcgTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTgc                                               >  6:161345/1‑149 (MQ=255)
                                                                                                                 gTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAg                                                                >  6:26316/1‑120 (MQ=255)
                                                                                                                   aCCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGTGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATg                                  >  4:50651/1‑148 (MQ=255)
                                                                                                                   aCCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGAGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATg                                  >  6:147033/1‑148 (MQ=255)
                                                                                                                    ccGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGGTAGGTGTGTTCATGCATGCGGGCCGTTTGCCGGATGGCTCGGGGTTCAAGGTGTTCGTCGAAGCTGCGCGGAaggggatt                                 >  1:96849/1‑142 (MQ=255)
                                                                                                                     cGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTa                                >  1:249686/1‑148 (MQ=255)
                                                                                                                     cGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTAc                               <  3:107464/149‑1 (MQ=255)
                                                                                                                       ttCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAg                             >  2:44713/1‑149 (MQ=255)
                                                                                                                           tCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCAcc                         >  2:108025/1‑149 (MQ=255)
                                                                                                                              gCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAg                             >  7:1323/1‑142 (MQ=255)
                                                                                                                              gCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAg                             <  8:1323/142‑1 (MQ=255)
                                                                                                                                aCCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGAtgt                    <  1:18291/149‑1 (MQ=255)
                                                                                                                                aCCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGAtgt                    <  5:125224/149‑1 (MQ=255)
                                                                                                                                  ccAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGc                            >  1:120620/1‑139 (MQ=255)
                                                                                                                                  ccAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGc                            <  2:120620/139‑1 (MQ=255)
                                                                                                                                   cAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGaa                 <  8:258641/149‑1 (MQ=255)
                                                                                                                                    aGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCATGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAAt                <  6:129605/149‑1 (MQ=255)
                                                                                                                                    aGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAAt                <  3:23343/149‑1 (MQ=255)
                                                                                                                                      gcgcTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAAtgt              <  2:249686/149‑1 (MQ=255)
                                                                                                                                        gcTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATgtgt            <  2:291260/149‑1 (MQ=255)
                                                                                                                                         ccggccAGCATGTTCGGCACGAACTTGCTGGTCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTc                                                                   >  2:219551/5‑93 (MQ=25)
                                                                                                                                         ccggccAGCATGTTCGGCACGAACTTGCTGGTCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTc                                                                   <  1:219551/89‑1 (MQ=25)
                                                                                                                                            gccAGCATGTTCGGCACGAACTTGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGa                                                     <  4:284091/103‑1 (MQ=38)
                                                                                                                                            gccAGCATGTTCGGCACGAACTTGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGa                                                     >  3:284091/2‑104 (MQ=38)
                                                                                                                                              cAGCATGTTCGGCACGAACTTGCTGGCCGGTTGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAaga                                      >  8:86637/1‑117 (MQ=255)
                                                                                                                                              cAGCATGTTCGGCACGAACTTGCTGGCCGGTTGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAaga                                      <  7:86637/117‑1 (MQ=255)
                                                                                                                                              cAGCATGTTCGGCACGAACTTGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGGTGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGt                                 <  8:243555/122‑1 (MQ=255)
                                                                                                                                              cAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTcggcgg      >  4:95058/1‑149 (MQ=255)
                                                                                                                                                gCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTcggcggcg    <  5:232431/149‑1 (MQ=255)
                                                                                                                                                gCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTcggcggcg    >  3:10692/1‑149 (MQ=255)
                                                                                                                                                 cAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTcggcggcgg   >  2:66451/1‑149 (MQ=255)
                                                                                                                                                  aaGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCggcggcggc  <  7:181468/149‑1 (MQ=255)
                                                                                                                                                   aGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCggcggcggc  >  8:70112/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGTGGCACCGTCGACCAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGC  >  NC_002947/74012‑74306

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: