Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,183,046 T→A 5.4% K246* (AAG→TAG) ‡ PP_4563 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,183,0460TA 5.4% 230.6 / 3.5 74K246* (AAG→TAG) ‡PP_4563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/35);  new base A (2/2);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGC  >  NC_002947/5182899‑5183189
                                                                                                                                                    |                                                                                                                                               
gTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGc                                                                                                                                               >  1:4764/1‑150 (MQ=255)
 tGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCt                                                                                                                                              >  2:212039/1‑150 (MQ=255)
 tGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCt                                                                                                                                              <  1:648877/150‑1 (MQ=255)
   gCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                            <  2:790020/150‑1 (MQ=255)
   cgcagcagcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                            <  1:96937/147‑1 (MQ=255)
     gAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATc                                                                                                                                          >  2:260920/1‑150 (MQ=255)
        cagcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                       <  1:647817/150‑1 (MQ=255)
           tccgatctCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAAT‑TCCTTGGGGTAATTTGCCCTTGCTGACGTTGAACAGTTTGCTGATCAACGCa                                                                                                                                    <  1:863660/143‑1 (MQ=255)
                 ttacGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCa                                                                                                                              >  1:175856/4‑150 (MQ=255)
                    cGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                            <  2:535016/133‑1 (MQ=255)
                    cGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                            >  1:535016/1‑133 (MQ=255)
                    cGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActc                                                                                                                           >  2:666678/1‑150 (MQ=255)
                     gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAAGAACCCGCCAATGGCATTCTTG‑AACTTGACATTGTTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActct                                                                                                                          >  1:725936/1‑150 (MQ=255)
                     gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActct                                                                                                                          >  1:1238631/1‑150 (MQ=255)
                     gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActct                                                                                                                          >  2:295109/1‑150 (MQ=255)
                       aaCTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACgtca                                                                                                                                   <  1:412065/139‑4 (MQ=255)
                       aaCTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACgtca                                                                                                                                   >  2:412065/1‑136 (MQ=255)
                         cTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTc                                                                                                                      <  2:852916/150‑1 (MQ=255)
                           tCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGc                                                                                                                    <  2:411264/150‑1 (MQ=255)
                                cAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACAcc                                                                                                               <  1:713139/150‑1 (MQ=255)
                                  ggCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACAcccc                                                                                                             <  1:1173158/150‑1 (MQ=255)
                                         ggCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACAcc                                                                                                      >  1:379908/1‑150 (MQ=255)
                                         ggCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACAcc                                                                                                      <  1:763647/150‑1 (MQ=255)
                                            cGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCg                                                                                                   >  2:951351/1‑150 (MQ=255)
                                              tCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCtt                                                                                                                         <  1:447868/126‑1 (MQ=255)
                                              tCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCtt                                                                                                                         >  2:447868/1‑126 (MQ=255)
                                               cAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCt                                                                                                <  2:4764/150‑1 (MQ=255)
                                                aGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCtt                                                                                               <  1:666678/150‑1 (MQ=255)
                                                aGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCtt                                                                                               >  1:1006736/1‑150 (MQ=255)
                                                  gCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCg                                                                                             <  2:175856/150‑1 (MQ=255)
                                                   cTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCgg                                                                                            >  2:292895/1‑150 (MQ=255)
                                                    tttCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGc                                                                                           <  1:546593/150‑1 (MQ=255)
                                                    tttCTCGAAGCGATGGTCATGGAAAACGATCAGGCTCTGGTTGAAAAACCGGCCAATGGCCTACTTG‑AAATTGCCGTTGCTGACGTTGAACAGCTTGATGATCAACGCAGGGTAACTCTTGTCGCACACCCCCCACACCCCTCCTTCCGc                                                                                           >  2:199287/1‑150 (MQ=255)
                                                    tgtgtcGATGCGATCGTCAAGGACGACGATCAGGCCCTGCTTAAACAACCCGCCAATGGCCCTCTTG‑AAGTTGCCCTTGCTGACGTTGCACAGCTTGCTGATCAACGCACTGTCACTCTTGTCGCACACCCCCAACACTCCGCCTTCGGc                                                                                           <  2:583657/146‑1 (MQ=255)
                                                        tcGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTcc                                                                                       >  2:334565/1‑150 (MQ=255)
                                                           gcctgATCGTCATCCATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActc                                                                                                                           >  2:285835/5‑111 (MQ=255)
                                                           gcctgATCGTCATCCATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActc                                                                                                                           <  1:285835/107‑1 (MQ=255)
                                                               gATCGTCATCCATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAggg                                                                                                                                 <  2:110167/101‑1 (MQ=255)
                                                               gATCGTCATCCATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAggg                                                                                                                                 >  1:110167/1‑101 (MQ=255)
                                                                  cGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCa                                                                             <  2:1006736/150‑1 (MQ=255)
                                                                    tCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATg                                                                           <  1:260920/150‑1 (MQ=255)
                                                                      aTGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGAt                                                                         <  1:951351/150‑1 (MQ=255)
                                                                           tGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGct                                                                    <  2:624153/150‑1 (MQ=255)
                                                                            tccGAGCAGTCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCACCGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGctc                                                                   <  1:1180321/148‑1 (MQ=255)
                                                                              cGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCtt                                                                 <  1:212039/150‑1 (MQ=255)
                                                                                   aGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAgg                                                            <  2:379908/150‑1 (MQ=255)
                                                                                    ggCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCt                                                                                          >  2:34987/1‑119 (MQ=255)
                                                                                    ggCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCt                                                                                          <  1:34987/119‑1 (MQ=255)
                                                                                     gCCCTTCTTGAACAACCCGCTAATGGTCTTCTTG‑AAGTTGCCCTTGCTGTCGTTTAACAGCTTGCTGATCAACGCATGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGgt                                                          <  1:78714/150‑1 (MQ=255)
                                                                                     gCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAAc                                                                                                                                       <  1:177519/73‑1 (MQ=255)
                                                                                     gCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAAc                                                                                                                                       >  2:177519/1‑73 (MQ=255)
                                                                                              gAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGACACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAg                                                 >  1:224853/1‑150 (MQ=255)
                                                                                                    cccGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGggcggc                                           >  1:787682/1‑150 (MQ=255)
                                                                                                    cccGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGggcggc                                           <  1:334565/150‑1 (MQ=255)
                                                                                                    cccGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGggcggc                                           <  1:982578/150‑1 (MQ=255)
                                                                                                       gCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGcc                                        >  1:1148243/1‑150 (MQ=255)
                                                                                                               ccTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCCGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCt                                <  1:921204/150‑1 (MQ=255)
                                                                                                                cttcttG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATAAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTATTGCAGGGGTTGCGCCAGGGGGGCGCCTACCGGCTg                               >  1:496812/1‑150 (MQ=255)
                                                                                                                    ttG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTa                                      <  2:768938/139‑1 (MQ=255)
                                                                                                                    ttG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTa                                      >  1:768938/1‑139 (MQ=255)
                                                                                                                        aacttGCCCTTGCTCACATTTAACAGCTTGTTGATCACCGCATGGTCACTCTTGTCTCACACCCCCAACACCCCCCCTTCGACTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTc                        <  2:496812/147‑1 (MQ=255)
                                                                                                                         aGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCa                       <  1:813409/150‑1 (MQ=255)
                                                                                                                          gTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAg                      >  1:77545/1‑150 (MQ=255)
                                                                                                                          gTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGAGGCGCCTACCGGCTGCAGGCTCAg                      >  1:371582/1‑150 (MQ=255)
                                                                                                                              cccTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGaa                  >  2:691880/1‑150 (MQ=255)
                                                                                                                               ccTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGaaa                 >  2:98374/1‑150 (MQ=255)
                                                                                                                                  accctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGTGCGCCATGATCTGCTCTTGCAGGTTGTCGGCCAGGGCGGGGCCTACCGGCTGCAGGCTCAGTGAActct              >  1:169228/7‑150 (MQ=255)
                                                                                                                                     tGACGTTGTACAGCTTGCTGATCAAAGCAGGGTCACTCTTGTCGCACACCCCCAAGACCCCGCCTTCGGCTTCCAGGCGCGCCATGATGTGCTCTTGCAGGGTGTCCGCCGGGGCGGAGCCTACCGGCTTAAGGCTCTGTGAAATCTTgg           >  2:162706/1‑149 (MQ=255)
                                                                                                                                     tGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCa                                                              <  2:423313/99‑1 (MQ=255)
                                                                                                                                     tGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCa                                                              >  1:423313/1‑99 (MQ=255)
                                                                                                                                       aCGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCg         >  2:334392/1‑150 (MQ=255)
                                                                                                                                         gTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTc       >  1:1177485/1‑150 (MQ=255)
                                                                                                                                             aaCAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTgg   >  2:1030999/1‑150 (MQ=255)
                                                                                                                                              aCAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGgc  >  1:1175956/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGC  >  NC_002947/5182899‑5183189

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑