Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,898,637 A→T 5.4% intergenic (+118/‑90) pyk → / → PP_4302 pyruvate kinase/urea transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,898,6370AT 5.4% 564.6 / 14.4 167intergenic (+118/‑90)pyk/PP_4302pyruvate kinase/urea transporter
Reads supporting (aligned to +/- strand):  ref base A (89/68);  new base T (7/2);  total (96/71)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 174 total aligned reads displayed.

TAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACG  >  NC_002947/4898489‑4898784
                                                                                                                                                    |                                                                                                                                                   
tAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGTCAAAATGCCGCTCCCCAAGGGGCACGCTGGCAGAAGGGGATTCAGTCCCTGTGGGTGGGGCCTTGTGTCGCGATTTGTGGGCAAAg                                                                                                                                                    >  1:1240646/1‑150 (MQ=255)
      gCTGCGGATCGAGAATTTGATCTGAGATGGCCGGGGGCGCGTTGCGCCCCTTTCGCGAAACACGGCCGCTCCCACCGGGCCACGCTGGCAAACGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGACGCGATTGGAGGGCAAAGccctcc                                                                                                                                              <  1:2191804/150‑1 (MQ=255)
          cGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCTTTGccctccaatc                                                                                                                                          <  2:5385128/150‑5 (MQ=255)
          cGGAGCGAGACTTTGATCTGAGATTGCGGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCCCCCACCGGGCCACGCTGGCACAAGGGAAGGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGa                                                                                                                                          <  1:450692/150‑1 (MQ=39)
              tCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACCCGGGCACAAGGGAAAGCAGTCCCGGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCAGAATAc                                                                                                                                      >  2:5182080/1‑150 (MQ=39)
               cGAGACTTTTATCTGAGATTGCGGGGGGCGCGTTGCGCCCCTTTCGCGACACACGGCCGCTCCCACTGGGCCACGCTGGCACATGGGATTGCATTCCCTGTGGGTGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATAcc                                                                                                                                     <  1:1748259/150‑1 (MQ=255)
                gagaCTTTGATGTGTGTGTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACACGGCCGCTCCCACCGGGCCACGCTGGCACAAGGGAATGCTGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCg                                                                                                                                    <  2:4241474/150‑1 (MQ=38)
                gagaCTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCGTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCg                                                                                                                                    >  2:1340548/1‑150 (MQ=255)
                 agaCTTAGATCTGAGAATGCCGGGGGCGCGTTGCGCCCCTAACGCGAAAAAAAGCCGCACCCACACGGCCACCCAGGCACCAGGGAATACAGTCCCGGTGGGCGCGGCCATGTGAACCGATTGGAGGGAAAATCCCACCCGGAAAAcccc                                                                                                                                   >  2:500736/1‑148 (MQ=255)
                   aCTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGctc                                                                                                                                 <  2:1713658/150‑1 (MQ=255)
                    cTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGctct                                                                                                                                <  2:1316077/150‑1 (MQ=255)
                     tttGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTa                                                                                                                               >  2:4931921/1‑150 (MQ=255)
                            tGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCACTACCCGAcc                                                                                                                        >  2:4559436/1‑150 (MQ=255)
                            tGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAAAGCAGTCCCTGTGGGAGCGGCCATGTGTCGCGATTGGAGGGCCAAGCCCACCCGGAATAACGCTCTACCCGAcc                                                                                                                        >  1:4659046/1‑150 (MQ=38)
                              agaTTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACcgc                                                                                                                      <  2:4333942/150‑1 (MQ=255)
                               gaTTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAAAACCGCTCTACCCGACcgcg                                                                                                                     >  2:3172846/1‑150 (MQ=255)
                                aTTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCgg                                                                                                                    <  2:5537249/150‑1 (MQ=255)
                                aTTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCGTGTGTCGCGATTGGAGGGCAAAGCCCTCCCGGAATACCGCTCTACCCGACCGCgg                                                                                                                    >  2:4348724/1‑150 (MQ=255)
                                      gggggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGACTGCAGTCCCTGTGGGAGCGGCCGTGTGTCGCGATTGGAGGGCCAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAAcc                                                                                                              >  2:2728337/1‑150 (MQ=255)
                                      gggggCGCGTTGCGCCCCATTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCAAAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAAcc                                                                                                              <  2:2322895/150‑1 (MQ=255)
                                      aggggCGCAACGCGCCCCTTTCGCGACACAAGGCCGCTCCCCCAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTTGAGGGCAAAGCCCTCCCTGAATAACGCTCTACCCGACCGCGGAAAAcc                                                                                                              >  2:563873/2‑150 (MQ=39)
                                        gggCGCGTTGCGCCCCTTTCGCGACACACGCCCGCTCCCACTGGGCCACGCTGGAACAAGGGAATGCTGTCCCTGTGGGCGCGGCCTTGTTACGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                                                            <  1:2943988/150‑1 (MQ=255)
                                        gggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACACGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                                                            <  2:873176/150‑1 (MQ=255)
                                        gggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGa                                                                                                            <  2:2302406/150‑1 (MQ=255)
                                         ggCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAg                                                                                                           <  2:5095459/150‑1 (MQ=255)
                                            gcgTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCGGGCACAAGGGAATGCAGTCCCGGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCCACCGCGGAAAACCGAgggc                                                                                                        >  2:5280558/1‑148 (MQ=255)
                                                ggcgcCCTTTTCGCGCCACACGGCCGCCCCCCCCGGGCCCCGCTGGAACATGGTAATGCTGTCCGTGTGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCaaa                                                                                                    <  1:4756384/149‑1 (MQ=255)
                                                cgcgcCCCTTGCGAGACCCAAGCCCCCACCCCAAGGCCCACCCTGGCACACCGGATTGCAGCCCCTCGGGCAGCGGCCTTGTGTCGCGATTGGAGGGCAAACCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCaaa                                                                                                    <  1:4880349/149‑1 (MQ=255)
                                                 gcgcCCCTTTCGCGACACACGGCCGCTCCCACTGGGCCCCGTGGGCAAACGGGATTGCTGGCCTTTTGGGAGGGGCCTTGTTTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAag                                                                                                   <  1:161506/150‑1 (MQ=255)
                                                    ccccTTTCGCGACACAATGCCGCTCCCACAGGGCCACGCTGGGACAAGGGAATGCAGTCCCTGTGGGTGCGGCCGTGTGTGGCGATTGGAGGGCAAAGCCCTCCCTGAATAACGCTATACCCCTCCGCGGAAAACCGAGGCCCAAATAgg                                                                                                >  2:3080654/1‑150 (MQ=255)
                                                      ccTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGTAGCGGCCTTGTGTCGCGATTGGAGGGCAAGGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGcc                                                                                              <  2:2275478/150‑1 (MQ=255)
                                                      ccTTTCGCGACACAAGGCCGCACCCACAGGGCCACGCGGGCACAAGGGAAAGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCACCCTGAATACCGCTCTACCCGAACGCGGAAAACCGAGGCCCAAAGAGGcc                                                                                              >  2:1078677/1‑150 (MQ=39)
                                                        tttCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCa                                                                                            >  1:4533373/1‑150 (MQ=255)
                                                         ttCGCGACACCCGGCCGCTCCCACCGGGCCACGCGGGCACATGGGAATGCAGTCCGTGTGGGCGCGGCCTTGTGGCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCAc                                                                                           <  1:3951779/150‑1 (MQ=255)
                                                          ccgcgacacaaggccgctcccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCAcc                                                                                          >  1:3232144‑M2/80‑150 (MQ=255)
                                                              cACACAGGGCCGCGCCCAGCGGGCCACGCTGTCACTCGGTATTGCCGCCCCTGTGGGTGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTc                                                                                      >  1:2351079‑M1/2‑150 (MQ=255)
                                                                   acGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                                                                 >  2:2949429‑M1/3‑150 (MQ=255)
                                                                   AAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACACGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                                                                 <  1:2833474‑M1/150‑1 (MQ=255)
                                                                   AAGGCCGCTCCCACAGGGCCACGCGGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCa                                                                                 >  2:4074603‑M1/1‑150 (MQ=255)
                                                                     GGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTTTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATc                                                                               >  2:2491545‑M1/1‑150 (MQ=255)
                                                                     GGCCGCTCCCACAGGGCCACGCTGGCAAAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATc                                                                               >  2:2491519‑M1/1‑150 (MQ=255)
                                                                     GGCCGCGCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCGTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATc                                                                               >  1:3202030‑M1/1‑150 (MQ=255)
                                                                         gctcccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAAcc                                                                           >  1:2552220‑M2/65‑150 (MQ=255)
                                                                           TCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCaaccac                                                                         <  2:3740679‑M1/150‑3 (MQ=255)
                                                                           TCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTc                                                                         <  1:3422522‑M1/150‑1 (MQ=255)
                                                                           TCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTc                                                                         <  1:1032711‑M1/150‑1 (MQ=255)
                                                                             ccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACTGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCCTCAACCTCCa                                                                       >  2:4691355‑M2/61‑150 (MQ=255)
                                                                             ccacagggactgcattcccttgtgccagcgtggccctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACTGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCa                                                                       >  2:4691462‑M2/61‑150 (MQ=255)
                                                                             CCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCa                                                                       <  1:4197554‑M1/150‑1 (MQ=255)
                                                                              CACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                                                      <  1:4931911‑M1/150‑1 (MQ=255)
                                                                              CACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCACTGAATACCGCTATACCCGACCGCGGACAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCcaa                                                                      <  1:1881360‑M1/150‑2 (MQ=255)
                                                                               ACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCact                                                                     <  2:3785641‑M1/150‑1 (MQ=255)
                                                                               ACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                                                      <  1:77718‑M1/149‑1 (MQ=255)
                                                                               ACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                                                      <  1:77701‑M1/149‑1 (MQ=255)
                                                                               ACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                                                      >  2:77701‑M1/1‑149 (MQ=255)
                                                                               ACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCac                                                                      >  2:77718‑M1/1‑149 (MQ=255)
                                                                                  GGGCCAGGCTGGCACATGGTTTTGCTGTCCGTGGGGGTGGGGCGTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCactgac                                                                  >  1:514636‑M1/1‑150 (MQ=255)
                                                                                  GGGCCACGCTGGCACTCGGGATTGCAGTCCGTGTGGGAGGGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCactgac                                                                  >  1:3092189‑M1/1‑150 (MQ=255)
                                                                                  GGGCCACGCTGCCACAAGGTAAGTCTGTCCGTGTTGGAGGGGCTTTTTTTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCactgac                                                                  >  1:193896‑M1/1‑150 (MQ=255)
                                                                                      cACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGcc                                                              >  1:4550554/1‑150 (MQ=255)
                                                                                       aCGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGccc                                                             <  2:5111945/150‑1 (MQ=255)
                                                                                         gCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCa                                                           >  1:4437145/1‑150 (MQ=255)
                                                                                             gCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGTCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                       >  2:1427029/1‑150 (MQ=255)
                                                                                             gCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTa                                                       <  2:3101198/150‑1 (MQ=255)
                                                                                              cacaAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAATCCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAAGCTCCACTGACTGCCCCATTTac                                                      >  2:5193331/1‑150 (MQ=255)
                                                                                               acaTGGGAATGCAGTCCGGGTGGGTGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     <  1:1585372/150‑1 (MQ=255)
                                                                                               acaAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     >  1:3345608/1‑150 (MQ=255)
                                                                                                 aaGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACAcc                                                   >  2:2740188/1‑150 (MQ=255)
                                                                                                  aGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCa                                                  >  2:4621665/1‑150 (MQ=255)
                                                                                                  aGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCa                                                  >  2:3598317/1‑150 (MQ=255)
                                                                                                  aGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCa                                                  >  2:1822965/1‑150 (MQ=255)
                                                                                                   gggAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCaa                                                 >  2:1957806/1‑150 (MQ=255)
                                                                                                    ggAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCaaa                                                >  1:5026981/1‑150 (MQ=255)
                                                                                                      aaTGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCaaaaa                                              >  1:3553104/1‑150 (MQ=255)
                                                                                                         gCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAAttt                                           <  2:4456602/150‑1 (MQ=255)
                                                                                                         gCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACAGCCCCATGTACACCAAAAAttt                                           >  2:2540523/1‑150 (MQ=255)
                                                                                                          cAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTc                                          <  2:1607600/150‑1 (MQ=255)
                                                                                                           aGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCg                                         >  1:4566004/1‑150 (MQ=255)
                                                                                                           aGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCg                                         >  1:3523408/1‑150 (MQ=255)
                                                                                                            gTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGt                                        >  1:5265988/1‑150 (MQ=255)
                                                                                                               cctgtgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCaa                                     >  1:932722‑M2/27‑150 (MQ=255)
                                                                                                                   tgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     >  2:3553812‑M2/23‑130 (MQ=255)
                                                                                                                   tgggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTaca                                                     <  1:3553818‑M2/108‑1 (MQ=255)
                                                                                                                    gggagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGt                                <  2:932722‑M2/129‑1 (MQ=255)
                                                                                                                      gagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCGAAAATTTCGTCAACCCGTGc                              <  2:3167235‑M2/131‑1 (MQ=255)
                                                                                                                      gagcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGc                              <  2:4721297‑M2/131‑1 (MQ=255)
                                                                                                                       agcggccttgtgtcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGcc                             >  2:5004843‑M2/19‑150 (MQ=255)
                                                                                                                                   tcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCa                 >  1:1809204‑M2/7‑150 (MQ=255)
                                                                                                                                   tcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCa                 >  1:1566557‑M2/7‑150 (MQ=255)
                                                                                                                                   tcgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCa                 >  1:2745877‑M2/7‑150 (MQ=255)
                                                                                                                                    cgcgaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCAc                <  2:630872‑M2/145‑1 (MQ=255)
                                                                                                                                       gaTTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAATCCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGc             >  1:5156093‑M2/3‑150 (MQ=255)
                                                                                                                                         agggAGGGCTTTGCCCTCCCTGAATACCTCTCAACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTACACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGt           >  1:2322126/3‑150 (MQ=255)
                                                                                                                                          gggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGtt          >  2:3654152/2‑150 (MQ=255)
                                                                                                                                           ggAGGGCTTTGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTa         >  2:5399291/1‑150 (MQ=255)
                                                                                                                                             aGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACt       >  2:242704/1‑150 (MQ=255)
                                                                                                                                               ggCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCaaa                                                                                                    <  2:3308818/55‑1 (MQ=255)
                                                                                                                                                  aaaGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACg  >  1:750369/1‑150 (MQ=255)
                                                                                                                                                  aaaGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACg  >  1:3156091/1‑150 (MQ=255)
                                                                                                                                                  aaaGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACg  >  1:2082493/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TAACACGCTGCGGATCGAGACTTTGATCTGAGATTGCCGGGGGCGCGTTGCGCCCCTTTCGCGACACAAGGCCGCTCCCACAGGGCCACGCTGGCACAAGGGAATGCAGTCCCTGTGGGAGCGGCCTTGTGTCGCGATTGGAGGGCAAAGCCCTCCCTGAATACCGCTCTACCCGACCGCGGAAAACCGAGGCCCAAAGAGGCCCACCGCTCCCCCATCAACCTCCACTGACTGCCCCATGTACACCAAAAATTTCGTCAACCCGTGCCCCGACTGGGCCACGGCGTTACTCAACG  >  NC_002947/4898489‑4898784

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: