Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 74,166 G→C 7.4% V214V (GTG→GTC) ‡ 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,1660GC 7.4% 214.7 / 5.4 67V214V (GTG→GTC) ‡PP_0063lipid A biosynthesis lauroyl acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base C (3/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

TCGACCAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAG  >  NC_002947/74023‑74312
                                                                                                                                               |                                                                                                                                                  
tCGACCAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGaa                                                                                                                                               <  5:103961/149‑1 (MQ=255)
   aCCAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACAt                                                                                                                                            >  8:330706/1‑149 (MQ=255)
     cAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGc                                                                                                                                          <  8:165642/149‑1 (MQ=255)
        ggaaggCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTgg                                                                                                                                       <  6:14897/149‑1 (MQ=255)
        ggaaggCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTgg                                                                                                                                       >  6:119102/1‑149 (MQ=255)
        ggaaggCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGCGGCCAGGTGGGTATTCATGCCAACCCGGAGCCGGCGGATTCGGCGGGCGTGTTCTTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAATTTCGTGCGGAAAATGATgg                                                                                                                                       >  7:51673/1‑149 (MQ=255)
          aaggCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGcc                                                                                                                                     <  5:316306/149‑1 (MQ=255)
                tCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTgg                                                                                                                                >  6:25352/1‑148 (MQ=255)
                     aGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAgg                                                                                                                           <  3:194746/148‑1 (MQ=255)
                           aTCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAg                                                                                                                     <  7:330706/148‑1 (MQ=255)
                           aTCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAg                                                                                                                     <  7:157321/148‑1 (MQ=255)
                           aTCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAg                                                                                                                     >  8:157321/1‑148 (MQ=255)
                            tCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAgg                                                                                                                           >  4:46902/1‑141 (MQ=255)
                            tCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAgg                                                                                                                           <  3:46902/141‑1 (MQ=255)
                            tCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgt                                                                                                                   <  5:25352/149‑1 (MQ=255)
                                ggAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgtgtgt                                                                                                               >  8:263330/1‑149 (MQ=255)
                                    gTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCt                                                                                                           <  3:172336/149‑1 (MQ=255)
                                    gTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCt                                                                                                           >  4:172336/1‑149 (MQ=255)
                                               tggCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgt                                                                                                                   <  5:42605/130‑1 (MQ=255)
                                               tggCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgt                                                                                                                   >  6:42605/1‑130 (MQ=255)
                                                  ccAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCATGTTCGGCACGAACTTGCTg                                                                                                                                        <  4:263052/106‑1 (MQ=255)
                                                  ccAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCATGTTCGGCACGAACTTGCTg                                                                                                                                        >  3:263052/1‑106 (MQ=255)
                                                    aGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTg                                                                                           >  6:139438/1‑149 (MQ=255)
                                                      gTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGcc                                                                                         <  6:137114/149‑1 (MQ=255)
                                                             ttCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATgg                                                                                   <  5:119102/148‑1 (MQ=255)
                                                                cTGCCGACCCGGATCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGAt                                                                                     <  7:12233/143‑1 (MQ=255)
                                                                cTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCg                                                                               <  2:47769/149‑1 (MQ=255)
                                                                cTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCCTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGAt                                                                                     >  8:12233/1‑143 (MQ=255)
                                                                 tGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGc                                                                                                                                      <  8:95085/93‑1 (MQ=255)
                                                                 tGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGc                                                                                                                                      >  7:95085/1‑93 (MQ=255)
                                                                   ccGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGc                                                                            <  5:139438/149‑1 (MQ=255)
                                                                      aCCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgt                                                                                                                   >  7:12560/1‑107 (MQ=255)
                                                                      aCCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgt                                                                                                                   <  8:12560/107‑1 (MQ=255)
                                                                        ccGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTc                                                                         >  8:186417/1‑147 (MQ=255)
                                                                        ccGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACAGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTc                                                                         >  8:182904/1‑147 (MQ=255)
                                                                         cGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGCCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAg                                                                      <  7:186417/149‑1 (MQ=255)
                                                                                     aaTCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCaa                                                                       >  1:250567/1‑136 (MQ=255)
                                                                                     aaTCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCaa                                                                       <  2:250567/136‑1 (MQ=255)
                                                                                         ggcgggcgTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCTGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTg                                                      >  8:171967/1‑149 (MQ=255)
                                                                                                   ttCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgt                                                                                                                   <  6:270819/78‑1 (MQ=255)
                                                                                                   ttCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGgt                                                                                                                   >  5:270819/1‑78 (MQ=255)
                                                                                                     cGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAagag                                           <  7:182904/148‑1 (MQ=255)
                                                                                                      gTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGAt                                         >  7:245968/1‑149 (MQ=255)
                                                                                                       tACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATg                                        >  1:67564/1‑149 (MQ=255)
                                                                                                           gTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTgc                                                     <  8:98620/132‑1 (MQ=255)
                                                                                                           gTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTgc                                                     >  7:98620/1‑132 (MQ=255)
                                                                                                             tCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGGATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTc                                                                >  1:111683/1‑119 (MQ=255)
                                                                                                             tCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTc                                                                <  2:111683/119‑1 (MQ=255)
                                                                                                              ccTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCAACCGCAGCGCATGCAGGaa                                                                                     <  6:342331/97‑2 (MQ=255)
                                                                                                              ccTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCAACCGCAGCGCATGCAGGaa                                                                                     >  5:342331/1‑96 (MQ=255)
                                                                                                               cTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCAc                                >  1:63881/1‑149 (MQ=255)
                                                                                                                tCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCAcc                               >  2:77505/1‑149 (MQ=255)
                                                                                                                tCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGACCGTGGCAAGCCGGTAGGTGAGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTAAAGCAcc                               >  2:250719/1‑149 (MQ=255)
                                                                                                                   gCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGAt                            <  2:63881/149‑1 (MQ=255)
                                                                                                                   gCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGAt                            <  8:245968/149‑1 (MQ=255)
                                                                                                                        cAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGaa                       <  1:77505/149‑1 (MQ=255)
                                                                                                                         aGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAAt                      <  7:263330/149‑1 (MQ=255)
                                                                                                                            caccggccAGCATGTTCGGCACGAACTTGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAAtgtg                   >  2:74918/7‑149 (MQ=255)
                                                                                                                                 aCCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGGGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGAt                            >  4:216856/1‑135 (MQ=255)
                                                                                                                                   cAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAAt                      >  7:240026/1‑139 (MQ=255)
                                                                                                                                   cAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTcggcgg            >  5:126357/1‑149 (MQ=255)
                                                                                                                                   cAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGGGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAAt                      <  8:240026/139‑1 (MQ=255)
                                                                                                                                    aGCATGTTCGGCACGAACTTGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGAt                            >  6:19782/1‑132 (MQ=255)
                                                                                                                                    aGCATGTTCGGCACGAACTTGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGAt                            <  5:19782/132‑1 (MQ=255)
                                                                                                                                      cAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTcggcggcgg         >  4:286956/1‑149 (MQ=255)
                                                                                                                                             gTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGa   <  3:286956/148‑1 (MQ=255)
                                                                                                                                              tGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAg  >  5:357332/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TCGACCAAGGAAGGCATCCTCAGCGTGATCAAGGAAGTGCGCCGTGGTGGCCAGGTGGGTATTCCTGCCGACCCGGAGCCGGCGGAATCGGCGGGCGTGTTCGTACCGTTCCTCGGCACCCAGGCGCTGACCAGCAAGTTCGTGCCGAACATGCTGGCCGGTGGCAAGGCGGTAGGTGTGTTCCTGCATGCGCTGCGGTTGCCGGATGGCTCGGGCTTCAAGGTGTTCCTCGAAGCTGCGCCGGAAGAGATGTACAGCACCGATGTGAATGTGTCGGCGGCGGCGATGAG  >  NC_002947/74023‑74312

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