Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,632,711 G→A R454C (CGT→TGT)  guaB ← inosine 5'‑monophosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,632,7110GA98.5% 214.8 / NA 69R454C (CGT→TGT) guaBinosine 5'‑monophosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/35);  minor base C (1/0);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACC  >  NC_000913/2632467‑2632957
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
aaCCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATGCCGAGGCAAGTGAAACAGATAATATAAAGCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTGCTGAATGCCCGCACCGCTTAGACGTACAAACGCCGCTTCAGTACACAGTtc                                                                                                                                                                                                                                                  <  2:72618/251‑1 (MQ=255)
aaCCAGTTGAGTGTAATGAGAACCGAAGTCCAGAATGAGTATGGGATTCTTATGAATGTTTTCCGTCTTTGACGCGTATTCCGAGGCAAGTGAATCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCCGTAGCCCAGACGGTAGTGCGGGGACTCTTTAGTAAGGGTCACGTCGTGAACGTGGCTTTCCTGAATTCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTtc                                                                                                                                                                                                                                                  <  1:624675/251‑1 (MQ=255)
    aGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTtcgtcg                                                                                                                                                                                                                                              >  1:306535/1‑251 (MQ=255)
          gtgtACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTa                                                                                                                                                                                                                                        <  2:510362/251‑1 (MQ=255)
            gtACTGAGAACAGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGAAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTAcc                                                                                                                                                                                                                                      >  2:453880/1‑251 (MQ=255)
              aCTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGAGAATATAAAGCGCCCGACATGAAGTCGGGCTAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACTTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACca                                                                                                                                                                                                                                     <  1:331559/250‑1 (MQ=255)
                ggagaACCGAAGTCCAAAATGGGGATGCGATGCTTATGAAAGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGAAAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAg                                                                                                                                                                                                                                  <  2:177834/250‑1 (MQ=255)
                   gaACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCg                                                                                                                                                                                                                               >  1:178637/1‑251 (MQ=255)
                              cAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCCTTGACGCTTATTCCGAGGCGAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTACGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGAcccct                                                                                                                                                                                                                    >  2:118379/1‑249 (MQ=255)
                                gAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATAc                                                                                                                                                                                                                  >  2:172484/1‑251 (MQ=255)
                                   tGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAgg                                                                                                                                                                                                               >  1:320334/1‑251 (MQ=255)
                                                 ttATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGccc                                                                                                                                                                                                 >  1:109491/1‑251 (MQ=255)
                                                    tGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATCATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGAGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGAGCAGGCCGCCCATc                                                                                                                                                                                              >  1:228847/1‑251 (MQ=255)
                                                         gTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATctgct                                                                                                                                                                                          <  2:113813/250‑1 (MQ=255)
                                                          ttttCCGTCATTGACGCTTATGCCGAGGAAAGGGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTTGTTAGTAATTGGCACGGCGTGAACGTGGCTTTCCTGAATGCCCGCGCCGCTGATACGTAGAAACTCCGCTTTAGTACACAGTTCGGCGATAGTTCCACAGCCGGTCAGGCCCATACAGGAGCGCATGCCGCCCGTctgctg                                                                                                                                                                                         <  2:318915/250‑1 (MQ=255)
                                                              ccGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGa                                                                                                                                                                                     >  1:167804/1‑250 (MQ=255)
                                                               cGTCATTGACGCTTATTCCGAGTCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCCCCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAAt                                                                                                                                                                                   <  2:204979/251‑1 (MQ=255)
                                                                 tCATTGTCGCTAATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCTGTCAGACCCAGACATGAGCGCAGGCCGCCCATCTGCTGGTGAATg                                                                                                                                                                                  <  2:113913/250‑1 (MQ=255)
                                                                  cATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGCAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCTCAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGa                                                                                                                                                                                 <  1:353602/250‑1 (MQ=255)
                                                                  cATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACGCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGa                                                                                                                                                                                 <  1:353632/250‑1 (MQ=255)
                                                                        cGCTTATTCCGAGGCACGTGAAACAGATACTATAAATCGCCCGACATGAAGTCGGGCGAAGAGCATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGCATGCCCGCACCGCTGATCCGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCttt                                                                                                                                                                          >  2:318306/1‑251 (MQ=255)
                                                                            aaatCCGAGGAAAGGGAAACAGAGAATATAAATCGCCCGACATGAAGGCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCGTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACTAACTCCGCTTTAGTACACAGTTCGTCGATATTTCCACAGCCGGTCAGACCCATTCAGGAGCGCAGGCCGCCCAGCTTCTGGTGAATGATCTCTTTCAg                                                                                                                                                                       <  2:426131/247‑1 (MQ=255)
                                                                              ttCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCg                                                                                                                                                                    >  2:538863/1‑251 (MQ=255)
                                                                                    ggCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTGTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCttt                                                                                                                                                              <  2:48776/251‑1 (MQ=255)
                                                                                     gCAAGTGAAACAGATAATATAAATCGCACGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGCCTCTTTAGTAAGGGTCACGGCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTGTCAGGCGACCTTta                                                                                                                                                             <  1:347160/251‑1 (MQ=255)
                                                                                      cAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATctt                                                                                                                                                                                            <  2:333116/219‑2 (MQ=255)
                                                                                      cAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATctg                                                                                                                                                                                            >  1:333116/1‑219 (MQ=255)
                                                                                       aaGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGCATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTCATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGCTACGAACAAACTCCGCTTTAGTACCCAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTtat                                                                                                                                                            >  1:603657/1‑250 (MQ=255)
                                                                                           gAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGc                                                                                                                                                        <  2:109491/250‑1 (MQ=255)
                                                                                           gAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCt                                                                                                                                                       <  2:599501/251‑1 (MQ=255)
                                                                                             aaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTAc                                                                                                                                                     <  1:417001/251‑1 (MQ=255)
                                                                                               cAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTAcgc                                                                                                                                                   <  1:4678/251‑1 (MQ=255)
                                                                                                aGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTAcgcg                                                                                                                                                  <  1:543268/251‑1 (MQ=255)
                                                                                                   taataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGAcc                                                                                                                                               >  1:6423/1‑251 (MQ=255)
                                                                                                    aataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCt                                                                                                                                              >  2:345888/1‑251 (MQ=255)
                                                                                                         aaaTCGCCCGACATGAAGGCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGAt                                                                                                                                         <  1:232250/251‑1 (MQ=255)
                                                                                                            gctcccGACATGAAGGCGGGCGAAGAGAATCAGGAGCCCAGGCGGTAGTTCGGGGCCTTCTTAGTAAGGGTCACGTCGTGAACGTGGCTTTCCTGAATGCGCGGACCGGTGATACGTACAAACTGCGCTTTAGTTCACAGTTCTTCGATAGTTCCACAGCCGGTCAGATCCGTACAGGAGCGCAGGCCGCCCATCTGCTGGTTAATGAGCTTTTTCAGGCGACCTTGATAGGCTACGCGGCATTCGATAcc                                                                                                                                      <  2:57828/248‑1 (MQ=255)
                                                                                                             cGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCt                                                                                                                                     <  2:613954/251‑1 (MQ=255)
                                                                                                               cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTc                                                                                                                                   <  2:167804/251‑1 (MQ=255)
                                                                                                                  gACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCgg                                                                                                                                >  1:466671/1‑251 (MQ=255)
                                                                                                                        aaGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATAAGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCAcc                                                                                                                            >  1:117564/1‑249 (MQ=255)
                                                                                                                            cGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATCCAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTg                                                                                                                      >  2:306107/1‑251 (MQ=255)
                                                                                                                                     aaaagcAGGAGCCCAGACGGTAGTTCGTGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAATCTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATATCTTCCGGCACCTGTTTGTCGGCatcg                                                                                                             <  2:341577/246‑4 (MQ=255)
                                                                                                                                        aTCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTa                                                                                                          <  1:484531/251‑1 (MQ=255)
                                                                                                                                                cccAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTg                                                                                                  >  2:185740/1‑251 (MQ=255)
                                                                                                                                                 ccAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTgg                                                                                                 >  2:206751/1‑251 (MQ=255)
                                                                                                                                                   aGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGaa                                                                                               <  1:140732/251‑1 (MQ=255)
                                                                                                                                                     aCGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCGGGTCCGACCCATACAGGAGCGCAGGCCGCCCCTCTGCTGGTGAATGATATCTTTCAGGCGACCTTTATAGGCTACGCGAACTTCGATCGCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGCAAt                                                                                             >  1:608960/1‑251 (MQ=255)
                                                                                                                                                       ggTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATaa                                                                                           <  1:448893/251‑1 (MQ=255)
                                                                                                                                                       ggTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGGTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATaa                                                                                           <  2:320334/251‑1 (MQ=255)
                                                                                                                                                        gTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTCCCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCCGCGTTATCGCTCTGGAAATAAc                                                                                          <  2:466671/251‑1 (MQ=255)
                                                                                                                                                              cGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCa                                                                                     <  1:555910/250‑1 (MQ=255)
                                                                                                                                                                      ctTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCt                                                                            >  1:91553/1‑251 (MQ=255)
                                                                                                                                                                         tAGTCATGGTCACGTCGTGGACGGGGCTTTCCTGAATGCCCGCAATGCTGATAAGACCAATCTCCGCTTTAAAACACAGTTCGGCGATAGTACCAAAGCAGGTCCGACTCATACAGAAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGACGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGATCTGGACCTAACGGTCAGAGGAACCTTTg                                                                         >  2:609829/1‑251 (MQ=255)
                                                                                                                                                                           gTAATGGTCACGTCGTGAACGTGGCTTTCCTTAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGAGACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGa                                                                       <  1:204145/251‑1 (MQ=255)
                                                                                                                                                                                      gtcgtGAACTTGGTTTTCCTGAATGCCCGCACCGCTGTTACGTACAAACTCCGCTTTAGTACTCAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATCCAGGAGCGCAGTCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATTCCTTCCGGCACCAGTTTGTCGGCTGCGTTATCGCTCTGGAAATAACGGTCAGAGGTACCTTTGTACATCGCGCCCa                                                            <  1:50885/251‑1 (MQ=255)
                                                                                                                                                                                            aaCGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAAc                                                      <  2:285342/251‑1 (MQ=255)
                                                                                                                                                                                               gTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGCTACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAAccc                                                    >  1:396756/1‑250 (MQ=255)
                                                                                                                                                                                               gTGGCTTTCCTGAAGGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCa                                                   <  2:178637/251‑1 (MQ=255)
                                                                                                                                                                                                          gAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTa                                         <  1:172484/250‑1 (MQ=255)
                                                                                                                                                                                                          gAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATAAAGGAGCGCAGGCAGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACGTTTATAGGCTACGCAACCTTCGATACCTTCCGGCACCAGTTTGTCGGCCGCGTTATCGCTCTGGAAATACCGGTCAGAGGAACCTTTGGACAACGCGCCCAGGGAACCCATACCACggca                                         >  2:523270/1‑248 (MQ=255)
                                                                                                                                                                                                           aaTGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAATTCGTAACTAGTACCACAGCCGGTCACACCCATACAGGAGCGCAGGACGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTCTCCGCTAAGCGCCCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTa                                         >  1:65519/1‑249 (MQ=255)
                                                                                                                                                                                                              gcccgcccCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTAGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGAtt                                    >  1:410458/1‑251 (MQ=255)
                                                                                                                                                                                                                 cGCACCGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGt                                 >  1:494369/1‑251 (MQ=255)
                                                                                                                                                                                                                     ccGCTGATACGTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATATCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTa                                >  1:29254/1‑248 (MQ=255)
                                                                                                                                                                                                                        cTGATACGTACAAACTCCGCTTTAGTACACCGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACg                          >  1:565511/1‑251 (MQ=255)
                                                                                                                                                                                                                               gTACAAACTCCGCTTTAGTACACAGTTCGTCGATAGTACCAAAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGAGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTgg                   >  1:216921/1‑251 (MQ=255)
                                                                                                                                                                                                                                            tttAGTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCCGTTTGTCGGCCGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGCAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGAttt      >  1:410172/1‑251 (MQ=255)
                                                                                                                                                                                                                                                gTACACAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCCGGGAACCCATACCACAGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCAcc  >  2:351572/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
AACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCCGGTCAGACCCATACAGGAGCGCAGGCCGCCCATCTGCTGGTGAATGATCTCTTTCAGGCGACCTTTATAGGCTACGCGACCTTCGATACCTTCCGGCACCAGTTTGTCGGCAGCGTTATCGCTCTGGAAATAACGGTCAGAGGAACCTTTGGACATCGCGCCCAGGGAACCCATACCACGGTAAGATTTGTAAGAACGGCCCTGGTAGAGTTCGATTTCACC  >  NC_000913/2632467‑2632957

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: