Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,736,589 Δ1 bp 7.2% coding (913/1674 nt) hutU ← urocanate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,736,5890C. 7.2% 259.0 / 15.3 83coding (913/1674 nt)hutUurocanate hydratase
Reads supporting (aligned to +/- strand):  ref base C (40/37);  new base . (4/2);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

GGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGTCCGGGCG  >  NC_002947/5736447‑5736733
                                                                                                                                               |                                                                                                                                                
gggCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATc                                                                                                                                           <  1:1368137/150‑1 (MQ=255)
 ggCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCg                                                                                                                                          >  2:612068/1‑150 (MQ=255)
    ggATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACt                                                                                                                                       >  1:625228/1‑150 (MQ=255)
             ccGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTggcgg                                                                                                                              >  2:1147477/1‑150 (MQ=255)
                 gACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCtt                                                                                                                          <  1:595199/150‑1 (MQ=255)
                      aGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCa                                                                                                                     >  2:180569/1‑150 (MQ=255)
                         cTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTg                                                                                                                  >  1:914094/1‑150 (MQ=255)
                          tGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGc                                                                                                                 >  1:1501925/1‑150 (MQ=255)
                            ggAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCGTTGACCACTGCAg                                                                                                               >  1:169585/1‑150 (MQ=255)
                                  cGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGtt                                                                                                         <  1:1147477/150‑1 (MQ=255)
                                   gAAGGCATTGGCCATGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTc                                                                                                        <  1:1074095/150‑1 (MQ=255)
                                   gAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCGTGCACGTGCAGGCCAtgct                                                                                                                                         <  1:1590551/117‑4 (MQ=255)
                                   gAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCGTGCACGTGCAGGCCAtgct                                                                                                                                         >  2:1590551/1‑114 (MQ=255)
                                   gAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGc                                                                                                              >  1:304976/1‑144 (MQ=255)
                                   gAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGc                                                                                                              <  2:304976/144‑1 (MQ=255)
                                   gAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTc                                                                                                        <  1:244972/150‑1 (MQ=255)
                                   gAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTc                                                                                                        <  1:1624383/150‑1 (MQ=255)
                                    aaGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCTTTGCCTGCACGTGCACGGCAATCGACTGCTGGGCGGCCTTGACCACTGCAGCCGGTTCg                                                                                                       >  1:837660/1‑150 (MQ=255)
                                      ggCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGt                                                                                                     <  1:1209768/150‑1 (MQ=255)
                                         aTTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTg                                                                                                  >  2:160390/1‑150 (MQ=255)
                                          ttGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATc                                                                                                                                           <  1:253236/108‑1 (MQ=255)
                                          ttGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATc                                                                                                                                           >  2:253236/1‑108 (MQ=255)
                                          ttGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCa                                                                                                                     <  1:465065/130‑1 (MQ=255)
                                          ttGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCa                                                                                                                     >  2:465065/1‑130 (MQ=255)
                                              gcACGCCCTCGTGGTTGGCTATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCAc                                                                                             <  2:1501925/149‑1 (MQ=255)
                                               cACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACg                                                                                            >  1:713406/1‑150 (MQ=255)
                                                     ctcctcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATcgcg                                                                                      <  1:503399/150‑1 (MQ=255)
                                                     ctcctcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATcgcg                                                                                      <  1:1054424/150‑1 (MQ=255)
                                                     ctcctcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATcgcg                                                                                      <  1:160390/150‑1 (MQ=255)
                                                      tcctcctTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCgg                                                                                     <  1:1595910/150‑1 (MQ=255)
                                                            ttGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGTGGCCTTGACCACTGCAGCGGGTTCGGTCTGCGCACGATCGCGGTACTGt                                                                               >  2:416583/1‑150 (MQ=255)
                                                            ttGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGt                                                                               >  1:434061/1‑150 (MQ=255)
                                                              ggCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCCGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTc                                                                             <  2:29407/150‑1 (MQ=255)
                                                                     tGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTc                                                                      <  2:713406/150‑1 (MQ=255)
                                                                       gCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTggcg                                                                                                                               <  1:679990/91‑1 (MQ=255)
                                                                       gCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTggcg                                                                                                                               >  2:679990/1‑91 (MQ=255)
                                                                         ggATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGc                                                                  >  1:1039967/1‑150 (MQ=255)
                                                                              ttgttgCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGt                                                                       >  2:1166013/1‑140 (MQ=255)
                                                                              ttgttgCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGt                                                                       <  1:1166013/140‑1 (MQ=255)
                                                                                gttgCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGGACTGTTCCCAGGTCCAGCCGGCTgg                                                           >  2:1442035/1‑150 (MQ=255)
                                                                                   gCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCAGGcgg                                                        <  2:313296/150‑3 (MQ=255)
                                                                                      gTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTa                                                     <  1:612068/150‑1 (MQ=255)
                                                                                       tAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGt                                                                                                     <  1:1556423/101‑1 (MQ=255)
                                                                                       tAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGt                                                                                                     >  2:1556423/1‑101 (MQ=255)
                                                                                       tAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAg                                                    >  2:1524608/1‑150 (MQ=255)
                                                                                       tAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAg                                                    >  1:772292/1‑150 (MQ=255)
                                                                                          tCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTcc                                                                            >  2:210899/1‑123 (MQ=255)
                                                                                          tCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTcc                                                                            <  1:210899/123‑1 (MQ=255)
                                                                                          tCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGc                                                                  >  2:745553/1‑133 (MQ=255)
                                                                                          tCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGc                                                                  <  1:745553/133‑1 (MQ=255)
                                                                                          tCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCg                                                 >  1:254484/1‑150 (MQ=255)
                                                                                            gAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGtt                                               >  2:538024/1‑150 (MQ=255)
                                                                                                gTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGt                                           <  1:1524608/150‑1 (MQ=255)
                                                                                                gTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGt                                           <  1:525909/150‑1 (MQ=255)
                                                                                                 tCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTg                                          >  1:1423373/1‑150 (MQ=255)
                                                                                                       tGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGt                                    <  1:180569/150‑1 (MQ=255)
                                                                                                        gCCCTGCTTCTGGAAGTCCAGCATGGCCtgcacgtgcacg                                                                                                                                                 <  1:1453414/40‑1 (MQ=255)
                                                                                                        gCCCTGCTTCTGGAAGTCCAGCATGGCCtgcacgtgcacg                                                                                                                                                 >  2:1453414/1‑40 (MQ=255)
                                                                                                          gcggCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTggg                                 >  1:1550304/4‑150 (MQ=255)
                                                                                                          ccTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTggg                                 >  2:244018/1‑150 (MQ=255)
                                                                                                           cTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGgc                                >  2:892445/1‑150 (MQ=255)
                                                                                                               ttCTGGAAGTCCCGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGctg                            <  2:914094/150‑1 (MQ=255)
                                                                                                               ttCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCa                                                                    <  2:815429/110‑1 (MQ=255)
                                                                                                               ttCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCa                                                                    >  1:815429/1‑110 (MQ=255)
                                                                                                               ttCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGctg                            <  1:695044/150‑1 (MQ=255)
                                                                                                               atCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGctg                            <  1:1170072/149‑1 (MQ=255)
                                                                                                                 cTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGctggt                          >  2:1356855/1‑150 (MQ=255)
                                                                                                                   ggAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGctg                            <  2:352352/146‑1 (MQ=255)
                                                                                                                   ggAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGctg                            >  1:352352/1‑146 (MQ=255)
                                                                                                                    gAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGt                                           <  2:416406/130‑1 (MQ=255)
                                                                                                                    gAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGt                                           >  1:416406/1‑130 (MQ=255)
                                                                                                                    gAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGgcgc                              <  2:1119471/143‑1 (MQ=255)
                                                                                                                    gAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGgcgc                              >  1:1119471/1‑143 (MQ=255)
                                                                                                                        tCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCg                   >  2:504155/1‑150 (MQ=255)
                                                                                                                         ccAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCgg                  <  2:625228/150‑1 (MQ=255)
                                                                                                                          cAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGt                 <  2:772292/150‑1 (MQ=255)
                                                                                                                          cAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGt                 >  2:251793/1‑150 (MQ=255)
                                                                                                                            gCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGa               >  1:10923/1‑150 (MQ=255)
                                                                                                                            gCAGGGCCTGCACGTGCACGGCCATCGACTGCTTGGTGGCCTTTACCACTGCAGTCGGTTCGGTCTGCGCACGATCCCGGTACTGTTCCCAGGTCCATCCGGCTGGCAGGTAGCCGTTCAGGGGGTCGTGGGTGCTGGTCTGGTCGGTGa               <  1:416583/150‑1 (MQ=255)
                                                                                                                             gaTGGCCGTGCACGTGCAGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGa               >  2:132847/2‑149 (MQ=255)
                                                                                                                             gaTGGCCGTGCACGTGCAGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGa               <  1:132847/148‑1 (MQ=255)
                                                                                                                             gaTGGCCGTGCACGTGCAGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGAc              >  2:940150/2‑150 (MQ=255)
                                                                                                                             gaTGGCCGTGCACGTGCAGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGAc              >  1:312184/2‑150 (MQ=255)
                                                                                                                                  ccTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTCCAGCGGGTCCGGTCTGCGCACGATCGCTGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGt         >  1:1265458/1‑150 (MQ=255)
                                                                                                                                   cTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGTc        >  1:1615396/1‑150 (MQ=255)
                                                                                                                                        acgtgcacgGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGTTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGTCCggg    <  1:766947/150‑1 (MQ=255)
                                                                                                                                          gtgcacgGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGTCCGGGCg  >  2:1072046/1‑150 (MQ=255)
                                                                                                                                             cacgGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAg                                                        >  1:995358/1‑93 (MQ=255)
                                                                                                                                             cacgGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAg                                                        <  2:995358/93‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GGGCGGATGTAGGCCGGGACGAAGCCTGGGAAGTCGAAGGCATTGGCCACGCCCTCCTCCTTGGCCATCTGGCGGATGTTGTTGCCGTAATCGAAGGTCGGGATGCCCTGCTTCTGGAAGTCCAGCATGGCCTGCACGTGCACGGCCATCGACTGCTTGGCGGCCTTGACCACTGCAGCCGGTTCGGTCTGCGCACGATCGCGGTACTGTTCCCAGGTCCAGCCGGCTGGCAGGTAGCCGTTCAGTGGGTCGTGGGCGCTGGTCTGGTCGGTGACCATGTCCGGGCG  >  NC_002947/5736447‑5736733

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