Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,183,051 A→T 5.9% I244N (ATC→AAC)  PP_4563 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,183,0510AT 5.9% 205.1 / 3.3 68I244N (ATC→AAC) PP_4563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base T (2/2);  total (35/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 = 9.55e-01

GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCT  >  NC_002947/5182902‑5183198
                                                                                                                                                      |                                                                                                                                                   
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                                                                                                                                                     >  1:535016/1‑133 (MQ=255)
                 cGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                                     <  2:535016/133‑1 (MQ=255)
                 cGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActc                                                                                                                                    >  2:666678/1‑150 (MQ=255)
                  gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAAGAACCCGCCAATGGCATTCTTG‑AACTTGACATTGTTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActct                                                                                                                                   >  1:725936/1‑150 (MQ=255)
                  gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActct                                                                                                                                   >  2:295109/1‑150 (MQ=255)
                  gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActct                                                                                                                                   >  1:1238631/1‑150 (MQ=255)
                    aaCTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACgtca                                                                                                                                            >  2:412065/1‑136 (MQ=255)
                    aaCTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACgtca                                                                                                                                            <  1:412065/139‑4 (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:763647/150‑1 (MQ=255)
                                      ggCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACAcc                                                                                                               >  1:379908/1‑150 (MQ=255)
                                         cGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCg                                                                                                            >  2:951351/1‑150 (MQ=255)
                                           tCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCtt                                                                                                                                  >  2:447868/1‑126 (MQ=255)
                                           tCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCtt                                                                                                                                  <  1:447868/126‑1 (MQ=255)
                                            cAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCt                                                                                                         <  2:4764/150‑1 (MQ=255)
                                             aGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCtt                                                                                                        >  1:1006736/1‑150 (MQ=255)
                                             aGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCtt                                                                                                        <  1:666678/150‑1 (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:334565/150‑1 (MQ=255)
                                                                                                 cccGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGggcggc                                                    <  1:982578/150‑1 (MQ=255)
                                                                                                 cccGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGggcggc                                                    >  1:787682/1‑150 (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)
                                                                                                                                                    tGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCt  >  1:1199454/1‑150 (MQ=255)
                                                                                                                                                      aTCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCATGGTGTcc                                                               <  2:1000966/87‑1 (MQ=255)
                                                                                                                                                      aTCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCATGGTGTcc                                                               >  1:1000966/1‑87 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTG‑AAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTCCT  >  NC_002947/5182902‑5183198

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